Equation of state of nitrogen (N2) at high pressures and high temperatures: molecular dynamics simulation

Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures above 800 K was studied by molecular dynamics (MD) simulations. The dynamics of the N(2) molecules is treated in hard rotor approximation, i.e., it accounts both translational and rotational degrees of freedom. The rotational motion of the N(2) molecule is treated assuming constant moment of inertia of the nitrogen molecule. The new MD program fully accounts anisotropic molecular nitrogen interaction. The N(2)-N(2) interaction potential has been derived by van der Avoird et al. [J. Chem. Phys. 84, 1629 (1986)] using the results of high precision Hartree-Fock ab initio quantum mechanical calculations. The potential, fully accounts rotational symmetry of the N(2)-N(2) system, by employing 6-j Wigner symbols, i.e., preserving full rotational symmetry of the system. Various numerical algorithms were tested, in order to achieve the energy preservation during the simulation. It has been demonstrated that the standard Verlet algorithm was not preserving the energy for the standard MD time step, equal to 5x10(-16) s. Runge-Kutta fourth order method was able to preserve the energy within 10(-4) relative error, but it requires calculation of the force four times for each time step and therefore it is highly inefficient. A predictor-corrector method of the fifth order (PC5) was found to be efficient and precise and was therefore adopted for the simulation of the molecular nitrogen properties at high pressure. Singer and Fincham algorithms were tested and were found to be as precise as PC5 algorithm and they were also used in the simulation of the equation of state. Results of MD simulations are in very good agreement with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars). For higher pressures, up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 85, 1262 (1984)], the obtained numerical results provide new data of the experimentally unexplored region. These data were formulated in the analytical form of pressure-density-temperature equation of state.     

A new generalization of the Carnahan-Starling equation of state to additive mixtures of hard spheres

We introduce an expansion of the equation of state for additive hard-sphere mixtures in powers of the total packing fraction with coefficients which depend on a set of weighted densities used in scaled particle theory and fundamental measure theory. We demand that the mixture equation of state recovers the quasiexact Carnahan-Starling [J. Chem. Phys. 51, 635 (1969)] result in the case of a one-component fluid and show from thermodynamic considerations and consistency with an exact scaled particle relation that the first and second orders of the expansion lead unambiguously to the Boublik-Mansoori-Carnahan-Starling-Leland [J. Chem. Phys. 53, 471 (1970); J. Chem. Phys. 54, 1523 (1971)] equation and the extended Carnahan-Starling equation introduced by Santos et al. [Mol. Phys. 96, 1 (1999)]. In the third order of the expansion, our approach allows us to define a new equation of state for hard-sphere mixtures which we find to be more accurate than the former equations when compared to available computer simulation data for binary and ternary mixtures. Using the new mixture equation of state, we calculate expressions for the surface tension and excess adsorption of the one-component fluid at a planar hard wall and compare its predictions to available simulation data.     

Recovering the Crooks equation for dynamical systems in the isothermal-isobaric ensemble: a strategy based on the equations of motion

The Crooks equation [Eq. (10) in J. Stat. Phys. 90, 1481 (1998)] relates the work done on a system during a nonequilibrium transformation to the free energy difference between the final and the initial state of the transformation. Recently, the authors have derived the Crooks equation for systems in the canonical ensemble thermostatted by the Nose-Hoover or Nose-Hoover chain method [P. Procacci et al., J. Chem. Phys. 125, 164101 (2006)]. That proof is essentially based on the fluctuation theorem by Evans and Searles [Adv. Phys. 51, 1529 (2002)] and on the equations of motion. Following an analogous approach, the authors derive here the Crooks equation in the context of molecular dynamics simulations of systems in the isothermal-isobaric (NPT) ensemble, whose dynamics is regulated by the Martyna-Tobias-Klein algorithm [J. Chem. Phys. 101, 4177 (1994)]. Their present derivation of the Crooks equation correlates to the demonstration of the Jarzynski identity for NPT systems recently proposed by Cuendet [J. Chem. Phys. 125, 144109 (2006)].     

Influence of high pressure on optical impurity spectra

