A new intermolecular potential energy surface for carbon dioxide from ab initio calculations

Quantum modeling, beyond secularity, of the collisional dissipation of molecular alignment using the energy-corrected sudden approximation

We propose a Markovian quantum model for the time dependence of the pressure-induced decoherence of rotational wave packets of gas-phase molecules beyond the secular approximation. It is based on a collisional relaxation matrix constructed using the energy-corrected sudden approximation, which improves the previously proposed infinite order sudden one by taking the molecule rotation during collisions into account. The model is tested by comparisons with time-domain measurements of the pressure-induced decays of molecular-axis alignment features (revivals and echoes) for HCl and CO2 gases, pure and diluted in He. For the Markovian systems HCl-He and CO2-He, the comparisons between computed and measured data demonstrate the robustness of our approach, even when the secular approximation largely breaks down. In contrast, significant differences are obtained in the cases of pure HCl and CO2, for which the model underestimates the decay rate of the alignment at short times. This result is attributed to the non-Markovianity of HCl-HCl and CO2-CO2 interactions and the important contribution of those collisions that are ongoing at the time when the system is excited by the aligning laser pulse.

Effect of the range of particle cohesion on the phase behavior and thermodynamic properties of fluids

Molecular simulations are performed for the (m + 1, m) potential to systematically investigate the effect of changing the range of particle cohesion on both vapor-liquid equilibria and thermodynamic properties of fluids. The results are reported for m = 4-11, which represent a progressive narrowing of the potential energy well. The conventional Lennard-Jones potential is used as a reference point for normal fluid behavior. Small values of m result in a broadening of the phase envelope compared with the Lennard-Jones potential, whereas a contraction is observed in other cases. The critical properties are reported, and a relationship between the critical temperature and the Boyle temperature is determined. The low values of the critical compressibility factor when m < 6 reflect the behavior observed for real fluids such as n-alkanes. The results for supercritical thermodynamic properties are much more varied. Properties such as pressure, potential energy, isochoric thermal pressure coefficient, and thermal expansion coefficient vary consistently with m, whereas other properties such as the Joule-Thomson coefficient exhibit much more nuanced behavior. Maximum and minimum values are reported for both the isochoric heat capacity and isothermal compressibility. A minimum in the speed of sound is also observed.

Ab Initio Intermolecular Potential Energy Surfaces and Cross Second Virial Coefficients for the Dimer N2-NO

Ab initio intermolecular potential energy surfaces (PES) of N2-NO have been constructed at the level of theory CCSD(T) with the augmented correlation-consistent basis sets aug-cc-pVmZ (with m = 2, 3, 4). The nitrogen in the closed-shell electronic configuration X1Σ+ and nitric oxide in the open-shell electronic configuration A2Σ+ were employed to calculate ab initio intermolecular interaction energies. The two new ab initio 5-site intermolecular pair potentials at the theoretical level CCSD(T)/aug-cc-pVmZ (with m = 4, 24) were developed appropriately and are suitable for N2-NO dimer by using the nonlinear least-squares fitting method combining MIO and Levenberg-Marquardt algorithms. The correlation quality of these two potentials was found to be very good with R 2 values in the range of 0.98372 to 0.99775. The cross second virial coefficients B 12(T) of the N2-NO dimer were calculated in the temperature range of 100 to 470 K using the two ab initio 5-site potentials. The discrepancies between the calculated results and the experimental data can be acceptable.

Diffusion in Binary Aqueous Solutions of Alcohols by Molecular Simulation

Based on the molecular dynamics method, the calculations for diffusion coefficients were carried out in binary aqueous solutions of three alcohols: ethanol, isopropanol, and tert-butanol. The intermolecular potential TIP4P/2005 was used for water; and five force fields were analyzed for the alcohols. The force fields providing the best accuracy of calculation were identified based on a comparison of the calculated self-diffusion coefficients of pure alcohols with the experimental data for internal (Einstein) diffusion coefficients of alcohols in solutions. The temperature and concentration dependences of the interdiffusion coefficients were determined using Darken’s Equation. Transport (Fickian) diffusion coefficients were calculated using a thermodynamic factor determined by the non-random two-liquid (NRTL) and Willson models. It was demonstrated that for adequate reproduction of the experimental data when calculating the transport diffusion coefficients, the thermodynamic factor has to be 0.64. Simple approximations were obtained, providing satisfactory accuracy in calculating the concentration and temperature dependences of the transport diffusion coefficients in the studied mixtures.

