Search

In two recent papers [A. D. Becke, J. Chem. Phys. 156, 214101 (2022) and A. D. Becke, J. Chem. Phys. 157, 234102 (2022)], we compared two Kohn–Sham density functionals based on physical modeling and theory with the best density-functional power-series fits in the literature. The best error statistics reported to date for a hybrid functional on the general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)] were obtained. In the present work, additional second-order perturbation-theory terms are considered. The result is a 12-parameter double-hybrid density functional with the lowest GMTKN55 WTMAD2 "weighted total mean absolute deviation" error (1.76 kcal/mol) yet seen for any hybrid or double-hybrid density-functional approximation. We call it "DH23." [ABSTRACT FROM AUTHOR]

In a previous paper [B. Poirier, J. Chem. Phys. 121, 4501 (2004)] a unique bipolar decomposition Ψ=Ψ1+Ψ2 was presented for stationary bound states Ψ of the one-dimensional Schrödinger equation, such that the components Ψ1 and Ψ2 approach their semiclassical WKB analogs in the large-action limit. The corresponding bipolar quantum trajectories, as defined in the usual Bohmian mechanical formulation, are classical-like and well behaved, even when Ψ has many nodes or is wildly oscillatory. A modification for discontinuous potential stationary scattering states was presented in a second, companion paper [C. Trahan and B. Poirier, J. Chem. Phys.124, 034115 (2006), previous paper], whose generalization for continuous potentials is given here. The result is an exact quantum scattering methodology using classical trajectories. For additional convenience in handling the tunneling case, a constant-velocity-trajectory version is also developed. [ABSTRACT FROM AUTHOR]

Light-induced orientation of gas phase molecules is a long-pursued goal in physics and chemistry. Here, we experimentally demonstrate a six-fold increase in the terahertz-induced orientation of iodomethane (CH3I) molecules at room temperature, provided by rotational pre-excitation with a moderately intense near-IR pulse. The paper highlights the underlying interference of multiple coherent transition pathways within the rotational coherence manifold and is analyzed accordingly. Our experimental and theoretical results provide desirable and practical means for all-optical experiments on oriented molecular ensembles. [ABSTRACT FROM AUTHOR]

In a previous paper [B. Poirier, J. Chem. Phys. 121, 4501 (2004)] a unique bipolar decomposition, Ψ=Ψ1+Ψ2, was presented for stationary bound states Ψ of the one-dimensional Schrödinger equation, such that the components Ψ1 and Ψ2 approach their semiclassical WKB analogs in the large action limit. Moreover, by applying the Madelung-Bohm ansatz to the components rather than to Ψ itself, the resultant bipolar Bohmian mechanical formulation satisfies the correspondence principle. As a result, the bipolar quantum trajectories are classical-like and well behaved, even when Ψ has many nodes or is wildly oscillatory. In this paper, the previous decomposition scheme is modified in order to achieve the same desirable properties for stationary scattering states. Discontinuous potential systems are considered (hard wall, step potential, and square barrier/well), for which the bipolar quantum potential is found to be zero everywhere, except at the discontinuities. This approach leads to an exact numerical method for computing stationary scattering states of any desired boundary conditions, and reflection and transmission probabilities. The continuous potential case will be considered in a companion paper [C. Trahan and B. Poirier, J. Chem. Phys. 124, 034116 (2006), following paper]. [ABSTRACT FROM AUTHOR]

The theory of adiabatic electron transport in a correlated condensed-matter system is rooted in a seminal paper by Niu and Thouless [J. Phys. A: Math. Gen. 17, 2453 (1984)]; I adopt here an analogous logic in order to retrieve the known expression for the adiabatic electronic flux in a molecular system [L. A. Nafie, J. Chem. Phys. 79, 4950 (1983)]. Its derivation here is considerably simpler than those available in the current quantum-chemistry literature; it also explicitly identifies the adiabaticity parameter, in terms of which the adiabatic flux and the electron density are both exact to first order. It is shown that the continuity equation is conserved to the same order. For the sake of completeness, I also briefly outline the relevance of the macroscopic electronic flux to the physics of solids and liquids. [ABSTRACT FROM AUTHOR]

In this paper, we introduce new integrators of Newton-Cotes type and investigate the connection between these new methods, differential methods, and symplectic integrators. From the literature, we can see that several one step symplectic integrators have been obtained based on symplectic geometry. However, the investigation of multistep symplectic integrators is very poor. In this paper, we introduce a new numerical method of closed Newton-Cotes type and we write it as a symplectic multilayer structure. We apply the symplectic schemes in order to solve Hamilton's equations of motion which are linear in position and momentum. We observe that the Hamiltonian energy of the system remains almost constant as integration proceeds. [ABSTRACT FROM AUTHOR]