Expressions are derived for inhomogeneous band shapes of impurity spectra in highly compressed glassy matrices. Intermolecular guest-host interactions are approximated to isotropic two-body Lennard-Jones 6-12 potentials having different parameters in the ground state and the excited state. Calculated shifts and widths are compared to published values for the absorption spectra of phenanthrene, anthracene [B. Y. Okamoto and H. G. Drickamer, J. Chem. Phys. 61, 2870 (1974)], and several polymethine dyes [G. A. Samara et al., J. Chem. Phys. 37, 1482 (1962)] embedded in polymer hosts and subject to pressures up to 140 kbars. The magnitudes of barochromic shifts of the band maxima and the inhomogeneous broadening suggest that the equilibrium coordinates of the excited state are typically less by 5+/-2%.     

Study of thermal properties of the metastable supersaturated vapor with the integral equation method

Pressure, excess chemical potential, and excess free energy data for different densities of the supersaturated argon vapor at reduced temperatures from 0.7 to 1.2 are obtained by solving the integral equation with perturbation correction to the radial distribution function [F. Lado, Phys. Rev. 135, A1013 (1964)]. For those state points where there is no solution, the integral equation is solved with the interaction between argon atoms modeled by Lennard-Jones potential plus a repulsive potential with one controlling parameter, alpha exp(-rsigma) and in the end, all the thermal properties are mapped back to the alpha=0 case. Our pressure data and the spinodal obtained from the current method are compared with a molecular dynamics simulation study [A. Linhart et al., J. Chem. Phys. 122, 144506 (2005)] of the same system.     

Molecular beam optical Stark study of rhodium mononitride

The optical Stark effect in the Q(1) and R(0) lines of the [15.1]1-X (1)Sigma+ (1,0) band of rhodium mononitride (RhN) were recorded and analyzed to determine the permanent electric dipole moments mu for the X (1)Sigma+(upsilon=0) and [15.1]1(upsilon=1) states to be 2.43(5) and 1.75(1) D, respectively. The determined dipole moments are compared to predicted values obtained from density functional theory [Stevens et al., Chem. Phys. Lett. 421, 281 (2006)] and an all-electron ab initio calculation [Shim et al., J. Mol. Struct. THEOCHEM 393, 127 (1997)]. A simple single configuration molecular orbital correlation diagram is used to rationalize the relative values of mu for the 4d mononitrides and RhO. An electronic configuration for the [15.1]1 state is proposed based on the interpretation of the (103)Rh and (14)N magnetic hyperfine interactions.     

Simulation of vibrational dephasing of I(2) in solid Kr using the semiclassical Liouville method

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.     

Molecular dynamics in the isothermal-isobaric ensemble: the requirement of a "shell" molecule. I. Theory and phase-space analysis

Current constant pressure molecular-dynamics (MD) algorithms are not consistent with the recent reformulation of the isothermal-isobaric (NpT) ensemble. The NpT ensemble partition function requires the use of a "shell" molecule to identify uniquely the volume of the system, thereby avoiding the redundant counting of configurations [e.g., G. J. M. Koper and H. Reiss, J. Phys. Chem. 100, 422 (1996); D. S. Corti, Phys. Rev. E, 64, 016128 (2001)]. So far, only the NpT Monte Carlo method has been updated to allow the system volume to be defined by a shell particle [D. S. Corti, Mol. Phys. 100, 1887 (2002)]. A shell particle has yet to be incorporated into MD simulations. The proper modification of the NpT MD algorithm is therefore the subject of this paper. Unlike Andersen's method [H. C. Andersen, J. Chem. Phys. 72, 2384 (1980)] where a piston of unknown mass serves to control the response time of volume fluctuations, the newly proposed equations of motion impose a constant external pressure via the introduction of a shell particle of known mass. Hence, the system itself sets the time scales for pressure and volume fluctuations. The new algorithm is subject to a number of fundamentally rigorous tests to ensure that the equations of motion sample phase space correctly. We also show that the Hoover NpT algorithm [W. G. Hoover, Phys. Rev. A. 31, 1695 (1985); 34, 2499 (1986)] does sample phase correctly, but only when periodic boundary conditions are employed.     