Lennard-Jones Intermolecular Potentials for the Description of 6-Membered Aromatic Heterocycles Interacting with the Isoelectronic CO2 and CS2

We have generated Lennard-Jones potentials for the interaction between CX₂ (X = O, S) and 11 nitrogen-doped benzene derivatives in different orientations at the M06-2X/def2-tzvpp level as tools to parametrize accurate force fields and to better understand the interaction of these greenhouse gases with heterocyclic building blocks used in the design of capture and detection systems. We find that the most favorable interactions are found between the carbon in CO₂ and the main heterocycle in the ring in a parallel orientation, whereas the preferred interaction mode of CS₂ is established between sulfur and the π density of the aromatic ring. The fact that the preferences for interaction sites and orientations of CO₂ and CS₂ are most of the times opposite helps in terms of ensuring the selectivity of these systems in front of these two isoelectronic compounds. The existence of very good linear correlations ( R² values very close to one) between the number of nitrogen atoms in the heterocyclic ring and the depth of the interaction potential wells opens the door to the use of these results in generating coarse-grained potentials or models with predictive power for use in the design of larger systems.

Rotational Echoes as a Tool for Investigating Ultrafast Collisional Dynamics of Molecules

We show that recently discovered rotational echoes of molecules provide an efficient tool for studying collisional molecular dynamics in high-pressure gases. Our study demonstrates that rotational echoes enable the observation of extremely fast collisional dissipation, at timescales of the order of a few picoseconds, and possibly shorter. The decay of the rotational alignment echoes in CO_{2} gas and CO_{2}-He mixture up to 50 bar was studied experimentally, delivering collision rates that are in good agreement with the theoretical expectations. The suggested measurement protocol may be used in other high-density media, and potentially in liquids.

Thermodynamic Properties of Vapor-Liquid Equilibria from Monte-Carlo Simulation using ab initio Intermolecular Potentials of Systems H2-H2 and F2-F2

In this work, we have been carried out GEMC-NVT simulations in the temperature range 18 K–32 K for fluid hydrogen and in range 60 K–140 K for fluid fluorine using four our developed ab initio 5-site intermolecular potentials for dimers H2-H2 and F2-F2, respectively. The thermodynamic properties of vapor-liquid equilibria and the critical points of fluids hydrogen and fluorine were calculated with the obtained densities of coexisting phases and vapor pressures. The simulation results drived from ab initio pair potentials were compared with those from ab initio potential plus three-body Axilrod-Teller potential and experimental data as well as those from Monte Carlo simulation using Lennard-Jones potentials, Deiters equation of state (D1-EOS) and Benedict-Webb-Rubin equation of state (EOS) reported in the literature.

Simulations of hydrogen, carbon dioxide, and small hydrocarbon sorption in a nitrogen-rich rht-metal–organic framework

Detailed theoretical insights into the gas-sorption mechanism of Cu-TDPAH are presented for the first time.

Five intermolecular vibrations of the CO2 dimer observed via infrared combination bands

The weakly bound van der Waals dimer (CO₂)₂ has long been of considerable theoretical and experimental interest. Here, we study its low frequency intermolecular vibrations by means of combination bands in the region of the CO₂ monomer ν₃ fundamental (≈2350 cm^-1), which are observed using a tunable infrared laser to probe a pulsed supersonic slit jet expansion. With the help of a recent high level ab initio calculation by Wang, Carrington, and Dawes, four intermolecular frequencies are assigned: the in-plane disrotatory bend (22.26 cm^-1); the out-of-plane torsion (23.24 cm^-1); twice the disrotatory bend (31.51 cm^-1); and the in-plane conrotatory bend (92.25 cm^-1). The disrotatory bend and torsion, separated by only 0.98 cm^-1, are strongly mixed by Coriolis interactions. The disrotatory bend overtone is well behaved, but the conrotatory bend is highly perturbed and could not be well fitted. The latter perturbations could be due to tunneling effects, which have not previously been observed experimentally for CO₂ dimer. A fifth combination band, located 1.3 cm^-1 below the conrotatory bend, remains unassigned.