In a previous paper [J. Theo. Comput. Chem. 2, 65 (2003)], one of the authors (B.P.) presented a method for solving the multidimensional Schrödinger equation, using modified Wilson-Daubechies wavelets, and a simple phase space truncation scheme. Unprecedented numerical efficiency was achieved, enabling a ten-dimensional calculation of nearly 600 eigenvalues to be performed using direct matrix diagonalization techniques. In a second paper [J. Chem. Phys. 121, 1690 (2004)], and in this paper, we extend and elaborate upon the previous work in several important ways. The second paper focuses on construction and optimization of the wavelength functions, from theoretical and numerical viewpoints, and also examines their localization. This paper deals with their use in representations and eigenproblem calculations, which are extended to 15-dimensional systems. Even higher dimensionalities are possible using more sophisticated linear algebra techniques. This approach is ideally suited to rovibrational spectroscopy applications, but can be used in any context where differential equations are involved. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We apply the Ornstein-Zernike equations together with the extended Zerah-Hansen approach, described in a previous paper [J. Chem. Phys. 109, 11074 (1998)], to determine the microstructure of a colloidal suspension at finite concentrations. Using these results as an input, the effective pair potential between two colloidal particles is also calculated following the general approach described in that paper. It is found that, for sufficiently large concentrations, this effective potential develops an attractive well with a minimum located at a distance of almost four macroparticle diameters. It is also shown that this attraction is generated by the inversion of the sign of the effective charged distribution surrounding each macroparticle involved in the electrostatic component of the effective pair potential. [ABSTRACT FROM AUTHOR]

It is well documented that conventional methods such as dynamic light scattering have encountered difficulties in characterizing charged macromolecules and, therefore, it is desirable that new methods and techniques are introduced. With the ultra-high sensitivity, single molecule fluorescence spectroscopy has successfully lowered the detection limit considerably and enabled measurement under extreme dilution conditions—around the concentration of 10−9M—at which the effect of inter-chain electrostatic repulsion is suppressed. Furthermore, the excellent spatial and temporal resolution as well as the capacity of molecular recognition of these methods help in obtaining rich information of charged macromolecules. This paper summarizes the applications of single molecule fluorescence spectroscopy, especially fluorescence correlation spectroscopy and photon counting histogram, in the studies on charged macromolecules in aqueous solutions and plenty of new information has been revealed on the molecular conformation, counterion distribution, and a few important governing factors. The powerfulness and effectiveness of single molecule fluorescence spectroscopy make it promising in the investigations of charged macromolecules. [ABSTRACT FROM AUTHOR]

In this paper, we present a new approach to study the probability distribution of the radius of gyration of freely jointed chains. In this approach, we study the joint distribution PN(Rg2,Rvector G), where Rg is the radius of gyration and Rvector G is the vector from the first bead to the center of the chain. We derive an iteration relation between the probability distributions of chains with lengths N and N+1. When the chain length is large, the iteration relation can be approximated by an evolution equation with integral determining condition. The equation is solved to obtain the probability distribution function, which yields the known result obtained by Fixman [J. Chem. Phys. 36, 306 (1962)] using the integral method. We hope the iteration method in this paper could be used to study the probability distribution of the structure factor of polymer molecules. [ABSTRACT FROM AUTHOR]

Nematic liquid crystals confined in a cylindrical capillary and subjected to strong homeotropic anchoring conditions is a long-studied fundamental problem that uniquely incorporates nonlinearity, topological stability, defects, and texture physics. The observed and predicted textures that continue to be investigated include escape radial, radial with a line defect, planar polar with two line defects, and periodic array of point defects. This paper presents theory and multiscale simulations of global and fine scale textures of nematic point defects, based on the Landau–de Gennes tensor order parameter equations. The aim of this paper is to further investigate the ringlike nature of point defect cores and its importance on texture transformation mechanisms and stability. The paper shows that the ringlike cores can be oriented either along the cylinder axis or along the radial direction. Axial rings can partially expand but are constrained by the capillary sidewalls. Radial rings can deform into elliptical structures whose major axis is along the capillary axis. The transformation between several families of textures under capillary confinement as well as their stability is discussed in terms of defect ring distortions. A unified view of nematic textures found in the cylindrical cavities is provided. [ABSTRACT FROM AUTHOR]

In this paper we improve the product simultaneous diagonalization (SD) basis method we previously proposed [J. Chem. Phys. 122, 134101 (2005)] and applied to solve the Schrödinger equation for the motion of nuclei on a potential surface. The improved method is tested using coupled complicated Hamiltonians with as many as 16 coordinates for which we can easily find numerically exact solutions. In a basis of sorted products of one-dimensional (1D) SD functions the Hamiltonian matrix is nearly diagonal. The localization of the 1D SD functions for coordinate qc depends on a parameter we denote αc. In this paper we present a trace minimization scheme for choosing αc to nearly block diagonalize the Hamiltonian matrix. Near-block diagonality makes it possible to truncate the matrix without degrading the accuracy of the lowest energy levels. We show that in the sorted product SD basis perturbation theory works extremely well. The trace minimization scheme is general and easy to implement. [ABSTRACT FROM AUTHOR]