Quantum and quasiclassical studies of the O(3P)+HCl-->OH+Cl(2P) reaction using benchmark potential surfaces

We have performed quantum mechanical (QM) dynamics calculations within the independent-state approximation with new benchmark triplet A" and A' surfaces [B. Ramachandran et al., J. Chem. Phys. 119, 9590 (2003)] for the rovibronic state-to-state measurements of the reaction O(3P)+HCl(v=2,j=1,6,9)-->OH(v'j')+Cl(2P) [Zhang et al., J. Chem. Phys. 94, 2704 (1991)]. The QM and experimental rotational distributions peak at similar OH(j') levels, but the QM distributions are significantly narrower than the measurements and previous quasiclassical dynamics studies. The OH(low j) populations observed in the measurements are nearly absent in the QM results. We have also performed quasiclassical trajectory with histogram binning (QCT-HB) calculations on these same benchmark surfaces. The QCT-HB rotational distributions, which are qualitatively consistent with measurements and classical dynamics studies using other surfaces, are much broader than the QM results. Application of a Gaussian binning correction (QCT-GB) dramatically narrows and shifts the QCT-HB rotational distributions to be in very good agreement with the QM results. The large QCT-GB correction stems from the special shape of the joint distribution of the classical rotational/vibrational action of OH products. We have also performed QM and QCT calculations for the transition, O+HCl(v=0,T=300 K)-->OH(v'j')+Cl from threshold to approximately 130 kcal mol(-1) collision energy as a guide for possible future hyperthermal O-atom measurements. We find in general a mixed energy release into translation and rotation consistent with a late barrier to reaction. Angular distributions at high collision energy are forward peaked, consistent with a stripping mechanism. Direct collisional excitation channel cross sections, O+HCl(v=0,T=300 K)-->O+HCl(v'=1), in the same energy range are large, comparable in magnitude to the reactive channel cross sections. Although the (3)A" state dominates most collision processes, above approximately 48 kcal mol(-1), the (3)A' state plays the major role in collisional excitation.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

Franck-Condon simulation of the single-vibronic-level emission spectra of HPCl/DPCl and the chemiluminescence spectrum of HPCl, including anharmonicity

Restricted-spin coupled-cluster single-double plus perturbative triple excitation [RCCSD(T)] potential energy functions (PEFs) were calculated for the X (2)A" and A (2)A' states of HPCl employing the augmented correlation-consistent polarized-valence-quadruple-zeta (aug-cc-pVQZ) basis set. Further geometry optimization calculations were carried out on both electronic states of HPCl at the RCCSD(T) level with all electron and quasirelativistic effective core potential basis sets of better than the aug-cc-pVQZ quality, and also including some core electrons, in order to obtain more reliable geometrical parameters and relative electronic energy of the two states. Anharmonic vibrational wave functions of the two states of HPCl and DPCl, and Franck-Condon (FC) factors of the A (2)A'-X (2)A" transition were computed employing the RCCSD(T)/aug-cc-pVQZ PEFs. Calculated FC factors with allowance for Duschinsky rotation and anharmonicity were used to simulate the single-vibronic-level (SVL) emission spectra of HPCl and DPCl reported by Brandon et al. [J. Chem. Phys. 119, 2037 (2003)] and the chemiluminescence spectrum reported by Bramwell et al. [Chem. Phys. Lett. 331, 483 (2000)]. Comparison between simulated and observed SVL emission spectra gives the experimentally derived equilibrium geometry of the A (2)A' state of HPCl of r(e)(PCl) = 2.0035 +/- 0.0015 A, theta(e) = 116.08 +/- 0.60 degrees, and r(e)(HP) = 1.4063+/-0.0015 A via the iterative Franck-Condon analysis procedure. Comparison between simulated and observed chemiluminescence spectra confirms that the vibrational population distribution of the A (2)A' state of HPCl is non-Boltzmann, as proposed by Baraille et al. [Chem. Phys. 289, 263 (2003)].     

Complexes of type C6H7(+)·L (L = N2 and CO2) studied by explicitly correlated coupled cluster theory