Transferable next-generation force fields from simple liquids to complex materials

Molecular simulations have had a transformative impact on chemists' understanding of the structure and dynamics of molecular systems. Simulations can both explain and predict chemical phenomena, and they provide a unique bridge between the microscopic and macroscopic regimes. The input for such simulations is the intermolecular interactions, which then determine the forces on the constituent atoms and therefore the time evolution and equilibrium properties of the system. However, in practice, accuracy and reliability are often limited by the fidelity of the description of those very same interactions, most typically embodied approximately in mathematical form in what are known as force fields. Force fields most often utilize conceptually simple functional forms that have been parametrized to reproduce existing experimental gas phase or bulk data. Yet, reliance on empirical parametrization can sometimes introduce limitations with respect to novel chemical systems or uncontrolled errors when moving to temperatures, pressures, or environments that differ from those for which they were developed. Alternatively, it is possible to develop force fields entirely from first principles, using accurate electronic structure calculations to determine the intermolecular interactions. This introduces a new set of challenges, including the transferability of the resulting force field to related chemical systems. In response, we recently developed an alternative approach to develop force fields entirely from first-principles electronic structure calculations based on intermolecular perturbation theory. Making use of an energy decomposition analysis ensures, by construction, that the resulting force fields contain the correct balance of the various components of intermolecular interaction (exchange repulsion, electrostatics, induction, and dispersion), each treated by a functional form that reflects the underlying physics. We therefore refer to the resulting force fields as physically motivated. We find that these physically motivated force fields exhibit both high accuracy and transferability, with the latter deriving from the universality of the fundamental physical laws governing intermolecular interactions. This basic methodology has been applied to a diverse set of systems, ranging from simple liquids to nanoporous metal-organic framework materials. A key conclusion is that, in many cases, it is feasible to account for nearly all of the relevant physics of intermolecular interactions within the context of the force field. In such cases, the structural, thermodynamic, and dynamic properties of the system become naturally emergent, even in the absence of explicit parameterization to bulk properties. We also find that, quite generally, the three-body contributions to the dispersion and exchange energies in bulk liquids are crucial for quantitative accuracy in a first-principles force field, although these contributions are almost universally neglected in existing empirical force fields.

CO2 isolated line shapes by classical molecular dynamics simulations: influence of the intermolecular potential and comparison with new measurements

Room temperature absorption spectra of various transitions of pure CO2 have been measured in a broad pressure range using a tunable diode-laser and a cavity ring-down spectrometer, respectively, in the 1.6 μm and 0.8 μm regions. Their spectral shapes have been calculated by requantized classical molecular dynamics simulations. From the time-dependent auto-correlation function of the molecular dipole, including Doppler and collisional effects, spectral shapes are directly computed without the use of any adjusted parameter. Analysis of the spectra calculated using three different anisotropic intermolecular potentials shows that the shapes of pure CO2 lines, in terms of both the Lorentz widths and non-Voigt effects, slightly depend on the used potential. Comparisons between these ab initio calculations and the measured spectra show satisfactory agreement for all considered transitions (from J = 6 to J = 46). They also show that non-Voigt effects on the shape of CO2 transitions are almost independent of the rotational quantum number of the considered lines.

Molecular dynamics simulations for CO2 spectra. IV. Collisional line-mixing in infrared and Raman bands

Ab initio calculations of the shapes of pure CO2 infrared and Raman bands under (pressure) conditions for which line-mixing effects are important have been performed using requantized classical molecular dynamics simulations. This approach provides the autocorrelation functions of the dipole vector and isotropic polarizability whose Fourier-Laplace transforms yield the corresponding spectra. For that, the classical equations of dynamics are solved for each molecule among several millions treated as linear rigid rotors and interacting through an anisotropic intermolecular potential. Two of the approximations used in the previous studies have been corrected, allowing the consideration of line-mixing effects without use of any adjusted parameters. The comparisons between calculated and experimental spectra under various conditions of pressure and temperature demonstrate the quality of the theoretical model. This opens promising perspectives for first principle ab initio predictions of line-mixing effects in absorption and scattering spectra of various systems involving linear molecules.

