High Accuracy Potential Energy Curves and Dipole Moment Functions for the Σ and Π Spin States of the CF Diatomic Molecule.

Potential energy curves and dipole moment functions constructed using high-accuracy methods allow for an in-depth examination of the electronic structure of diatomic molecules. computations serve as a valuable complement to experimental data, offering insights into the nature of short-lived molecules such as those encountered within the interstellar medium (ISM). While laboratory experiments provide critical groundwork, the ISM's conditions often permit longer lifetimes for lower stability molecules, enabling unique observations.

High accuracy ab initio potential energy surface for the HO-H van der Waals dimer.

A representation of the three-dimensional potential energy surface (PES) of the HO-H van der Waals dimer is presented. The HO molecule is treated as a rigid body held at its experimentally determined equilibrium geometry, with the OH bond length set to 1.809 650 34 a and the HOH bond angle set to 1.

Selecting Quantum-Chemical Methods for Lanthanide-Containing Molecules: A Balance between Accuracy and Efficiency.

An analysis of how different density functionals, basis sets, and relativistic approximations affect the computed properties of lanthanide-containing molecules allows one to determine which method provides the highest accuracy. Historically, many different density functional methods have been employed to perform calculations on lanthanide complexes and so herein is a detailed analysis of how different methodological combinations change the computed properties of three different families of lanthanide-bearing species: lanthanide diatomic molecules (fluorides and oxides) and their dissociation energies; larger, molecular complexes and their geometries; and lanthanide bis(2-ethylhexyl)phosphate structures and their separation free energies among the lanthanide series. The B3LYP/Sapporo/Douglas-Kroll-Hess (DKH) method was shown to most accurately reproduce dissociation energies calculated at the CCSDT(Q) level of theory with a mean absolute deviation of 1.

Three-body interactions and the elastic constants of hcp solid He.

The effect of three-body interactions on the elastic properties of hexagonal close packed solid He is investigated using variational path integral (VPI) Monte Carlo simulations. The solid's nonzero elastic constants are calculated, at T = 0 K and for a range of molar volumes from 7.88 cm/mol to 20.

Highly Dynamic Anion-Quadrupole Networks in Proteins.

The dynamics of anion-quadrupole (or anion-π) interactions formed between negatively charged (Asp/Glu) and aromatic (Phe) side chains are for the first time computationally characterized in RmlC (Protein Data Bank entry 1EP0 ), a homodimeric epimerase. Empirical force field-based molecular dynamics simulations predict anion-quadrupole pairs and triplets (anion-anion-π and anion-π-π) are formed by the protein during the simulated trajectory, which suggests that the anion-quadrupole interactions may provide a significant contribution to the overall stability of the protein, with an average of -1.6 kcal/mol per pair.

Atomic Spectral Methods for Ab Initio Molecular Electronic Energy Surfaces: Transitioning From Small-Molecule to Biomolecular-Suitable Approaches.

Continuing attention has addressed incorportation of the electronically dynamical attributes of biomolecules in the largely static first-generation molecular-mechanical force fields commonly employed in molecular-dynamics simulations. We describe here a universal quantum-mechanical approach to calculations of the electronic energy surfaces of both small molecules and large aggregates on a common basis which can include such electronic attributes, and which also seems well-suited to adaptation in ab initio molecular-dynamics applications. In contrast to the more familiar orbital-product-based methodologies employed in traditional small-molecule computational quantum chemistry, the present approach is based on an "ex-post-facto" method in which Hamiltonian matrices are evaluated prior to wave function antisymmetrization, implemented here in the support of a Hilbert space of orthonormal products of many-electron atomic spectral eigenstates familiar from the van der Waals theory of long-range interactions.

Search for anisotropy in the Debye-Waller factor of HCP solid (4)He.

The properties of hexagonal close packed (hcp) solid (4)He are dominated by large atomic zero point motions. An accurate description of these motions is therefore necessary in order to accurately calculate the properties of the system, such as the Debye-Waller (DW) factors. A recent neutron scattering experiment reported significant anisotropy in the in-plane and out-of-plane DW factors for hcp solid (4)He at low temperatures, where thermal effects are negligible and only zero-point motions are expected to contribute.

Estimation of the Ground State Energy of an Atomic Solid by Employing Quantum Trajectory Dynamics with Friction.