Complexes of the benzenium ion (C(6)H(7)(+)) with N(2) or CO(2) have been studied by explicitly correlated coupled cluster theory at the CCSD(T)-F12x (x = a, b) level [T. B. Adler et al., J. Chem. Phys. 127, 221106 (2007)] and the double-hybrid density functional B2PLYP-D [T. Schwabe and S. Grimme, Phys. Chem. Chem. Phys. 9, 3397 (2007)]. Improved harmonic vibrational wavenumbers for C(6)H(7)(+) have been obtained by CCSD(T?)-F12a calculations with the VTZ-F12 basis set. Combining them with previous B2PLYP-D anharmonic contributions we arrive at anharmonic wavenumbers which are in excellent agreement with recent experimental data from p-H(2) matrix isolation IR spectroscopy [M. Bahou et al., J. Chem. Phys. 136, 154304 (2012)]. The energetically most favourable conformer of C(6)H(7)(+)·N(2) shows a π-bonded structure similar to C(6)H(7)(+)·Rg (Rg = Ne, Ar) [P. Botschwina and R. Oswald, J. Phys. Chem. A 115, 13664 (2011)] with D(e) ≈ 870 cm(-1). For C(6)H(7)(+)·CO(2), a slightly lower energy is calculated for a conformer with the CO(2) ligand lying in the ring-plane of the C(6)H(7)(+) moiety (D(e) ≈ 1508 cm(-1)). It may be discriminated from other conformers through a strong band predicted at 1218 cm(-1), red-shifted by 21 cm(-1) from the corresponding band of free C(6)H(7)(+).     

Quantum mechanical and quasiclassical trajectory scattering calculations for the C(1D) + H2 reaction on the second excited 1 1A" potential energy surface

Time-independent quantum mechanical (QM) and quasiclassical trajectory (QCT) scattering calculations have been carried out for the C(1D) + H2 --> CH + H reaction at a collision energy of 80 meV on a newly developed ab initio potential energy surface [B. Bussery-Honvault et al., Phys. Chem. Chem. Phys. 7, 1476 (2005)] of 1 1A" symmetry, corresponding to the second singlet state 1 1B1 of CH2. A general good agreement has been found between the QM and QCT rotational distributions and differential cross sections (DCSs). In both cases, DCSs are strongly peaked in the forward direction with a small contribution in the backward direction in contrast with those obtained on the 1 1A' surface, which are nearly symmetric. Rotational distributions obtained on the 1 1A" surface are somewhat colder than those calculated on the 1 1A' surface. The specific dynamics and the contribution of the 1 1A" surface to the overall reactivity of this system are discussed.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Regarding convergence curve of virial expansion for the Lennard-Jones system

We calculate the convergence curve of the virial expansion for a Lennard-Jones system in the density-temperature plane using the approximate method based on the density expansion of the Ornstein-Zernike equation and the condition of thermodynamic consistency [J. Chem. Phys. 106, 6095 (1997)]. At subcritical temperatures, this curve is close to the binodal. At supercritical temperatures, the curve does not coincide with the freezing curve. In the latter case, the densities along the convergence line are distinctly smaller than the densities corresponding to the condensed state.     

Broad shape resonance effects in CaF Rydberg states

Results of ab initio R-matrix calculations [S. N. Altunata et al., J. Chem. Phys. 123, 084319 (2005)] indicate the presence of a broad shape resonance in electron-CaF(+) scattering for the (2)Sigma(+) electronic symmetry near the ionization threshold. The properties of this shape resonance are analyzed using the adiabatic partial-wave expansion of the scattered electron wave function introduced by Le Dourneuf et al. [J. Phys. B 15, L685 (1982)]. The qualitative aspects of the shape resonance are explained by an adiabatic approximation on the electronic motion. Mulliken's rule for the structure of the Rydberg state wave functions [R. S. Mulliken, J. Am. Chem. Soc. 86, 3183 (1964)] specifies that, except for an (n*)(-32) amplitude scale factor, every excited state wave function within one Rydberg series is built on an innermost lobe that remains invariant in shape and nodal position as a function of the excitation energy. Mulliken's rule implies a weak energy dependence of the quantum defects for an unperturbed molecular Rydberg series, which is given by the Rydberg-Ritz formula. This zero-order picture is violated by a single (2)Sigma(+) CaF Rydberg series at all Rydberg state energies (n*=5-->infinity, more so with increasing n*) below the ionization threshold, under the broad width of the shape resonance. Such a violation is diagnostic of a global "scarring" of the Rydberg spectrum, which is distinct from the more familiar local level perturbations.     