Accurate calculation of the energies and widths of the resonances of HOCl—an important intermediate in the O(1D)HCl reactive system—poses a challenging benchmark for computational methods. The need for very large direct product basis sets, combined with an extremely high density of states, results in difficult convergence for iterative methods. A recent calculation of the highly excited OH stretch mode resonances using the filter diagonalization method, for example, required 462 000 basis functions, and 180 000 iterations. In contrast, using a combination of new methods, we are able to compute the same resonance states to higher accuracy with a basis less than half the size, using only a few hundred iterations—although the CPU cost per iteration is substantially greater. Similar performance enhancements are observed for calculations of the high-lying bound states, as reported in a previous paper [J. Theo. Comput. Chem. 2, 583 (2003)]. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper, we extend and elaborate upon a wavelet method first presented in a previous publication [B. Poirier, J. Theo. Comput. Chem. 2, 65 (2003)]. In particular, we focus on construction and optimization of the wavelet functions, from theoretical and numerical viewpoints, and also examine their localization properties. The wavelets used are modified Wilson-Daubechies wavelets, which in conjunction with a simple phase space truncation scheme, enable one to solve the multidimensional Schrödinger equation. This approach is ideally suited to rovibrational spectroscopy applications, but can be used in any context where differential equations are involved. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper, we compare and contrast basis set sampling techniques recently developed for use in the ab initio multiple cloning method, a direct dynamics extension to the multiconfigurational Ehrenfest approach, used recently for the quantum simulation of ultrafast photochemistry. We demonstrate that simultaneous use of basis set cloning and basis function trains can produce results which are converged to the exact quantum result. To demonstrate this, we employ these sampling methods in simulations of quantum dynamics in the spin boson model with a broad range of parameters and compare the results to accurate benchmarks. [ABSTRACT FROM AUTHOR]

A recent quantum computing paper (G. S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007)) analytically derived optimal pulse spacings for a multiple spin echo sequence designed to remove decoherence in a two-level system coupled to a bath. The spacings in what has been called a “Uhrig dynamic decoupling (UDD) sequence” differ dramatically from the conventional, equal pulse spacing of a Carr–Purcell–Meiboom–Gill (CPMG) multiple spin echo sequence. The UDD sequence was derived for a model that is unrelated to magnetic resonance, but was recently shown theoretically to be more general. Here we show that the UDD sequence has theoretical advantages for magnetic resonance imaging of structured materials such as tissue, where diffusion in compartmentalized and microstructured environments leads to fluctuating fields on a range of different time scales. We also show experimentally, both in excised tissue and in a live mouse tumor model, that optimal UDD sequences produce different T2-weighted contrast than do CPMG sequences with the same number of pulses and total delay, with substantial enhancements in most regions. This permits improved characterization of low-frequency spectral density functions in a wide range of applications. [ABSTRACT FROM AUTHOR]

This paper presents dielectric relaxation data for organic glass-forming liquids compiled from different groups and supplemented by new measurements. The main quantity of interest is the “minimum slope” of the α dielectric loss plotted as a function of frequency in a log-log plot, i.e., the numerically largest slope above the loss peak frequency. The data consisting of 347 spectra for 53 liquids show prevalence of minimum slopes close to -1/2, corresponding to approximate t dependence of the dielectric relaxation function at short times. The paper studies possible correlations between minimum slopes and (1) temperature (quantified via the loss peak frequency); (2) how well an inverse power-law fits data above the loss peak; (3) degree of time-temperature superposition; (4) loss peak half width; (5) deviation from non-Arrhenius behavior; (6) loss strength. For the first three points we find correlations that show a special status of liquids with minimum slopes close to -1/2. For the last three points only fairly insignificant correlations are found, with the exception of large-loss liquids that have minimum slopes that are numerically significantly larger than 1/2. We conclude that—excluding large-loss liquids—approximate t relaxation appears to be a generic property of the α relaxation of organic glass formers. [ABSTRACT FROM AUTHOR]

We have constructed free-time and fixed end-point optimal control theory for quantum systems and applied it to entanglement generation between rotational modes of two polar molecules coupled by dipole-dipole interaction. The motivation of the present work is to solve optimal control problems more flexibly by extending the popular fixed time and fixed end-point optimal control theory for quantum systems to free-time and fixed end-point optimal control theory. As a demonstration, the theory that we have constructed in this paper will be applied to entanglement generation in rotational modes of NaCl–NaBr polar molecular systems that are sensitive to the strength of entangling interactions. Our method will significantly be useful for the quantum control of nonlocal interaction such as entangling interaction, which depends crucially on the strength of the interaction or the distance between the two molecules, and other general quantum dynamics, chemical reactions, and so on. [ABSTRACT FROM AUTHOR]