A molecular H2 potential for heterogeneous simulations including polarization and many-body van der Waals interactions

A highly accurate aniostropic intermolecular potential for diatomic hydrogen has been developed that is transferable for molecular modeling in heterogeneous systems. The potential surface is designed to be efficacious in modeling mixed sorbates in metal-organic materials that include sorption interactions with charged interfaces and open metal sites. The potential parameters are compatible for mixed simulations but still maintain high accuracy while deriving dispersion parameters from a proven polarizability model. The potential includes essential physical interactions including: short-range repulsions, dispersion, and permanent and induced electrostatics. Many-body polarization is introduced via a point-atomic polarizability model that is also extended to account for many-body van der Waals interactions in a consistent fashion. Permanent electrostatics are incorporated using point partial charges on atomic sites. However, contrary to expectation, the best potentials are obtained by permitting the charges to take on values that do not reproduce the first non-vanishing moment of the electrostatic potential surface, i.e., the quadrupole moment. Potential parameters are fit to match ab initio energies for a representative range of dimer geometries. The resulting potential is shown to be highly effective by comparing to electronic structure calculations for a thermal distribution of trimer geometries, and by reproducing experimental bulk pressure-density isotherms. The surface is shown to be superior to other similarly portable potential choices even in tests on homogeneous systems without strong polarizing fields. The present streamlined approach to developing such potentials allows for a simple adaptation to other molecules amenable to investigation by high-level electronic structure methods.

High resolution infrared spectroscopy of carbon dioxide clusters up to (CO2)13

Thirteen specific infrared bands in the 2350 cm(-1) region are assigned to carbon dioxide clusters, (CO(2))(N), with N = 6, 7, 9, 10, 11, 12 and 13. The spectra are observed in direct absorption using a tuneable infrared laser to probe a pulsed supersonic jet expansion of a dilute mixture of CO(2) in He carrier gas. Assignments are aided by cluster structure calculations made using two reliable CO(2) intermolecular potential functions. For (CO(2))(6), two highly symmetric isomers are observed, one with S(6) symmetry (probably the more stable form), and the other with S(4) symmetry. (CO(2))(13) is also symmetric (S(6)), but the remaining clusters are asymmetric tops with no symmetry elements. The observed rotational constants tend to be slightly (≈2%) smaller than those from the predicted structures. The bands have increasing vibrational blueshifts with increasing cluster size, similar to those predicted by the resonant dipole-dipole interaction model but significantly larger in magnitude.

Molecular dynamics simulations for CO2 spectra. III. Permanent and collision-induced tensors contributions to light absorption and scattering

Classical molecular dynamics simulations have been performed for gaseous CO(2) starting from an accurate anisotropic intermolecular potential. Through calculations of the evolutions of the positions and orientations of a large number of molecules, the time evolutions of the permanent and collision-induced electric dipole vector and polarizability tensor are obtained. These are computed from knowledge of static molecular parameters taking only the leading induction terms into account. The Laplace transforms of the auto-correlation functions of these tensors then directly yield the light absorption and scattering spectra. These predictions are, to our knowledge, the first in which the contributions of permanent and collision-induced tensors are simultaneously taken into account for gaseous CO(2), without any adjusted parameter. Comparisons between computations and measurements are made for absorption in the region of the ν(3) infrared band and for depolarized Rayleigh scattering in the roto-translational band. They demonstrate the quality of the model over spectral ranges from the band center to the far wings where the spectrum varies by several orders of magnitude. The contributions of the permanent and interaction-induced (dipole and polarizability) tensors are analyzed for the first time, through the purely permanent (allowed), purely induced, and cross permanent∕induced components of the spectra. It is shown that, while the purely induced contribution is negligible when compared to the collision-broadened allowed component, the cross term due to interferences between permanent and induced tensors significantly participates to the wings of the bands. This successfully clarifies the long lasting, confusing situation for the mechanisms governing the wings of the CO(2) spectra considered in this work.