Spectroscopy of the UO+2 cation and the delayed ionization of UO2

Vibronically resolved spectra for the UO+2 cation have been recorded using the pulsed field ionization zero electron kinetic energy (PFI-ZEKE) technique. For the ground state, long progressions in both the bending and symmetric stretch vibrations were observed. Bend and stretch progressions of the first electronically excited state were also observed, and the origin was found at an energy of 2678 cm(-1) above the ground state zero-point level. This observation is consistent with a recent theoretical prediction [Infante et al., J. Chem. Phys. 127, 124308 (2007)]. The ionization energy for UO2, derived from the PFI-ZEKE spectrum, namely, 6.127(1) eV, is in excellent agreement with the value obtained from an earlier photoionization efficiency measurement. Delayed ionization of UO2 in the gas phase has been reported previously [Han et al., J. Chem. Phys. 120, 5155 (2004)]. Here, we extend the characterization of the delayed ionization process by performing a quantitative study of the ionization rate as a function of the energy above the ionization threshold. The ionization rate was found to be 5 x 10(6) s(-1) at threshold, and increased linearly with increasing energy in the range investigated (0-1200 cm(-1)).     

A time-dependent wave-packet quantum scattering study of the reaction H2+(v = 0-2,4,6;j = 1) + He--> HeH+ + H

The quantum scattering dynamics calculation was carried out for the titled reaction in the collision energy range of 0.0-2.4 eV with reactant H(2) (+) in the rotational state j = 1 and vibrational states v = 0-2, 4, and 6. The present time-dependent wave-packet calculation takes into account the Coriolis coupling (CC) and uses the accurate ab initio potential-energy surface of Palmieri et al. [Mol. Phys. 98, 1835 (2000)]. The importance of including the CC quantum scattering calculation has been revealed by the comparison between the CC calculation and the previous coupled state (CS) calculation. The CC total cross sections for the v = 2, 4, and 6 states show collision energy-dependent behaviors different from those based on the CS calculation. Furthermore, the collision energy dependence of the total cross sections obtained in the present CC calculation only exhibits minor oscillations, indicating that the chance is slim for reactive resonances in total cross sections to survive through the partial-wave averaging. The magnitude and profile of the CC total cross sections for v = 0-2 in the collision energy range of 0.0-2.5 eV are found to be consistent with experimental cross sections obtained recently by Tang et al. [J. Chem. Phys. 122, 164301 (2005)] after taking into account the experimental uncertainties.     

Accuracy of recent potential energy surfaces for the He-N2 interaction. I. Virial and bulk transport coefficients

A new exchange-Coulomb semiempirical model potential energy surface for the He-N2 interaction has been developed. Together with two recent high-level ab initio potential energy surfaces, it has been tested for the reliability of its predictions of second-virial coefficients and bulk transport phenomena in binary mixtures of He and N2. The agreement with the relevant available measurements is generally within experimental uncertainty for the exchange-Coulomb surface and the ab initio surface of Patel et al. [J. Chem. Phys. 119, 909 (2003)], but with slightly poorer agreement for the earlier ab initio surface of Hu and Thakkar [J. Chem. Phys. 104, 2541 (1996)].     

Surface relaxation in water clusters: evidence from theoretical analysis of the oxygen 1s photoelectron spectrum

We present a theoretical interpretation of the oxygen 1s photoelectron spectrum published by Ohrwall et al. [J. Chem. Phys. 123, 054310 (2005)]. A water cluster that contains 200 molecules was simulated at 215 K using the polarizable AMOEBA force field. The force field predicts longer O...O distances at the cluster surface than in the bulk. Comparisons to ab initio molecular dynamics (MD) simulations indicate that the force field underestimates the degree of surface relaxation. By comparing cluster lineshape models, computed from MD simulations, to the experimental spectrum we find further evidence of surface relaxation.