With the advent of nanotechnology, nano-electro-osmosis flow (nano-EOF) has shown great promise in the next generation of lab-on-a-chip systems such as attoliter or picoliter syringes or pipettors. In order to optimize the design of such systems and to precisely control analysis processes, it is important to obtain better fundamental understanding of EOF at nanoscales. Therefore, a more comprehensive electric double layer (EDL) theory is in need to improve upon the conventional EDL theory based on the Poisson–Boltzmann (PB) equation. In this paper, the modified PB theory is utilized to investigate the flow behavior of EOF at micro- and nanoscales. The effect of ion size and the image effect are particularly emphasized. Both effects remarkably influence nano-EOF. More importantly, this study predicts a new phenomenon: two “peaks” may appear in the velocity profile of nano-EOF for some specific solid-electrolyte systems. [ABSTRACT FROM AUTHOR]

Local coupled-cluster theory provides an algorithm for measuring electronic correlation quickly, using only the spatial locality of localized electronic orbitals. Previously, we showed [J. Subotnik et al., J. Chem. Phys. 125, 074116 (2006)] that one may construct a local coupled-cluster singles-doubles theory which (i) yields smooth potential energy surfaces and (ii) achieves near linear scaling. That theory selected which orbitals to correlate based only on the distances between the centers of different, localized orbitals, and the approximate potential energy surfaces were characterized as smooth using only visual identification. This paper now extends our previous algorithm in three important ways. First, locality is now based on both the distances between the centers of orbitals as well as the spatial extent of the orbitals. We find that, by accounting for the spatial extent of a delocalized orbital, one can account for electronic correlation in systems with some electronic delocalization using fast correlation methods designed around orbital locality. Second, we now enforce locality on not just the amplitudes (which measure the exact electron-electron correlation), but also on the two-electron integrals themselves (which measure the bare electron-electron interaction). Our conclusion is that we can bump integrals as well as amplitudes, thereby gaining a tremendous increase in speed and paradoxically increasing the accuracy of our LCCSD approach. Third and finally, we now make a rigorous definition of chemical smoothness as requiring that potential energy surfaces not support artificial maxima, minima, or inflection points. By looking at first and second derivatives from finite difference techniques, we demonstrate complete chemical smoothness of our potential energy surfaces (bumping both amplitudes and integrals). These results are significant both from a theoretical and from a computationally practical point of view. [ABSTRACT FROM AUTHOR]

After many unsuccessful efforts, the structure of solid CD4 III finally has been solved. In this paper, we examine if the known tunneling spectra are consistent with the orientational potentials at the two sites with different symmetries. To this end, we study the rotational kinetic energy of the molecules, construct appropriate pocket states for the tunneling problem, and set up a model potential. Approximate energy levels are obtained from the Ritz variational principle. The agreement between the experimentally determined tunneling frequencies and the calculations is rather good, corroborating the findings of the structural analysis. A continuation of this paper will deal with the partly deuterated methanes. [ABSTRACT FROM AUTHOR]

The authors have recorded the ν3 infrared spectrum of methane in helium nanodroplets using our cw infrared optical parametric oscillator. In a previous paper, Nauta and Miller [Chem. Phys. Lett. 350, 225 (2001)] reported the observation of the monomer rovibrational transitions of methane in helium nanodroplets. Here, they report the observation of additional absorption bands in the frequency range between 2990 and 3070 cm-1 blueshifted compared to the monomer transitions. They attribute these absorption features to phonon wings of individual rovibrational transitions, i.e., the simultaneous excitation of collective excitation modes of the quantum fluid and the rovibrational excitation of the methane monomer in the helium nanodroplet. [ABSTRACT FROM AUTHOR]

Protonation and subsequent intramolecular hydrogen bonding as methods to control chain structure and tune luminescence in heteroatomic conjugated polymers were reported experimentally [A. P. Monkman et al., J. Am. Chem. Soc. 124, 6049 (2002)]. In this paper, the structure and photophysical properties of the model teraryl compound of phenylene-pyridylene copolymer before and after protonation are theoretically studied with quantum chemistry methods. From the optimized ground states, intramolecular hydrogen bonding to the adjacent oxygen atom in the alkoxy substituent planarizes the backbone of the molecules, and the optimized detailed results of compound 9 before and after protonation, such as the dihedral angles between the central benzene and the two pyridyl rings, the bond lengths, and the bond angles, are consistent with the experimental results. From the results of the calculated excited states, the protonation and subsequent intramolecular hydrogen bonding result in the redshifts of the absorption, the increase of the ionization energy, the increase of the electron affinity, the decrease of the energy difference of the highest occupied molecular orbital and lowest unoccupied molecular orbital, the decrease of the binding gap, and the delocalization of the electron-hole coherence. The photophysical properties of compound 9 before and after protonation are further studied with a three-dimensional real-space analysis method of transition and charge difference densities (study transition dipole moment and charge transfer in the absorption and fluorescence processes) and two-dimensional real-space analysis method of transition density matrices (study the electron-hole coherence and the excitation delocalization). The calculated results show theoretically an insight understanding on the influence of the protonation and subsequent intramolecular hydrogen bonding to chain structure and photophysical properties. [ABSTRACT FROM AUTHOR]