Evolution with energy dissipation can be used to obtain the ground state of a quantum-mechanical system. This dissipation is introduced in the quantum trajectory framework by adding an empirical friction force to the equations of motion for the trajectories, which, as an ensemble, represent a wave function. The quantum effects in dynamics are incorporated via the quantum force derived from the properties of this ensemble.

Infrared-active spin-orbit transitions of halogen atom dopants in solid parahydrogen: the role of trapping site geometry.

We present theoretical calculations of the (2)P(1/2) ← (2)P(3/2) spin-orbit transition of Cl dopants embedded as substitutional impurities in solid parahydrogen (pH2) matrices. In the lower-energy (2)P(3/2) spin-orbit level, the Cl atom's electron density distribution is anisotropic, and slightly distorts the geometry of the atom's trapping site. This distortion leads to a blue shift in the spin-orbit transition energy; the blue shift is enhanced when we account for the large-amplitude zero point motions of the pH2 molecules surrounding the Cl dopant.

STAAR: statistical analysis of aromatic rings.

The statistical analysis of aromatic rings program allows for an automated search for anion-π interactions between phenylalanine residues and carboxylic acid moieties of neighboring aspartic acid or glutamic acid residues in protein data bank (PDB) structures. The program is written in C++ and is available both as a standalone code and through a web implementation that allows users to upload and analyze biomolecular structures in PDB format. The program outputs lists of Phe/Glu or Phe/Asp pairs involved in potential anion-π interactions, together with geometrical (distance and angle between the Phe's center of mass and Glu or Asp's center of charge) and energetic (quantum mechanical Kitaura-Morokuma interaction energy between the residues) descriptions of each anion-π interaction.

Vibrationally averaged isotropic dispersion energy coefficients of the parahydrogen dimer.

We compare the sum-over-states and coupled cluster linear response formalisms for the determination of imaginary-frequency polarizabilities of H(2). Using both approaches, we compute isotropic dispersion energy coefficients C(n) (n = 6, 8, 10) for H(2)-H(2) molecular pairs over a wide range of H(2) bond lengths. We present vibrationally averaged dispersion energy coefficients for H(2)-H(2), H(2)-D(2), and D(2)-D(2) molecular pairs and examine the coefficients' convergence with respect to basis set.

Pairwise additive model for the He-MgO(100) interaction.

We develop a model, based on pairwise additive He-Mg and He-O interactions, for the potential energy of He adsorbates above a rigid MgO(100) surface. The attractive long-range He-Mg and He-O interactions are assumed to have the form C(6)/r(6), with the C(6) coefficients determined from atomic data within the context of the Slater-Kirkwood approximation. The repulsive short-range He-Mg and He-O interactions are assumed to have the form C(p)/r(p), with the exponent p and the C(p) coefficients taken as adjustable parameters.

A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs.

