Showing papers in "arXiv: Chemical Physics in 2010"
TL;DR: Results of similar quality than the best available MCTDH benchmark, which is based on a wavepacket with 4.6×10(7)time-dependent coefficients, are obtained with a much more compact wavefunction consisting of only 4.5× 10(5) coefficients and requiring a shorter computation time.
Abstract: The multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method is discussed and a fully general implementation for any number of layers based on the recursive ML-MCTDH algorithm given by Manthe [J. Chem. Phys. {\bf 128}, 164116 (2008)] is presented. The method is applied first to a generalized Henon-Heiles (HH) Hamiltonian. For 6D HH the overhead of ML-MCTDH makes the method slower than MCTDH, but for 18D HH ML-MCTDH starts to be competitive. We report as well 1458D simulations of the HH Hamiltonian using a seven layer scheme. The photoabsorption spectrum of pyrazine computed with the 24D Hamiltonian of Raab {\em et. al.} [J. Chem. Phys. {\bf 110}, 936 (1999)] provides a realistic molecular test case for the method. Quick and small ML-MCTDH calculations needing a fraction of the time and resources of reference MCTDH calculations provide already spectra with all the correct features. Accepting slightly larger deviations, the calculation can be accelerated to take only 7 minutes. When pushing the method towards convergence, results of similar quality than the best available MCTDH benchmark, which is based on a wavepacket with $4.6\times 10^7$ time-dependent coefficients, are obtained with a much more compact wavefunction consisting of only $4.5\times 10^5$ coefficients and requiring a shorter computation time.
295 citations
TL;DR: It is shown that conventional hybrid functionals over-delocalize excitons and underestimate quasiparticle energy gaps in the acene systems, and the importance of both a range-separated and asymptotically correct contribution of exchange in TDDFT for investigating optoelectronic and excitonic properties, even for these simple valence excitations.
Abstract: The optoelectronic and excitonic properties in a series of linear acenes (naphthalene up to heptacene) are investigated using range-separated methods within time-dependent density functional theory (TDDFT). In these rather simple systems, it is well-known that TDDFT methods using conventional hybrid functionals surprisingly fail in describing the low-lying La and Lb valence states, resulting in large, growing errors for the La state and an incorrect energetic ordering as a function of molecular size. In this work, we demonstrate that the range-separated formalism largely eliminates both of these errors and also provides a consistent description of excitonic properties in these systems. We further demonstrate that re-optimizing the percentage of Hartree-Fock exchange in conventional hybrids to match wavefunction-based benchmark calculations still yields serious errors, and a full 100% Hartree-Fock range separation is essential for simultaneously describing both of the La and Lb transitions. Based on an analysis of electron-hole transition density matrices, we finally show that conventional hybrid functionals overdelocalize excitons and underestimate quasiparticle energy gaps in the acene systems. The results of our present study emphasize the importance of both a range-separated and asymptotically-correct contribution of exchange in TDDFT for investigating optoelectronic and excitonic properties, even for these simple valence excitations.
226 citations
TL;DR: The results of the present study emphasize the importance of long-range exchange corrections in TDDFT for investigating the excited-state properties in solar cell dyes and find that the LC technique provides a consistent picture of charge-transfer excitations as a function of molecular size.
Abstract: The excited-state properties in a series of coumarin solar cell dyes are investigated with a long-range-corrected (LC) functional which asymptotically incorporates Hartree-Fock exchange. Using time-dependent density functional theory (TDDFT), we calculate excitation energies, oscillator strengths, and excited-state dipole moments in each of the dyes as a function of the range-separation paramenter, mu. To investigate the acceptable range of mu and assess the quality of the LC-TDDFT formalism, an extensive comparison is made between LC-BLYP excitation energies and approximate coupled cluster singles and doubles (CC2) calculations. When using a properly-optimized value of mu, we find that the LC technique provides a consistent picture of charge-transfer excitations as a function of molecular size. In contrast, we find that the widely-used B3LYP hybrid functional severely overestimates excited-state dipole moments and underestimates vertical excitations energies, especially for larger dye molecules. The results of the present study emphasize the importance of long-range exchange corrections in TDDFT for investigating the excited-state properties in solar cell dyes.
150 citations
TL;DR: First principles computational techniques suggest that the proposed reaction mechanism is viable if electron transfer to the catalyst is sufficiently fast and the variation in cobalt ion charge and spin states during bond breaking, DFT+U treatment of cobalt 3d orbitals, and the need for computing electrochemical potentials are emphasized.
Abstract: We apply first principles computational techniques to analyze the two-electron, multi-step, electrochemical reduction of CO2 to CO in water using cobalt porphyrin as a catalyst. Density Functional Theory calculations with hybrid functionals and dielectric continuum solvation are used to determine the steps at which electrons are added. This information is corroborated with ab initio molecular dynamics simulations in an explicit aqueous environment which reveal the critical role of water in stabilizing a key intermediate formed by CO2 bound to cobalt. Using potential of mean force calculations, the intermediate is found to spontaneously accept a proton to form a carboxylate acid group at pH<9.0, and the subsequent cleavage of a C-OH bond to form CO is exothermic and associated with a small free energy barrier. These predictions suggest that the proposed reaction mechanism is viable if electron transfer to the catalyst is sufficiently fast. The variation in cobalt ion charge and spin states during bond breaking, DFT+U treatment of cobalt 3d orbitals, and the need for computing electrochemical potentials are emphasized.
120 citations
TL;DR: In this article, double core hole electron spectroscopy for chemical analysis in terms of x-ray two-photon photoelectron spectroscope (XTPPS) was explored.
Abstract: We explore the potential of double core hole electron spectroscopy for chemical analysis in terms of x-ray two-photon photoelectron spectroscopy (XTPPS). The creation of deep single and double core vacancies induces significant reorganization of valence electrons. The corresponding relaxation energies and the interatomic relaxation energies are evaluated by CASSCF calculations. We propose a method how to experimentally extract these quantities by the measurement of single and double core-hole ionization potentials (IPs and DIPs). The influence of the chemical environment on these DIPs is also discussed for states with two holes at the same atomic site and states with two holes at two different atomic sites. Electron density difference between the ground and double core-hole states clearly shows the relaxations accompanying the double core-hole ionization. The effect is also compared with the sensitivity of single core hole ionization potentials (IPs) arising in single core hole electron spectroscopy. We have demonstrated the method for a representative set of small molecules LiF, BeO, BF, CO, N2, C2H2, C2H4, C2H6, CO2 and N2O. The scalar relativistic effect on IPs and on DIPs are briefly addressed.
89 citations
TL;DR: In this paper, an ensemble of molecules in the rotationally ground state trapped on an optical lattice exhibits collective spin excitations that can be controlled by applying superimposed electric and magnetic fields.
Abstract: We show that an ensemble of $^2\Sigma$ molecules in the rotationally ground state trapped on an optical lattice exhibits collective spin excitations that can be controlled by applying superimposed electric and magnetic fields. In particular, we show that the lowest energy excitation of the molecular ensemble at certain combinations of electric and magnetic fields leads to the formation of a magnetic Frenkel exciton. The exciton bandwidth can be tuned by varying the electric or magnetic fields. We show that the exciton states can be localized by creating vacancies in the optical lattice. The localization patterns of the magnetic exciton states are sensitive to the number and distribution of vacancies, which can be exploited for engineering many-body entangled spin states. We consider the dynamics of magnetic exciton wavepackets and show that the spin excitation transfer between molecules in an optical lattice can be accelerated or slowed down by tuning an external magnetic or electric field.
76 citations
TL;DR: In this paper, the ultraviolet photodissociation of amorphous water ice at different ice temperatures is investigated using molecular dynamics simulations and analytical potentials, showing that desorption dominates in the top layers, while trapping occurs deeper in the ice.
Abstract: The ultraviolet (UV) photodissociation of amorphous water ice at different ice temperatures is investigated using molecular dynamics (MD) simulations and analytical potentials. Previous MD calculations of UV photodissociation of amorphous and crystalline water ice at 10 K [S. Andersson et al., J. Chem. Phys. 124, 064715 (2006)] revealed -for both types of ice- that H atom, OH, and H2O desorption are the most important processes after photoexcitation in the uppermost layers of the ice. Water desorption takes place either by direct desorption of recombined water, or when, after dissociation, an H atom transfers part of its kinetic energy to one of the surrounding water molecules which is thereby kicked out from the ice. We present results of MD simulations of UV photodissociation of amorphous ice at 10, 20, 30, and 90 K in order to analyze the effect of ice temperature on UV photodissociation processes. Desorption and trapping probabilities are calculated for photoexcitation of H2O in the top four monolayers and the main conclusions are in agreement with the 10 K results: desorption dominates in the top layers, while trapping occurs deeper in the ice. The hydrogen atom photodesorption probability does not depend on ice temperature, but OH and H2O photodesorption probabilities tend to increase slightly (~30%) with ice temperature. We have compared the total photodesorption probability (OH+H2O) with the experimental total photodesorption yield, and in both cases the probabilities rise smoothly with ice temperature. The experimental yield is on average 3.8 times larger than our theoretical results, which can be explained by the different time scales studied and the approximations in our model.
76 citations
TL;DR: In this article, a tensor tensor network (CGTN) is proposed to capture the electron correlation of a molecule of arbitrary structure, which is applied to the energy splitting of states of different spin for methylene and the strongly correlated ozone molecule at a transition state structure.
Abstract: We present a new class of tensor network states that are specifically designed to capture the electron correlation of a molecule of arbitrary structure. In this ansatz, the electronic wave function is represented by a Complete-Graph Tensor Network (CGTN) ansatz which implements an efficient reduction of the number of variational parameters by breaking down the complexity of the high-dimensional coefficient tensor of a full-configuration-interaction (FCI) wave function. We demonstrate that CGTN states approximate ground states of molecules accurately by comparison of the CGTN and FCI expansion coefficients. The CGTN parametrization is not biased towards any reference configuration in contrast to many standard quantum chemical methods. This feature allows one to obtain accurate relative energies between CGTN states which is central to molecular physics and chemistry. We discuss the implications for quantum chemistry and focus on the spin-state problem. Our CGTN approach is applied to the energy splitting of states of different spin for methylene and the strongly correlated ozone molecule at a transition state structure. The parameters of the tensor network ansatz are variationally optimized by means of a parallel-tempering Monte Carlo algorithm.
67 citations
TL;DR: The performance of nine density functionals demonstrates that Zhao and Truhlar's MPWB1K and M05‐2X functionals provide energies similar to that obtained at the SCS‐MP2 level, and further support the experimental interpretation of pure π–π interactions between a convex fullerene and the concave surfaces of two corannulene subunits.
Abstract: Stimulated by the recent observation of {\pi}-{\pi} interactions between C60 and corannulene subunits in a molecular tweezer arrangement (J. Am. Chem. Soc. 2007, 129, 3842), a density functional theory study was performed to analyze the electronic structure and properties of various noncovalent corannulene complexes. The theoretical approach is first applied to corannulene complexes with a series of benchmark molecules (CH4, NH3, and H2O) using several new-generation density functionals. The performance of nine density functionals, illustrated by computing binding energies of the corannulene complexes, demonstrates that Zhao and Truhlar's MPWB1K and M05-2X functionals provide energies similar to that obtained at the SCS-MP2 level. In contrast, most of the other popular density functionals fail to describe this noncovalent interaction or yield purely repulsive interactions. Further investigations with the M05-2X functional show that the binding energy of C60 with corannulene subunits in the relaxed molecular receptor clip geometry is -20.67 kcal/mol. The results of this calculation further support the experimental interpretation of pure {\pi}-{\pi} interactions between a convex fullerene and the concave surfaces of two corannulene subunits.
60 citations
TL;DR: In this paper, the properties of finite gold nanocones as optical antennas for enhancing molecular fluorescence were investigated, and it was shown that the maximum field and fluorescence enhancements do not occur for the same nanocone parameters.
Abstract: We investigate the properties of finite gold nanocones as optical antennas for enhancing molecular fluorescence. We compute the modification of the excitation rate, spontaneous emission rate, and quantum efficiency as a function of the nanocone base and length, showing that the maximum field and fluorescence enhancements do not occur for the same nanocone parameters. We compare the results with those for nanorods and nanospheroids and find that nanocones perform better.
54 citations
TL;DR: In this article, the static hyperpolarizability of CHCl_3 was calculated using three different kinds of basis sets: Gaussian-type orbitals, numerical basis sets, and real-space grids.
Abstract: Calculations of the hyperpolarizability are typically much more difficult to converge with basis set size than the linear polarizability. In order to understand these convergence issues and hence obtain accurate ab initio values, we compare calculations of the static hyperpolarizability of the gas-phase chloroform molecule (CHCl_3) using three different kinds of basis sets: Gaussian-type orbitals, numerical basis sets, and real-space grids. Although all of these methods can yield similar results, surprisingly large, diffuse basis sets are needed to achieve convergence to comparable values. These results are interpreted in terms of local polarizability and hyperpolarizability densities. We find that the hyperpolarizability is very sensitive to the molecular structure, and we also assess the significance of vibrational contributions and frequency dispersion.
TL;DR: In this article, the results of a combined metadynamics-umbrella sampling investigation of the puckered conformers of pyranoses described using the gromos 45a4 force field are presented.
Abstract: We present the results of a combined metadynamics--umbrella sampling investigation of the puckered conformers of pyranoses described using the gromos 45a4 force field. The free energy landscape of Cremer--Pople puckering coordinates has been calculated for the whole series of alpha and beta aldohexoses, showing that the current force field parameters fail in reproducing proper puckering free energy differences between chair conformers. We suggest a modification to the gromos 45a4 parameter set which improves considerably the agreement of simulation results with theoretical and experimental estimates of puckering free energies. We also report on the experimental measurement of altrose conformers populations by means of NMR spectroscopy, which show good agreement with the predictions of current theoretical models.
TL;DR: In this article, a novel, highly flexible, conductive polymer-based fiber with high electric capacitance is reported, which is fabricated using fiber drawing method, where a multi-material macroscopic preform is drawn into a sub-millimeter capacitor fiber in a single fabrication step.
Abstract: A novel, highly flexible, conductive polymer-based fiber with high electric capacitance is reported. In its crossection the fiber features a periodic sequence of hundreds of conductive and isolating plastic layers positioned around metallic electrodes. The fiber is fabricated using fiber drawing method, where a multi-material macroscopic preform is drawn into a sub-millimeter capacitor fiber in a single fabrication step. Several kilometres of fibers can be obtained from a single preform with fiber diameters ranging between 500um -1000um. A typical measured capacitance of our fibers is 60-100 nF/m and it is independent of the fiber diameter. For comparison, a coaxial cable of the comparable dimensions would have only ~0.06nF/m capacitance. Analysis of the fiber frequency response shows that in its simplest interrogation mode the capacitor fiber has a transverse resistance of 5 kOhm/L, which is inversely proportional to the fiber length L and is independent of the fiber diameter. Softness of the fiber materials, absence of liquid electrolyte in the fiber structure, ease of scalability to large production volumes, and high capacitance of our fibers make them interesting for various smart textile applications ranging from distributed sensing to energy storage.
Posted Content•
TL;DR: In this article, a new electronic structure principle, viz., the principle of electrophilicity equalization, is proposed and an analytical justification as well as a numerical support for the same is provided.
Abstract: A new electronic structure principle, viz., the principle of electrophilicity equalization is proposed. An analytical justification as well as a numerical support for the same is provided.
TL;DR: In this article, an approach combining Becke and Johnson's [J. Phys. 127, 154108 (2007)] method for the evaluation of dispersion energy corrections and a Hirshfeld method for partitioning of molecular polarizability tensors into atomic contributions is presented.
Abstract: This paper presents an approach for obtaining accurate interaction energies at the DFT level for systems where dispersion interactions are important. This approach combines Becke and Johnson's [J. Chem. Phys. 127, 154108 (2007)] method for the evaluation of dispersion energy corrections and a Hirshfeld method for partitioning of molecular polarizability tensors into atomic contributions. Due to the availability of atomic polarizability tensors, the method is extended to incorporate anisotropic contributions, which prove to be important for complexes of lower symmetry. The method is validated for a set of eighteen complexes, for which interaction energies were obtained with the B3LYP, PBE and TPSS functionals combined with the aug-cc-pVTZ basis set and compared with the values obtained at CCSD(T) level extrapolated to a complete basis set limit. It is shown that very good quality interaction energies can be obtained by the proposed method for each of the examined functionals, the overall performance of the TPSS functional being the best, which with a slope of 1.00 in the linear regression equation and a constant term of only 0.1 kcal/mol allows to obtain accurate interaction energies without any need of a damping function for complexes close to their exact equilibrium geometry.
TL;DR: In this article, a high resolution spectroscopy of asymptotic vibrational levels and the use of available photo-association data were performed by high-resolution spectrograms of the ground state X^1\Sigma^+_g of Sr_2.
Abstract: An experimental study of the long range behavior of the ground state X^1\Sigma^+_g of Sr_2 is performed by high resolution spectroscopy of asymptotic vibrational levels and the use of available photoassociation data Ground state levels as high as v"=60 (outer turning point at 23 \AA and 01 cm^-1 below the asymptote) could be observed by Fourier-transform spectroscopy of fluorescence progressions induced by single frequency laser excitation of the v'=4, J'=9 rovibrational level of the state 2^1\Sigma^+_u A precise value of the scattering length for the isotopologue ^88Sr_2 is derived and transferred to all other isotopic combinations by mass scaling with the given potential The derived potential together with already published information about the state 2^1\Sigma^+_u directs to promising optical paths for producing cold molecules in the electronic ground state from an ultracold ensemble of Sr atoms
TL;DR: In this article, the generalized mass action law together with the basic relations between kinetic factors are proved for the positivity of the entropy production but hold even without microreversibility, when the detailed balance is not applicable.
Abstract: We study chemical reactions with complex mechanisms under two assumptions: (i) intermediates are present in small amounts (this is the quasi-steady-state hypothesis or QSS) and (ii) they are in equilibrium relations with substrates (this is the quasiequilibrium hypothesis or QE). Under these assumptions, we prove the generalized mass action law together with the basic relations between kinetic factors, which are sufficient for the positivity of the entropy production but hold even without microreversibility, when the detailed balance is not applicable. Even though QE and QSS produce useful approximations by themselves, only the combination of these assumptions can render the possibility beyond the "rarefied gas" limit or the "molecular chaos" hypotheses. We do not use any a priori form of the kinetic law for the chemical reactions and describe their equilibria by thermodynamic relations. The transformations of the intermediate compounds can be described by the Markov kinetics because of their low density ({\em low density of elementary events}). This combination of assumptions was introduced by Michaelis and Menten in 1913. In 1952, Stueckelberg used the same assumptions for the gas kinetics and produced the remarkable semi-detailed balance relations between collision rates in the Boltzmann equation that are weaker than the detailed balance conditions but are still sufficient for the Boltzmann $H$-theorem to be valid. Our results are obtained within the Michaelis-Menten-Stueckelbeg conceptual framework.
TL;DR: In this article, the potential distribution framework was used to calculate the excess free energy of liquid water simulated with the BLYP-D density functional, and the calculated free energy is in fair agreement with experiments but the excess internal energy and hence also the excess entropy are not.
Abstract: We regularize the potential distribution framework to calculate the excess free energy of liquid water simulated with the BLYP-D density functional. The calculated free energy is in fair agreement with experiments but the excess internal energy and hence also the excess entropy are not. Our work emphasizes the importance of thermodynamic characterization in assessing the quality of electron density functionals in describing liquid water and hydration phenomena.
TL;DR: Three isostructural cyanido-bridged heptanuclear complexes, with unique octameric water cluster, are obtained by reacting the dinuclear copper(II) complex.
Abstract: Three isostructural cyano-bridged heptanuclear complexes, [{CuII(saldmen)(H2O)}6{MIII(CN)6}](ClO4)3$\cdotp$8H2O (M = FeIII 2; CoIII, 3; CrIII 4), have been obtained by reacting the binuclear copper(II) complex, [Cu2(saldmen)2(mu-H2O)(H2O)2](ClO4)2$\cdotp$2H2O 1, with K3[Co(CN)6], K4[Fe(CN)6], and, respectively, K3[Cr(CN)6] (Hsaldmen is the Schiff base resulted from the condensation of salicylaldehyde with N,N-dimethylethylenediamine). A unique octameric water cluster, with bicyclo[2,2,2]octane-like structure, is sandwiched between the heptanuclear cations in 2, 3 and 4. The cryomagnetic investigations of compounds 2 and 4 reveal ferromagnetic couplings of the central FeIII or CrIII ions with the CuII ions (JCuFe = +0.87 cm-1, JCuCr = +30.4 cm-1). The intramolecular Cu-Cu exchange interaction in 3, across the diamagnetic cobalt(III) ion, is -0.3 cm-1. The solid-state1H-NMR spectra of compounds 2 and 3 have been investigated.
TL;DR: In this article, the authors apply modern surface and bulk techniques to the study of paraffin coatings, in order to characterize the properties that enable the effective preservation of alkali spin polarization.
Abstract: Many technologies based on cells containing alkali-metal atomic vapor benefit from the use of anti-relaxation surface coatings in order to preserve atomic spin polarization. In particular, paraffin has been used for this purpose for several decades and has been demonstrated to allow an atom to experience up to 10,000 collisions with the walls of its container without depolarizing, but the details of its operation remain poorly understood. We apply modern surface and bulk techniques to the study of paraffin coatings, in order to characterize the properties that enable the effective preservation of alkali spin polarization. These methods include Fourier transform infrared spectroscopy, differential scanning calorimetry, atomic force microscopy, near-edge X-ray absorption fine structure spectroscopy, and X-ray photoelectron spectroscopy. We also compare the light-induced atomic desorption yields of several different paraffin materials. Experimental results include the determination that crystallinity of the coating material is unnecessary, and the detection of C=C double bonds present within a particular class of effective paraffin coatings. Further study should lead to the development of more robust paraffin anti-relaxation coatings, as well as the design and synthesis of new classes of coating materials.
TL;DR: In this article, the authors present a spatial version of the partitioned-leaping algorithm (PLA), a multiscale accelerated-stochastic simulation approach built upon the tau-leaving framework of Gillespie.
Abstract: There is a great need for accurate and efficient computational approaches that can account for both the discrete and stochastic nature of chemical interactions as well as spatial inhomogeneities and diffusion. This is particularly true in biology and nanoscale materials science, where the common assumptions of deterministic dynamics and well-mixed reaction volumes often break down. In this article, we present a spatial version of the partitioned-leaping algorithm (PLA), a multiscale accelerated-stochastic simulation approach built upon the tau-leaping framework of Gillespie. We pay special attention to the details of the implementation, particularly as it pertains to the time step calculation procedure. We point out conceptual errors that have been made in this regard in prior implementations of spatial tau-leaping and illustrate the manifestation of these errors through practical examples. Finally, we discuss the fundamental difficulties associated with incorporating efficient exact-stochastic techniques, such as the next-subvolume method, into a spatial-leaping framework and suggest possible solutions.
TL;DR: In this article, two key theoretical approaches to investigate the excited state properties of organic light emitting diodes (OLEDs) and organic photovoltaics (OPV) are discussed: first principles atomistic models and effective Hamiltonian models.
Abstract: Organometallic complexes have potential applications as the optically active components of organic light emitting diodes (OLEDs) and organic photovoltaics (OPV). Development of more effective complexes may be aided by understanding their excited state properties. Here we discuss two key theoretical approaches to investigate these complexes: first principles atomistic models and effective Hamiltonian models. We review applications of these methods, such as, determining the nature of the emitting state, predicting the fraction of injected charges that form triplet excitations, and explaining the sensitivity of device performance to small changes in the molecular structure of the organometallic complexes.
TL;DR: In this paper, the interaction potential between N atoms and NH molecules is calculated and the dependence of the scattering properties on the reduced mass and explain the resonant effects observed using angular-momentum-insensitive quantum defect theory.
Abstract: We calculate the interaction potential between N atoms and NH molecules and use it to investigate cold and ultracold collisions important for sympathetic cooling. The ratio of elastic to inelastic cross sections is large over a wide range of collision energy and magnetic field for most isotopic combinations, so that sympathetic cooling of NH molecules by N atoms is a good prospect. However, there are important effects due to a p-wave shape resonance that may inhibit cooling in some cases. We show that scaling the reduced mass used in the collision is approximately equivalent to scaling the interaction potential. We then explore the dependence of the scattering properties on the reduced mass and explain the resonant effects observed using angular-momentum-insensitive quantum defect theory
TL;DR: In this paper, an approximate method for evaluating the importance of non-Born-Oppenheimer effects on the quantum dynamics of nuclei is proposed, which can be implemented easily into any molecular dynamics program and also can utilize on-the-fly ab initio electronic structure information.
Abstract: We propose an approximate method for evaluating the importance of non-Born-Oppenheimer effects on the quantum dynamics of nuclei. The method uses a generalization of the dephasing representation (DR) of quantum fidelity to several diabatic potential energy surfaces and its computational cost is the cost of dynamics of a classical phase space distribution. It can be implemented easily into any molecular dynamics program and also can utilize on-the-fly ab initio electronic structure information. We test the methodology on three model problems introduced by Tully and on the photodissociation of NaI. The results show that for dynamics close to the diabatic limit the decay of fidelity due to nondiabatic effects is described accurately by the DR. In this regime, unlike the mixed quantum-classical methods such as surface hopping or Ehrenfest dynamics, the DR can capture more subtle quantum effects than the population transfer between potential energy surfaces. Hence we propose using the DR to estimate the dynamical importance of diabatic, spin-orbit, or other couplings between potential energy surfaces. The acquired information can help reduce the complexity of a studied system without affecting the accuracy of the quantum simulation.
TL;DR: In this article, the authors derive structure-property relationships for methine ("Brooker") dyes relating the color of the dye and its symmetric parents to its bond alternation in the ground state and also to the dipole properties associated with its low-lying charge-resonance (or charge-transfer) transition.
Abstract: We derive structure-property relationships for methine ("Brooker") dyes relating the color of the dye and its symmetric parents to its bond alternation in the ground state and also to the dipole properties associated with its low-lying charge-resonance (or charge-transfer) transition. We calibrate and test these relationships on an array of different protonation states of the green fluorescent protein chromophore motif (an asymmetric halochromic methine dye) and its symmetric parent dyes. The relationships rely on the assumption that the diabatic states that define the Platt model for methine dye color [J.R. Platt, J. Chem. Phys. 25 80 (1956)] can also be distinguished by their single-double bond alternation and by their charge localization character. These assumptions are independent of the primary constraint that defines the diabatic states in the Platt model - specifically, the Brooker deviation rule for methine dyes [L.G.S. Brooker, Rev. Mod. Phys. 14 275 (1942)]. Taking these assumptions, we show that the Platt model offers an alternate route to known structure-property relationships between the bond length alternation and the quadratic nonlinear polarizability {\beta}. We show also that the Platt model can be parameterized without the need for synthesis of the symmetric parents of a given dye, using dipole data obtained through spectroscopic measurements. This suggests that the Platt model parameters may be used as independent variables in free-energy relationships for chromophores whose symmetric parents cannot be synthesized or chromophores strongly bound to biomolecular environments. The latter category includes several recently characterized biomolecular probe constructs. We illustrate these concepts by an analysis of previously reported electroabsorption and second-harmonic generation experiments on green fluorescent proteins.
TL;DR: The kinetics of the self-assembly of nanocomponents into a virus, nanocapsule, or other composite structure is analyzed via a multiscale approach to achieve predictability and to preserve key atomic-scale features that underlie the formation and stability of the composite structures.
Abstract: The kinetics of the self-assembly of nanocomponents into a virus, nanocapsule, or other composite structure is analyzed via a multiscale approach. The objective is to achieve predictability and to preserve key atomic-scale features that underlie the formation and stability of the composite structures. We start with an all-atom description, the Liouville equation, and the order parameters characterizing nanoscale features of the system. An equation of Smoluchowski type for the stochastic dynamics of the order parameters is derived from the Liouville equation via a multiscale perturbation technique. The self-assembly of composite structures from nanocomponents with internal atomic structure is analyzed and growth rates are derived. Applications include the assembly of a viral capsid from capsomers, a ribosome from its major subunits, and composite materials from fibers and nanoparticles. Our approach overcomes errors in other coarse-graining methods which neglect the influence of the nanoscale configuration on the atomistic fluctuations. We account for the effect of order parameters on the statistics of the atomistic fluctuations which contribute to the entropic and average forces driving order parameter evolution. This approach enables an efficient algorithm for computer simulation of self-assembly, whereas other methods severely limit the timestep due to the separation of diffusional and complexing characteristic times. Given that our approach does not require recalibration with each new application, it provides a way to estimate assembly rates and thereby facilitate the discovery of self-assembly pathways and kinetic dead-end structures.
TL;DR: The values of molecular polarizabilities and softnesses of the 20 amino acids were computed ab initio (MP2) and it is demonstrated that the global softness is not simply proportional to the global polarizability in disagreement with the intuition that "a softer moiety is also more polarizable".
Abstract: The values of molecular polarizabilities and softnesses of the twenty amino acids were computed ab initio (MP2). By using the iterative Hirshfeld scheme to partition the molecular electronic properties, we demonstrate that the values of the softness of the side chain of the twenty amino acid are clustered in groups reflecting their biochemical classification, namely: aliphatic, basic, acidic, sulfur containing, and aromatic amino acids . The present findings are in agreement with previous results using different approximations and partitioning schemes [P. Senet and F. Aparicio, J. Chem. Phys. 126,145105 (2007)]. In addition, we show that the polarizability of the side chain of an amino acid depends mainly on its number of electrons (reflecting its size) and consequently cannot be used to cluster the amino acids in different biochemical groups, in contrast to the local softness. Our results also demonstrate that the global softness is not simply proportional to the global polarizability in disagreement with the intuition that "a softer moiety is also more polarizable". Amino acids with the same softness may have a polarizability differing by a factor as large as 1.7. This discrepancy can be understood from first principles as we show that the molecular polarizability depends on a "softness dipole vector" and not simply on the global softness.
TL;DR: In this paper, first-principle calculations using a van der Waals density functional (vdW-DF), [Phys. Rev. Lett. $\mathbf{92}$, 246401 (2004)] permits determination of molecular crystal structure.
Abstract: In this paper, we show that first-principle calculations using a van der Waals density functional (vdW-DF), [Phys. Rev. Lett. $\mathbf{92}$, 246401 (2004)] permits determination of molecular crystal structure. We study the crystal structures of hexamine and the platonic hydrocarbons (cubane and dodecahedrane). The calculated lattice parameters and cohesion energy agree well with experiments. Further, we examine the asymptotic accounts of the van der Waals forces by comparing full vdW-DF with asymptotic atom-based pair potentials extracted from vdW-DF. The character of the binding differ in the two cases, with vdW-DF giving a significant enhancement at intermediate and relevant binding separations. We analyze consequences of this result for methods such as DFT-D, and question DFT-D's transferability over the full range of separations.
TL;DR: In this paper, photostop has been used to produce high densities of trappable molecules with a narrow velocity distribution centered at zero in the laboratory frame, where NO2, initially cooled in a pulsed molecular beam, is photodissociated such that the recoil velocity of the NO photofragments cancels out the velocity in the beam.
Abstract: We have demonstrated a new, accessible and economical technique, dubbed photostop, for producing high densities of trappable molecules. Direct measurements are presented of NO molecules produced with a narrow velocity distribution centered at zero in the laboratory frame. NO2, initially cooled in a pulsed molecular beam, is photodissociated such that the recoil velocity of the NO photofragments cancels out the velocity of the beam. NO(X^2Pi_3/2, v=0, J=1.5) molecules are observed up to 10 mircoseconds after the dissociation event in the probe volume at an estimated density of 1E7 cm-3 per quantum state and at a translational temperature of 1.6 K. Through the choice of suitable precursors, photostop has the potential to extend the list atoms and molecules that can be slowed or trapped. It should be possible to accumulate density in a trap through consecutive loading of multiple pulses.
TL;DR: In this paper, the authors used the hyperspherical quantum reactive scattering method developed by Launay et al. for the collision of open shell S(1D) atoms with ortho- and para-hydrogen, in the kinetic energy range 1-120 K.
Abstract: Reactive cross-sections for the collision of open shell S(^1D) atoms with ortho- and para-hydrogen, in the kinetic energy range 1-120 K, have been calculated using the hyperspherical quantum reactive scattering method developed by Launay et al. [Chem. Phys. Lett. 169, 473 (1990)]. Short-range interactions, described using the {\em ab initio} potential energy surface by Ho et al. were complemented with an accurate description of the long-range interactions, where the main electrostatic (~R^{-5}) and dispersion (~ R^{-6}) contributions were considered. This allows the comparison with recent experimental measurements of rate constants and excitation functions for the title reaction at low temperatures [Berteloite et al., Phys. Rev. Lett., 2010]. The agreement is fairly good. The behavior in the considered energy range can be understood on the average in terms of a classical Langevin (capture) model, where the centrifugal barriers determine the amount of reactive flux which reaches the barrierless transition state. Additionally, the structure of the van der Waals well provides temporal trapping at short distances thus allowing the system to find its way to the reaction at some classically-forbidden energies. Interestingly, the cross-section for para-hydrogen shows clearly oscillating features associated to the opening of new partial waves and to shape resonances which may be amenable to experimental detection.