In a recent paper [D. Babikov, J. Chem. Phys. 121, 7577 (2004)], quantum optimal control theory was applied to analyze the accuracy of quantum gates in a quantum computer based on molecular vibrational eigenstates. The effects of the anharmonicity parameter of the molecule, the target time of the pulse, and the penalty function on the accuracy of the qubit transformations were investigated. We demonstrate that the effects of all the molecular and laser-pulse parameters can be explained utilizing the analytical pulse area theorem, which originates from the standard two-level model. Moreover, by analyzing the difference between the optimal control theory results and those obtained using the pulse area theorem, it is shown that extremely high quantum gate fidelity can be achieved for a qubit system based on vibrational eigenstates. [ABSTRACT FROM AUTHOR]

For proton transfer reactions, the tunneling contributions to the rates are often much larger than thermally activated rates at temperatures of interest. A number of separable tunneling corrections have been proposed that capture the dependence of tunneling rates on barrier height and imaginary frequency size. However, the effects of reaction pathway curvature and barrier anharmonicity are more difficult to quantify. The nonseparable semiclassical transition state theory (TST) of Hernandez and Miller [Chem. Phys. Lett. 214, 129 (1993)] accounts for curvature and barrier anharmonicity, but it requires prohibitively expensive cubic and quartic derivatives of the potential energy surface at the transition state. This paper shows how the reaction path Hamiltonian can be used to approximate the cubic and quartic derivatives used in nonseparable semiclassical transition state theory. This enables tunneling corrections that include curvature and barrier anharmonicity effects with just three frequency calculations as required by a conventional harmonic transition state theory calculation. The tunneling correction developed here is nonseparable, but can be expressed as a thermal average to enable efficient Monte Carlo calculations. For the proton exchange reaction NH2+CH4⇔NH3+CH3, the nonseparable rates are very accurate at temperatures from 300 K up to about 1000 K where the TST rate itself begins to diverge from the experimental results. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

This paper presents a density functional study of mixed amphiphile solutions. The amphiphiles are modeled using four to seven fused hard spheres where the connectivity is such that the amphiphiles have a proper head group. We present calculations of lamellar and vesicle density distributions, paying particular attention to the morphology of the bilayer structures. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface–ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations. [ABSTRACT FROM AUTHOR]

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The above comment (in the following called “the Comment”) questions our recent paper on diffusion measurements with nonlinear spin echoes.1 These doubts are unfounded and may even be interpreted in a misleading way. There is no unjustified or “undocumented” approximation in our papers on diffusion measurements with nonlinear spin echoes. © 2002 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In their comment, Sanz and Miret-Artés (SMA) describe previous trajectory-based formalisms based on the quantum Hamilton-Jacobi (QHJ) formalism. In this reply, we highlight our unique contributions: the identification of the smallness of the quantum force in the complex QHJ and its solution using complex trajectories. SMA also raise the question of how the term locality should be used in quantum mechanics. We suggest that at least certain aspects of nonlocality can depend on the method used to solve the problem. [ABSTRACT FROM AUTHOR]

Comparisons are presented of experimental and theoretical studies of the rotationally inelastic scattering of CD3 radicals with H2 and D2 collision partners at respective collision energies of 680 ± 75 and 640 ± 60 cm-1. Close-coupling quantum-mechanical calculations performed using a newly constructed ab initio potential energy surface (PES) provide initial-to-final CD3 rotational level (n, k → n', k') integral and differential cross sections (ICSs and DCSs). The DCSs are compared with crossed molecular beam and velocity map imaging measurements of angular scattering distributions, which serve as a critical test of the accuracy of the new PES. In general, there is very good agreement between the experimental measurements and the calculations. The DCSs for CD3 scattering from both H2 and D2 peak in the forward hemisphere for n' = 2-4 and shift more to sideways and backward scattering for n' = 5. For n' = 6-8, the DCSs are dominated by backward scattering. DCSs for a particular CD3 n→n' transition have a similar angular dependence with either D2 or H2 as collision partner. Any differences between DCSs or ICSs can be attributed to mass effects because the PES is unchanged for CD3-H2 and CD3-D2 collisions. Further comparisons are drawn between the CD3-D2 scattering and results for CD3-He presented in our recent paper [O. Tkáč, A. G. Sage, S. J. Greaves, A. J. Orr-Ewing, P. J. Dagdigian, Q. Ma, and M. H. Alexander, Chem. Sci. 4, 4199 (2013)]. These systems have the same reduced mass, but are governed by different PESs. [ABSTRACT FROM AUTHOR]