Protein structures are stabilized using noncovalent interactions. In addition to the traditional noncovalent interactions, newer types of interactions are thought to be present in proteins. One such interaction, an anion-π pair, in which the positively charged edge of an aromatic ring interacts with an anion, forming a favorable anion-quadrupole interaction, has been previously proposed [Jackson, M.

Culture evolves.

Culture pervades human lives and has allowed our species to create niches all around the world and its oceans, in ways quite unlike any other primate. Indeed, our cultural nature appears so distinctive that it is often thought to separate humanity from the rest of nature and the Darwinian forces that shape it. A contrary view arises through the recent discoveries of a diverse range of disciplines, here brought together to illustrate the scope of a burgeoning field of cultural evolution and to facilitate cross-disciplinary fertilization.

An ab initio study of van der Waals potential energy parameters for silver clusters.

We employ ab initio calculations of van der Waals complexes to study the potential energy parameters (C(6) coefficients) of van der Waals interactions for modeling of the adsorption of silver clusters on the graphite surface. Electronic structure calculations of the (Ag(2))(2), Ag(2)-H(2), and Ag(2)-C(6)H(6) complexes are performed using a coupled-cluster approach that includes single, double, and perturbative triple excitations (CCSD(T)), Møller-Plesset second-order perturbation theory (MP2), and spin-component-scaled MP2 (SCS-MP2) methods. Using the atom pair approximation, the C(6) coefficients for silver-silver, silver-hydrogen, and silver-carbon atom systems are obtained after subtracting the energies of quadrupole-quadrupole interactions from the total electronic energy.

Accurate Computation of Electric Field Enhancement Factors for Metallic Nanoparticles Using the Discrete Dipole Approximation.

We model the response of nanoscale Ag prolate spheroids to an external uniform static electric field using simulations based on the discrete dipole approximation, in which the spheroid is represented as a collection of polarizable subunits. We compare the results of simulations that employ subunit polarizabilities derived from the Clausius-Mossotti relation with those of simulations that employ polarizabilities that include a local environmental correction for subunits near the spheroid's surface [Rahmani et al. Opt Lett 27: 2118 (2002)].

A six-dimensional H(2)-H(2) potential energy surface for bound state spectroscopy.

We present a six-dimensional potential energy surface for the (H(2))(2) dimer based on coupled-cluster electronic structure calculations employing large atom-centered Gaussian basis sets and a small set of midbond functions at the dimer's center of mass. The surface is intended to describe accurately the bound and quasibound states of the dimers (H(2))(2), (D(2))(2), and H(2)-D(2) that correlate with H(2) or D(2) monomers in the rovibrational levels (v,j)=(0,0), (0,2), (1,0), and (1,2). We employ a close-coupled approach to compute the energies of these bound and quasibound dimer states using our potential energy surface, and compare the computed energies for infrared and Raman transitions involving these states with experimentally measured transition energies.

A preference for edgewise interactions between aromatic rings and carboxylate anions: the biological relevance of anion-quadrupole interactions.

Noncovalent interactions are quite important in biological structure-function relationships. To study the pairwise interaction of aromatic amino acids (phenylalanine, tyrosine, tryptophan) with anionic amino acids (aspartic and glutamic acids), small molecule mimics (benzene, phenol or indole interacting with formate) were used at the MP2 level of theory. The overall energy associated with an anion-quadrupole interaction is substantial (-9.

Quasiclassical trajectory study of the vibrational quenching of hydroxyl radicals through collision with O atoms.

The collisional removal of vibrationally excited OH radicals by O atoms is studied by the quasiclassical trajectory method. To evaluate the effect of different topological features on the scattering processes two different global potential energy surfaces, DMBE IV and TU, are used. Results for reactive, exchange, and inelastic scattering probabilities are reported for central collisions (with zero total angular momentum) with a fixed relative translational energy for vibrational levels of OH ranging from nu=1 to v=8.

Population size bias in descendant-weighted diffusion quantum Monte Carlo simulations.

We consider the influence of population size on the accuracy of diffusion quantum Monte Carlo simulations that employ descendant weighting or forward walking techniques to compute expectation values of observables that do not commute with the Hamiltonian. We show that for a simple model system, the d-dimensional isotropic harmonic oscillator, the population size must increase rapidly with d in order to ensure that the simulations produce accurate results. When the population size is too small, expectation values computed using descendant-weighted diffusion quantum Monte Carlo simulations exhibit significant systematic biases.

Interaction-induced dipole moment of the Ar-H2 dimer: dependence on the H2 bond length.

We present ab initio calculations of the interaction-induced dipole moment of the Ar-H2 van der Waals dimer. The primary focus of our calculations is on the H2 bond length dependence of the dipole moment, which determines the intensities of both the collision-induced H2 upsilon = 1 <-- 0 fundamental band in gaseous Ar-H2 mixtures and the dopant-induced H2 upsilon = 1 <-- 0 absorption feature in Ar-doped solid H2 matrices. Our calculations employ large atom-centered basis sets, diffuse bond functions positioned between the two monomers, and a coupled cluster treatment of valence electron correlation; core-valence correlation effects appear to make negligible contributions to the interaction-induced dipole moment for the Ar-H2 configurations considered here.

Vibrational dependence of the H2-H2 C6 dispersion coefficients.

We use the sum-over-states formalism to compute the imaginary-frequency dipole polarizabilities for H2, as a function of the H-H bond length, at the full configuration interaction level of theory using atom-centered d-aug-cc-pVQZ basis sets. From these polarizabilities, we obtain isotropic and anisotropic C6 dispersion coefficients for a pair of H2 molecules as functions of the two molecules' bond lengths.