Pairwise and many-body contributions to interaction potentials in He(n) clusters

High level ab initio calculations have been carried out to assess the pairwise additivity of potentials in the attractive or well regions of the potential surfaces of clusters of helium atoms. A large basis set was employed and calculations were done at the Brueckner orbital coupled cluster level. Differences between calculated potentials for several interacting atoms and the corresponding summed pair potentials reveal the three‐body and certain higher order contributions to the interaction strengths. Attraction between rare gas atoms develops from dispersion, and so helium clusters provide the most workable systems for analyzing nonadditivity of dispersion. The results indicate that the many‐body or nonpairwise contributions tend to be less than a few percent of the attractive interaction across regions around the minima of the potential energy surfaces of small clusters. Dipole–dipole–dipole dispersion and dipole–dipole–quadrupole dispersion are noticeable parts of the small three‐body terms

Three-body analytical potential for interacting helium atoms

Large basis set ab initio calculations have been carried out for a dense grid of points on the He, potential energy surface. Three-body contributions were extracted at every point, and a number of concise functional representations for the three-body potential surface were then examined. Three-body multipolar dispersion terms and other radial and angular terms were used in the representations, and an assessment of relative importance of the different terms is presented. Combined with a two-body He-He potential, the results of this work should offer a high quality interaction potential for simulations of aggregated helium

An Extended Multireference Study of the Electronic States of Para-benzyne

A state-averaged, multireference complete active space (CAS) approach was used for the determination of the vertical excitation energies of valence and Rydberg states of para-benzyne. Orbitals were generated with a 10- and 32-state averaged multiconfigurational self-consistent field approach. Electron correlation was included using multireference configuration interaction with singles and doubles, including the Pople correction for size extensivity, multireference averaged quadratic coupled cluster (MR-AQCC) and MR-AQCC based on linear response theory. There is a very high density of electronic states in this diradical system—there are more than 17 states within 7 eV of the ground state including two 3s Rydberg states. All excitations, except 2 1Ag. are from the π system to the δδ*system. Of the 32 states characterized, 15 were multiconfigurational, including the ground 1Ag state, providing further evidence for the necessity of a multireference approach for p-benzyne. The vertical singlet-triplet splitting was also characterized using a two-state averaged approach. A CAS (2,2) calculation was shown to be inadequate due to interaction with the π orbitals

Molecular Dynamics Study of the Opening Mechanism for DNA Polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although Xray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme:DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics

A comparison of the Low Mode and Monte Carlo conformational search methods

The Low Mode (LM) and Monte Carlo (MC) conformational search methods were compared on three diverse molecular systems; (4R, 5S, 6S, 7R)-hexahydro-5,6-dihydroxy-1,3,4,7-tetrakis(phenylmethyl)-2H-1,3-diazapin-2-one (1), 2-methoxy-2-phenyl-2-triflouromethyl-N-α-methyl benzyl propanamide (2) and a trimeric 39-membered polyazamacrolide (3). We find that either method, or a combination of the methods, is equally efficient at searching the conformational space of the smaller molecular systems while a 50:50 hybrid of Low Mode and Monte Carlo is most efficient at searching the space of the larger molecular system

Halogen Bonding Interactions for Aromatic and Non-Aromatic Explosive Detection

Improved sensing strategies are needed for facile, accurate and rapid detection of aromatic and nonaromatic explosives. Density functional theory was used to evaluate the relative binding interaction energies between halogen-containing sensor model molecules and nitro-containing explosives. Interaction energies ranged from –18 to –14 kJ/mol and highly directional halogen bonding interactions were observed with bond distances ranging between 3.0 and 3.4 Å. In all geometry optimized structures, the sigma-hole of electropositive potential on the halogen aligned with a lone pair of electrons on the nitro-moiety of the explosive. The computational results predict that the strongest interactions will occur with iodine-based sensors as, of all the halogens studied, iodine is the largest, most polarizable halogen with the smallest electronegativity. Based on these promising proof-of-concept results, synthetically accessible sensors were designed using1, 4-dihalobenzene (X= Cl, Br and I) with and without tetra-fluoro electron withdrawing groups attached to the benzene ring. These sensing molecules were embedded onto single walled carbon nanotubes that were mechanically abraded onto interdigitated array electrodes and these were used to measure the responses to explosive model compounds cyclohexanone and dimethyl-dinitro-benzene in nitrogen gas. Amperometric current-time curves for selectors and control molecules, including concentration correlated signal enhancement, as well as response and recovery times, indicate selector responsiveness to these model compounds, with the largest response observed for iodo-substituted sensors

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calcu- lations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of dia- batization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully vari- ational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on stan- dard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview

Challenging Disciplinary Boundaries in the First Year: A New Introductory Integrated Science Course for STEM Majors

To help undergraduates make connections among disciplines so they are able to approach, evaluate, and contribute to the solutions of important global problems, our campus has been focused on interdisciplinary research and education opportunities across the science, technology, engineering, and mathematics (STEM) disciplines. This paper describes the mobilization, planning, and implementation of a first-year interdisciplinary course for STEM majors that integrates key concepts found in traditional first-semester biology, chemistry, computer science, mathematics, and physics courses. This team-taught course, Integrated Quantitative Science (IQS), is half of a first-year student’s schedule in both semesters and is composed of a double lecture and a weekly lab and workshop. Features of this first-year course, including themes and concepts covered each semester, some of the materials developed, lessons learned, challenges, and preliminary measures of success are described

Challenging Disciplinary Boundaries in the First Year: A New Introductory Integrated Science Course for STEM Majors

To help undergraduates make connections among disciplines so they are able to approach, evaluate, and contribute to the solutions of important global problems, our campus has been focused on interdisciplinary research and education opportunities across the science, technology, engineering, and mathematics (STEM) disciplines. This paper describes the mobilization, planning, and implementation of a first-year interdisciplinary course for STEM majors that integrates key concepts found in traditional first-semester biology, chemistry, computer science, mathematics, and physics courses. This team-taught course, Integrated Quantitative Science (IQS), is half of a first-year student’s schedule in both semesters and is composed of a double lecture and a weekly lab and workshop. Features of this first-year course, including themes and concepts covered each semester, some of the materials developed, lessons learned, challenges, and preliminary measures of success are described