In this paper we investigate the dynamic properties of the minimal Bell-Lavis (BL) water model and their relation to the thermodynamic anomalies. The BL model is defined on a triangular lattice in which water molecules are represented by particles with three symmetric bonding arms interacting through van der Waals and hydrogen bonds. We have studied the model diffusivity in different regions of the phase diagram through Monte Carlo simulations. Our results show that the model displays a region of anomalous diffusion which lies inside the region of anomalous density, englobed by the line of temperatures of maximum density. Further, we have found that the diffusivity undergoes a dynamic transition which may be classified as fragile-to-strong transition at the critical line only at low pressures. At higher densities, no dynamic transition is seen on crossing the critical line. Thus evidence from this study is that relation of dynamic transitions to criticality may be discarded. [ABSTRACT FROM AUTHOR]

A full-configuration-interaction nuclear orbital treatment has been recently developed as a benchmark quantum-chemistry-like method to study small doped 3He clusters [M. P. de Lara-Castells et al., J. Chem. Phys. 125, 221101 (2006)]. Our objective in this paper is to extend our previous study on (3He)N–Cl2(B) clusters, using an enhanced implementation that allows employing very large one-particle basis sets [M. P. de Lara-Castells et al., J. Chem. Phys. 131, 194101 (2009)], and apply the method to the (3He)N–Cl2(X) case, using both a semiempirical T-shaped and an ab initio He-dopant potential with minima at both T-shaped and linear conformations. Calculations of the ground and low-lying excited solvent states stress the key role played by the anisotropy of the He-dopant interaction in determining the global energies and the structuring of the 3He atoms around the dopant. Whereas 3He atoms are localized in a broad belt around the molecular axis in ground-state N-sized complexes with N=1–3, irrespective of using the T-shaped or the ab initio He-dopant potential function, the dopant species becomes fully coated by just four 3He atoms when the He-dopant potential also has a minimum at linear configurations. However, excited solvent states with a central ring-type clustering of the host molecule are found to be very close in energy with the ground state by using the ab initio potential function. A microscopic analysis of this behavior is provided. Additional simulations of the molecular rovibrational Raman spectra, also including excited solvent states, provide further insights into the importance of proper modeling the anisotropy of the He-dopant interaction in these weakly bound systems and of taking into account the low-lying excitations. [ABSTRACT FROM AUTHOR]

The nonlocal van der Waals density functional (vdW-DF) captures the essential physics of the dispersion interaction not only in the asymptotic regime but also for a general case of overlapping fragment densities. A balanced treatment of other energetic contributions, such as exchange, is crucial if we aim for accurate description of various properties of weakly bound systems. In this paper, the vdW-DF correlation functional is modified to make it better compatible with accurate exchange functionals. We suggest a slightly simplified construction of the nonlocal correlation, yielding more accurate asymptotic C6 coefficients. We also derive a gradient correction, containing a parameter that can be adjusted to suit an exchange functional of choice. We devise a particularly apt combination of exchange and correlation terms, which satisfies many important constraints and performs well for our benchmark tests. [ABSTRACT FROM AUTHOR]

In a molecular system of energy donors and acceptors, resonance energy transfer is the primary mechanism by means of which electronic energy is redistributed between molecules, following the excitation of a donor. Given a suitable geometric configuration it is possible to completely inhibit this energy transfer in such a way that it can only be activated by application of an off-resonant laser beam: this is the principle of optically controlled resonance energy transfer, the basis for an all-optical switch. This paper begins with an investigation of optically controlled energy transfer between a single donor and acceptor molecule, identifying the symmetry and structural constraints and analyzing in detail the dependence on molecular energy level positioning. Spatially correlated donor and acceptor arrays with linear, square, and hexagonally structured arrangements are then assessed as potential configurations for all-optical switching. Built on quantum electrodynamical principles the concept of transfer fidelity, a parameter quantifying the efficiency of energy transportation, is introduced and defined. Results are explored by employing numerical simulations and graphical analysis. Finally, a discussion focuses on the advantages of such energy transfer based processes over all-optical switching of other proposed forms. [ABSTRACT FROM AUTHOR]

The partially linearized (pl), fully size-extensive multireference (MR) coupled-cluster (CC) method, fully accounting for singles (S) and doubles (D) and approximately for a subset of primary higher than doubles, referred to as plMR CCSD, as well as its plMR CCSD(T) version corrected for secondary triples, as described in Part I of this paper [X. Li and J. Paldus, J. Chem. Phys. 128, 144118 (2008)], are applied to the problem of bond breaking in the HF, F2, H2O, and N2 molecules, as well as to the H4 model, using basis sets of a DZ or a cc-pVDZ quality that enable a comparison with the full configuration interaction (FCI) exact energies for a given ab initio model. A comparison of the performance of the plMR CCSD/CCSD(T) approaches with those of the reduced MR (RMR) CCSD/CCSD(T) methods, as well as with the standard single reference (SR) CCSD and CCSD(T) methods, is made in each case. For the H4 model and N2 we also compare our results with the completely renormalized (CR) CC(2,3) method [P. Piecuch and M. Włoch, J. Chem. Phys. 123, 224105 (2005)]. An important role of a proper choice of the model space for the MR-type methods is also addressed. The advantages and shortcomings of all these methods are pointed out and discussed, as well as their size-extensivity characteristics, in which case we distinguish supersystems involving noninteracting SR and MR subsystems from those involving only MR-type subsystems. Although the plMR-type approaches render fully size-extensive results, while the RMR CCSD may slightly violate this property, the latter method yields invariably superior results to the plMR CCSD ones and is more easy to apply in highly demanding cases, such as the triple-bond breaking in the nitrogen molecule. [ABSTRACT FROM AUTHOR]

We describe a fully size-extensive alternative of the reduced multireference (RMR) coupled-cluster (CC) method with singles (S) and doubles (D) that generates a subset of higher-than-pair cluster amplitudes, using linearized CC equations from the full CC chain, projected onto the corresponding higher-than-doubly excited configurations. This approach is referred to as partially linearized (pl) MR CCSD method and characterized by the acronym plMR CCSD. In contrast to a similar CCSDT-1 method [Y. S. Lee et al., J. Chem. Phys. 81, 5906 (1984)] this approach also considers higher than triples (currently up to hexuples), while focusing only on a small subset of such amplitudes, referred to as the primary ones. These amplitudes are selected using similar criteria as in RMR CCSD. An extension considering secondary triples via the standard (T)-type corrections, resulting in the plMR CCSD(T) method, is also considered. The relationship of RMR and plMR CCSD and CCSD(T) approaches is discussed, and their performance and characteristics are the subject of the subsequent Part II of this paper. [ABSTRACT FROM AUTHOR]

Nonequilibrium molecular dynamics simulations are performed on a dense simple dipolar fluid under a planar Couette shear flow. Shear generates heat, which is removed by thermostatting terms added to the equations of motion of the fluid particles. The spatial structure of simple fluids at high shear rates is known to depend strongly on the thermostatting mechanism chosen. Kinetic thermostats are either biased or unbiased: biased thermostats neglect the existence of secondary flows that appear at high shear rates superimposed upon the linear velocity profile of the fluid. Simulations that employ a biased thermostat produce a string phase where particles align in strings with hexagonal symmetry along the direction of the flow. This phase is known to be a simulation artifact of biased thermostatting, and has not been observed by experiments on colloidal suspensions under shear flow. In this paper, we investigate the possibility of using a suitably directed electric field, which is coupled to the dipole moments of the fluid particles, to stabilize the string phase. We explore several thermostatting mechanisms where either the kinetic or configurational fluid degrees of freedom are thermostated. Some of these mechanisms do not yield a string phase, but rather a shear-thickening phase; in this case, we find the influence of the dipolar interactions and external field on the packing structure, and in turn their influence on the shear viscosity at the onset of this shear-thickening regime. [ABSTRACT FROM AUTHOR]

In this paper, ice nanotubes confined in carbon nanotubes are investigated by molecular dynamics. The trigonal, square, pentagonal, and hexagonal water tubes are obtained, respectively. The current-voltage (I-V) curves of water nanotubes are found to be nonlinear, and fluctuations of conductance spectra of these ice nanotubes show that the transport properties of ice nanotubes are quite different from those of bulk materials. Our studies indicate that the conductance gap of ice nanotube is related to the difference value from the Fermi energy EF to the nearest molecular energy level E0. Increasing the diameter of a water molecular nanostructure results in the increase of the conductance. [ABSTRACT FROM AUTHOR]

The satellite transitions acquired in real time by magic angle spinning (STARTMAS) NMR experiment combines a train of pulses with sample rotation at the magic angle to refocus the first- and second-order quadrupolar broadening of spin I=3/2 nuclei in a series of echoes, while allowing the isotropic chemical and quadrupolar shifts to evolve. The result is real-time isotropic NMR spectra at high spinning rates using conventional MAS equipment. In this paper we describe in detail how STARTMAS data can be acquired and processed with ease on commercial equipment. We also discuss the advantages and limitations of the approach and illustrate the discussion with numerical simulations and experimental data from four different powdered solids. [ABSTRACT FROM AUTHOR]

In this paper we establish the temperature dependence of the kinetic coefficient associated with crystal growth into the supercooled liquid for a wide range of organic and inorganic materials. We show that the kinetic coefficient for crystal growth scales with the shear viscosity η as η-ξ and that the exponent depends systematically on the fragility of the liquid. The greater the fragility (i.e., deviation away from an Arrhenius temperature dependence for η), the larger the difference 1-ξ. We argue that this breakdown in scaling between the crystal growth kinetics and the viscosity is a manifestation of heterogeneous dynamics in supercooled liquids. In addition, we show that the absolute growth rate at intermediate viscosities is correlated with the entropy difference between the liquid and the crystal. [ABSTRACT FROM AUTHOR]

Presents corrections for errors in the article 'Simulation of the coexistence of a shearing liquid and a strained crystal,' published in the 'Journal of Chemical Physics,' 118, 4115, 2003 issue. Identification of typographical errors; Expressions used in the calculations presented in the paper.

Hybrid functionals are responsible for much of the utility of modern Kohn-Sham density functional theory. When rigorously applied to solid-state metallic and small band gap systems, however, the slow decay of their nonlocal Hartree-Fock-type exchange makes hybrids computationally challenging and introduces unphysical effects. This can be remedied by using a range-separated hybrid which only keeps short-range nonlocal exchange, as in the functional of Heyd et al. [J. Chem. Phys. 118, 8207 (2003)]. On the other hand, many molecular properties require full long-range nonlocal exchange, which can also be included by means of a range-separated hybrid such as the recently introduced LC-ωPBE functional [O. A. Vydrov and G. E. Scuseria, J. Chem. Phys. 125, 234109 (2006)]. In this paper, we show that a three-range hybrid which mainly includes middle-range Hartree-Fock-type exchange and neglects long- and short-range Hartree-Fock-type exchange yields excellent accuracy for thermochemistry, barrier heights, and band gaps, emphasizing that the middle-range part of the 1/r potential seems crucial to accurately model these properties. [ABSTRACT FROM AUTHOR]

According to linear response theory, all relaxation functions in the linear regime can be obtained using time correlation functions calculated under equilibrium. In this paper, we demonstrate that the cross correlations make a significant contribution to the partial stress relaxation functions in polymer melts. We present two illustrations in the context of polymer rheology using (1) Brownian dynamics simulations of a single chain model for entangled polymers, the slip-spring model, and (2) molecular dynamics simulations of a multichain model. Using the single chain model, we analyze the contribution of the confining potential to the stress relaxation and the plateau modulus. Although the idea is illustrated with a particular model, it applies to any single chain model that uses a potential to confine the motion of the chains. This leads us to question some of the assumptions behind the tube theory, especially the meaning of the entanglement molecular weight obtained from the plateau modulus. To shed some light on this issue, we study the contribution of the nonbonded excluded-volume interactions to the stress relaxation using the multichain model. The proportionality of the bonded/nonbonded contributions to the total stress relaxation (after a density dependent “colloidal” relaxation time) provides some insight into the success of the tube theory in spite of using questionable assumptions. The proportionality indicates that the shape of the relaxation spectrum can indeed be reproduced using the tube theory and the problem is reduced to that of finding the correct prefactor. [ABSTRACT FROM AUTHOR]

We report in this paper a quantum dynamics study for the reaction H+NH3→NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement. [ABSTRACT FROM AUTHOR]

This paper is concerned with the structural transition dynamics of the six-atom Morse cluster with zero total angular momentum, which serves as an illustrative example of the general reaction dynamics of isolated polyatomic molecules. It develops a methodology that highlights the interplay between the effects of the potential energy topography and those of the intrinsic geometry of the molecular internal space. The method focuses on the dynamics of three coarse variables, the molecular gyration radii. By using the framework of geometric mechanics and hyperspherical coordinates, the internal motions of a molecule are described in terms of these three gyration radii and hyperangular modes. The gyration radii serve as slow collective variables, while the remaining hyperangular modes serve as rapidly oscillating “bath” modes. Internal equations of motion reveal that the gyration radii are subject to two different kinds of forces: One is the ordinary force that originates from the potential energy function of the system, while the other is an internal centrifugal force. The latter originates from the dynamical coupling of the gyration radii with the hyperangular modes. The effects of these two forces often counteract each other: The potential force generally works to keep the internal mass distribution of the system compact and symmetric, while the internal centrifugal force works to inflate and elongate it. Averaged fields of these two forces are calculated numerically along a reaction path for the structural transition of the molecule in the three-dimensional space of gyration radii. By integrating the sum of these two force fields along the reaction path, an effective energy curve is deduced, which quantifies the gross work necessary for the system to change its mass distribution along the reaction path. This effective energy curve elucidates the energy-dependent switching of the structural preference between symmetric and asymmetric conformations. The present methodology should be of wide use for the systematic reduction of dimensionality as well as for the identification of kinematic barriers associated with the rearrangement of mass distribution in a variety of molecular reaction dynamics in vacuum. [ABSTRACT FROM AUTHOR]

In this paper, we further explore the applicability of the accelerated molecular dynamics simulation method using a bias potential. The method is applied to both simple model systems and real multidimensional systems. The method is also compared to replica exchange simulations in folding a small protein, Trp cage, using an all atom potential for the protein and an implicit model for the solvent. We show that the bias potential method allows quick searches of folding pathways. We also show that the choice of the bias potential has significant influence on the efficiency of the bias potential method. [ABSTRACT FROM AUTHOR]