ALBERT

Description: Using velocity map ion imaging, the photodissociation of chlorobromomethane (CH 2 BrCl) at 233–234 nm has been studied. The total translational energy distributions and the anisotropy parameters have been determined from the ion images of the photofragments Br ( 2 P 1/2 ) (denoted as Br * ) and Br ( 2 P 3/2 ) (denoted as Br) for the dominant CH 2 BrCl + hv → CH 2 Cl + Br * and CH 2 BrCl + hv → CH 2 Cl + Br channels. Using an impulsive model invoking angular momentum conservation, the vibrational energy distributions of the chloromethyl radicals have been derived from the total translational energy distributions for the two channels. The study suggests that there are a number of vibrational modes of the chloromethyl radical to be excited in both of the two photodissociation channels. In the Br* channel, the CH 2 s-stretch mode v 1 has the most probability of excitation. While in the Br channel, the CH 2 scissors mode ν 2 is attributed to the highest peak of the vibrational energy curve of the chloromethyl radical. The results further imply that, following absorption of one UV photon of 234 nm, other vibrational modes besides v 5 (C–Br stretch mode) are also excited in the parent molecule.

Description: A definition of the triplet-triplet energy transfer reaction coordinate within the very weak electronic coupling limit is proposed, and a novel theoretical formalism is developed for its quantitative determination in terms of internal coordinates The present formalism permits (i) the separation of donor and acceptor contributions to the reaction coordinate, (ii) the identification of the intrinsic role of donor and acceptor in the triplet energy transfer process, and (iii) the quantification of the effect of every internal coordinate on the transfer process. This formalism is general and can be applied to classical as well as to nonvertical triplet energy transfer processes. The utility of the novel formalism is demonstrated here by its application to the paradigm of nonvertical triplet-triplet energy transfer involving cis -stilbene as acceptor molecule. In this way the effect of each internal molecular coordinate in promoting the transfer rate, from triplet donors in the low and high-energy limit, could be analyzed in detail.

Description: The fluorescence excitation spectra of jet-cooled benzocyclobutane have been recorded and together with its ultraviolet absorption spectra have been used to assign the vibrational frequencies for this molecule in its S 1 (π,π * ) electronic excited state. Theoretical calculations at the CASSCF(6,6)/aug-cc-pVTZ level of theory were carried out to compute the structure of the molecule in its excited state. The calculated structure was compared to that of the molecule in its electronic ground state as well as to the structures of related molecules in their S 0 and S 1 (π,π * ) electronic states. In each case the decreased π bonding in the electronic excited states results in longer carbon-carbon bonds in the benzene ring. The skeletal vibrational frequencies in the electronic excited state were readily assigned and these were compared to the ground state and to the frequencies of five similar molecules. The vibrational levels in both S 0 and S 1 (π,π * ) states were remarkably harmonic in contrast to the other bicyclic molecules. The decreases in the frequencies of the out-of-plane skeletal modes reflect the increased floppiness of these bicyclic molecules in their S 1 (π,π * ) excited state.

Description: Unimolecular decay of sample ions imposes a limit on the usable laser fluence in matrix-assisted laser desorption/ionization (MALDI) ion sources. Traditionally, some modest degree of collisional sample ion cooling has been achieved by connecting MALDI ion sources directly to gas-filled radio frequency (RF) multipoles. It was also discovered in the early 1990s that gas-filled RF multipoles exhibit increased ion transmission efficiency due to collisional ion focusing effects. This unexpected experimental finding was later supported by elementary Monte Carlo simulations. Both experiments and simulations assumed a resting background gas with typical pressures of the order of 1 Pa. However, considerable additional improvements can be achieved if laser desorbed sample ions are introduced immediately after desorption, still within the ion source, in an axisymmetric rarefied supersonic gas jet with peak pressure of the order of 100 Pa and flow velocities 〉300 m/s, and under weak electric fields. We describe here the design principle and report performance data of an ion source coined “MALDI-2,” which incorporates elements of both rarefied aerodynamics and particle optics. Such a design allows superb suppression of metastable fragmentation due to rapid collisional cooling in

Description: The influence of the soft interfacial volume fraction on physical properties of composite materials has been found to be significant. However, the soft interfacial volume fraction is difficultly determined by traditional experimental methods and simple models proposed so far. This article addresses the problem by means of theoretical and numerical approaches that start at a microscopic scale of composite materials, which are regarded as a three-phase composite structure with polydisperse convex particles, soft interfaces, and a matrix. A theoretical scheme for the soft interfacial volume fraction is proposed by a theory of the nearest-surface distribution functions and geometrical configurations of polydisperse convex particles. The theoretical scheme represents a generalized model for the soft interfacial volume fraction in that it cannot only determine the interfacial volume fraction around convex polyhedral particles but also to derive that around ellipsoidal and spherical particles. In order to test the theoretical scheme, a numerical model that adopts the three-phase composite structure and a numerical Monte Carlo integration scheme is presented. Also, theoretical and numerical results of the soft interfacial volume fraction around ellipsoidal and spherical particles in the literature are further compared. By way of application, it is shown that the developed model provides a quantitative means to evaluate the dependence of the soft interfacial volume fraction on various factors, such as geometrical configurations of particles and the interfacial thickness.

Description: We present resonance Raman measurements of crystalline trishomocubane and diamantane dimers containing a C=C double bond. Raman spectra were recorded with excitation energies between 2.33 eV and 5.42 eV. The strongest enhancement is observed for the C=C stretch vibration and a bending mode involving the two carbon atoms of the C=C bond, corresponding to the B 2 g wagging mode of ethylene. This is associated with the localization of the π-HOMO and LUMO and the elongation of the C=C bond length and a pyramidalization of the two sp 2 -hybridized carbon atoms at the optical excitation. The observed Raman resonance energies of the trishomocubane and diamantane dimers are significantly lower than the HOMO-LUMO gaps of the corresponding unmodified diamondoids.

Description: The state-to-state dynamics of the H * ( n ) + HD → D * ( n ′ ) + H 2 reactive scattering at the collision energy of 0.5 eV have been carried out for the first time by using H-atom Rydberg tagging time-of-flight technique. Experimental results show that the angular distribution of the total H 2 products presents clearly forward-backward asymmetric, which considerably differs from that of the corresponding H + + HD → D + + H 2 reaction predicted by previously theoretical calculations. Such disagreement between these two processes suggests that the Fermi independent-collider model is also not valid in describing the dynamics of isotopic variants of the H * + H 2 reaction. The rotational state distribution of the H 2 products demonstrates a saw-toothed distribution with odd- j ′ 〉 even- j ′ . This interesting observation is strongly influenced by nuclear spin statistics.

Description: We report high-resolution anion photoelectron spectra of vanadium dioxide (VO 2 − ) obtained by slow electron velocity-map imaging of trapped and cryogenically cooled anions. Vibrationally resolved spectra are obtained for photodetachment to the first three neutral electronic states, giving an electron affinity of 1.8357(5) eV for the X ̃ A 1 2 ground state and term energies of 0.1845(8) eV and 0.8130(5) eV for the A ̃ B 1 2 and B ̃ A 1 2 excited states, respectively. The vibrational fundamentals ν 1 and ν 2 are obtained for all three states. Experimental assignments are confirmed by energies from electronic structure calculations and Franck-Condon spectral simulations. These simulations support assigning the anion ground state as the X ̃ B 1 3 state. With this assignment, photodetachment to the B ̃ A 1 2 state involves a nominally forbidden two-electron transition, suggesting extensive configuration interaction in neutral VO 2 .

Description: Rotational diffusion processes are correlated with nanoparticle visualization and manipulation techniques, widely used in nanocomposites, nanofluids, bioscience, and so on. However, a systematical methodology of deriving this diffusivity is still lacking. In the current work, three molecular dynamics (MD) schemes, including equilibrium (Green-Kubo formula and Einstein relation) and nonequilibrium (Einstein–Smoluchowski relation) methods, are developed to calculate the rotational diffusion coefficient, taking a single rigid carbon nanotube in fluid argon as a case. We can conclude that the three methods produce same results on the basis of plenty of data with variation of the calculation parameters (tube length, diameter, fluid temperature, density, and viscosity), indicative of the validity and accuracy of the MD simulations. However, these results have a non-negligible deviation from the theoretical predictions of Tirado et al. [J. Chem. Phys.81, 2047 (1984)], which may come from several unrevealed factors of the theory. The three MD methods proposed in this paper can also be applied to other situations of calculating rotational diffusion coefficient.

Description: A general theory of electronic excitations in aggregates of molecules coupled to intramolecular vibrations and the harmonic environment is developed for simulation of the third-order nonlinear spectroscopy signals. It is applied in studies of the time-resolved two-dimensional coherent spectra of four characteristic model systems: weakly/strongly vibronically coupled molecular dimers interacting with high/low frequency intramolecular vibrations. The results allow us to (i) classify and define the typical spectroscopic features of vibronically coupled molecules, (ii) separate the cases, when the long-lived quantum coherences due to vibrational lifetime borrowing should be expected, (iii) define when the complete exciton–vibrational mixing occurs, and (iv) when separation of excitonic and vibrational coherences is possible.

Description: An analytical equation of state (EoS) is derived to describe the isotropic (I) and nematic (N) phase of linear- and partially flexible tangent hard-sphere chain fluids and their mixtures. The EoS is based on an extension of Onsager's second virial theory that was developed in our previous work [T. van Westen, B. Oyarzún, T. J. H. Vlugt, and J. Gross, J. Chem. Phys.139, 034505 (2013)]. Higher virial coefficients are calculated using a Vega-Lago rescaling procedure, which is hereby generalized to mixtures. The EoS is used to study (1) the effect of length bidispersity on the I-N and N-N phase behavior of binary linear tangent hard-sphere chain fluid mixtures, (2) the effect of partial molecular flexibility on the binary phase diagram, and (3) the solubility of hard-sphere solutes in I- and N tangent hard-sphere chain fluids. By changing the length bidispersity, two types of phase diagrams were found. The first type is characterized by an I-N region at low pressure and a N-N demixed region at higher pressure that starts from an I-N-N triphase equilibrium. The second type does not show the I-N-N equilibrium. Instead, the N-N region starts from a lower critical point at a pressure above the I-N region. The results for the I-N region are in excellent agreement with the results from molecular simulations. It is shown that the N-N demixing is driven both by orientational and configurational/excluded volume entropy. By making the chains partially flexible, it is shown that the driving force resulting from the configurational entropy is reduced (due to a less anisotropic pair-excluded volume), resulting in a shift of the N-N demixed region to higher pressure. Compared to linear chains, no topological differences in the phase diagram were found. We show that the solubility of hard-sphere solutes decreases across the I-N phase transition. Furthermore, it is shown that by using a liquid crystal mixture as the solvent, the solubility difference can by maximized by tuning the composition. Theoretical results for the Henry's law constant of the hard-sphere solute are in good agreement with the results from molecular simulation.

Description: Equilibrium and dynamical characteristics pertaining to the solvation of the fluorescent probe Coumarin 480 within different confining environments are investigated using molecular dynamics simulations. Three kinds of confining systems are examined: (i) the cetyltrimethylammonium bromide (CTAB)/isooctane/1-hexanol/water; cationic inverse micelle (IM) (ii) a CTAB/water direct micelle (DM), and (iii) a silica-surfactant nanocomposite, comprising a cylindrical silica pore (SP) containing small amounts of water and CTAB species adsorbed at the pore walls. The solvation structures in the three environments differ at a qualitative level: an exchange between bulk- and interface-like solvation states was found in the IM, whereas in the DM, the solvation states of the probe are characterized by its embedding at the interface, trapped among the surfactant heads and tails. Within the SP structure, the coumarin exhibits alternations between internal and interfacial solvation states that occur on a ∼20 ns time scale and operate via 90° rotations of its molecular plane. The solvation responses of the environment following a vertical excitation of the probe are also investigated. Solvation times resulted between 2 and 1000 times longer than those found in bulk water, with a fast-to-slow trend IM→DM→SP, which can be interpreted in terms of the solvation structures that prevail in each case.

Description: A fully alignment-resolved O 2 sticking experiment on a single domain Si(100)-(2×1) surface is presented. This provides the first experimental evidence that the reactivity of O 2 depends on both the polar and azimuthal angles of the molecular axis relative to a surface. It has been found that, in case of side-on collision, an O 2 molecule perpendicular to the dimer on Si(100) is about 40% more reactive than that parallel to the dimer. Comparison of the O 2 sticking on flat and vicinal Si(100) surfaces indicates that barrierless dissociation channels exist at the double layer step.

Description: A unified view on linear response of interacting systems utilizing multiconfigurational time-dependent Hartree methods is presented. The cases of one-particle and two-particle response operators for identical particles and up to all-system response operators for distinguishable degrees-of-freedom are considered. The working equations for systems of identical bosons and identical fermions, as well for systems of distinguishable particles are explicitly derived. These linear-response theories – applicable for discrete excitation spectra – provide numerically exact excitation energies and system's properties, when numerical convergence is achieved in the calculations.

Description: We present the synthesis and characterization of a benzodithiophene/thiophene alternating copolymer decorated with rigid, singly branched pendant side chains. We characterize exciton migration and recombination dynamics in these molecules in tetrahydrofuran solution, using a combination of static and time-resolved spectroscopies. As control experiments, we also measure electronic relaxation dynamics in isolated molecular analogues of both the side chain and polymer moieties. We employ semi-empirical and time-dependent density functional theory calculations to show that photoexcitation of the decorated copolymer using 395 nm laser pulses results in excited states primarily localized on the pendant side chains. We use ultrafast transient absorption spectroscopy to show that excitations are transferred to the polymer backbone faster than the instrumental response function, ∼250 fs.

Description: The dynamics properties of a new “next generation” model asphalt system that represents SHRP AAA-1 asphalt using larger molecules than past models is studied using molecular simulation. The system contains 72 molecules distributed over 12 molecule types that range from nonpolar branched alkanes to polar resins and asphaltenes. Molecular weights range from 290 to 890 g/mol. All-atom molecular dynamics simulations conducted at six temperatures from 298.15 to 533.15 K provide a wealth of correlation data. The modified Kohlrausch-Williams-Watts equation was regressed to reorientation time correlation functions and extrapolated to calculate average rotational relaxation times for individual molecules. The rotational relaxation rate of molecules decreased significantly with increasing size and decreasing temperature. Translational self-diffusion coefficients followed an Arrhenius dependence. Similar activation energies of ∼42 kJ/mol were found for all 12 molecules in the model system, while diffusion prefactors spanned an order of magnitude. Viscosities calculated directly at 533.15 K and estimated at lower temperatures using the Debye-Stokes-Einstein relationship were consistent with experimental data for asphalts. The product of diffusion coefficient and rotational relaxation time showed only small changes with temperature above 358.15 K, indicating rotation and translation that couple self-consistently with viscosity. At lower temperatures, rotation slowed more than diffusion.

Description: A new general model for describing intersecting multidimensional potential energy surfaces when motions of large amplitude are involved is presented. This model can be seen as an extension of the vibronic coupling models of Köppel et al. [“Multimode molecular dynamics beyond the Born-Oppenheimer approximation,” Adv. Chem. Phys.57, 59 (1984)]. In contrast to the original vibronic coupling models, here the number of diabatic states is larger than the number of adiabatic states and curvilinear coordinates are used in a systematic way. Following general considerations, the approach is applied to the fitting of the potential energy surfaces for the very complex nonadiabatic photodynamics of benzopyran. Preliminary results are presented at the complete active space self-consistent field level of theory and with up to 12 active degrees of freedom. Special emphasis is placed on the physical interpretation of the diabatic states and on the influence of the various degrees of freedom on the fit.

Description: The claim that Grote-Hynes theory (GHT), when it provides accurate rates, is equivalent to multidimensional variational transition state theory (VTST) has been debated for decades with convincing arguments on both sides. For the two theories to be equivalent a perfect dividing surface with no recrossing must exist. We describe an easily implemented test employing deterministic microcanonical (NVE) trajectories which can identify situations where no perfect dividing surface exists and thereby potentially falsify the claim of equivalence. We use this test to reach data-supported conclusions about the relationship between GHT and VTST.

Description: We report thermal rate coefficients for 12 reactions of rare gas cations (Ne + , Ar + , Kr + , Xe + ) with halide anions (Cl − , Br − , I − ), comprising both mutual neutralization (MN) and transfer ionization. No rate coefficients have been previously reported for these reactions; however, the development of the Variable Electron and Neutral Density Attachment Mass Spectrometry technique makes it possible to measure the difference of the rate coefficients for pairs of parallel reactions in a Flowing Afterglow-Langmuir Probe apparatus. Measurements of 18 such combinations of competing reaction pairs yield an over-determined data set from which a consistent set of rate coefficients of the 12 MN reactions can be deduced. Unlike rate coefficients of MN reactions involving at least one polyatomic ion, which vary by at most a factor of ∼3, those of the atom-atom reactions vary by at least a factor 60 depending on the species. It is found that the rate coefficients involving light rare-gas ions are larger than those for the heavier rare-gas ions, but the opposite trend is observed in the progression from Cl − to I − . The largest rate coefficient is 6.5 × 10 −8 cm 3  s −1 for Ne + with I − . Rate coefficients for Ar + , Kr + , and Xe + reacting with Br 2 − are also reported.

Description: Following the recent quantum dynamics investigation of the charge transfer at an oligothiophene-fullerene heterojunction by the multi-configuration time dependent Hartree method [H. Tamura, R. Martinazzo, M. Ruckenbauer and I. Burghardt, J. Chem. Phys.137, 22A540 (2012)], we revisit the transfer process by a perturbative non-Markovian master equation treated by the time local auxiliary density matrix approach. We compare the efficiency of the spin-boson model calibrated by quantum chemistry with the effective mode representation. A collective mode is extracted from the spin-boson spectral density. It is weakly coupled to a residual bath of vibrational modes, allowing second-order dynamics. The electron transfer is analyzed for a sampling of inter-fragment distances showing the fine interplay of the electronic coupling and energy gap on the relaxation. The electronic coherence, expected to play a role in the process, is preserved during about 200 fs.

Description: For discrete classical Molecular dynamics obtained by the “Verlet” algorithm (VA) with the time increment h there exists a shadow Hamiltonian H ̃ with energy E ̃ ( h ) , for which the discrete particle positions lie on the analytic trajectories for H ̃ . Here, we proof that there, independent of such an analytic analogy, exists an exact hidden energy invariance E * for VA dynamics. The fact that the discrete VA dynamics has the same invariances as Newtonian dynamics raises the question, which of the formulations that are correct, or alternatively, the most appropriate formulation of classical dynamics. In this context the relation between the discrete VA dynamics and the (general) discrete dynamics investigated by Lee [Phys. Lett. B122, 217 (1983)] is presented and discussed.

Description: The electrostatic interaction between finite charge distributions, ρ(r), in a neutralizing background is considered as an extension of the one component plasma (OCP) model of point charges. A general form for the interaction potential is obtained which can be applied to molecular theories of many simple charged fluids and mixtures and to the molecular dynamics (MD) simulation of such systems. The formalism is applied to the study of a fluid of Gaussian charges in a neutralizing background by MD simulation and using hypernetted-chain integral equation theory. The treatment of these interactions is extended to a periodic system using a Fourier Transform formulation and, for a rapidly decaying charge distribution, an application of the Ewald method. The contributions of the self-energy and neutralizing background to the system's energy are explicitly included in the formulation. Calculations reveal differences in behavior from the OCP model when the Wigner-Seitz radius is of order and less than the Gaussian charge density decay length. For certain parameter values these systems can exhibit a multiple occupancy crystalline phase at high density which undergoes re-entrant melting at higher density. An exploration of the effects of the various length scales of the system on the equation of state and radial distribution function is made.

Description: We have investigated the inhomogeneous interior of confined spherical cavities as capsules containing encapsulated binary hard sphere mixtures for different compositions and cavity wall rigidity. Such a greatly simplified model manifests the effects of macromolecular crowding arising from excluded volume interactions in a tiny cell or a cellular nucleus. By fixing the number of large particles, the level of crowding is adjusted by changing the amount of small hard spheres in the cavity. For a rigid cavity, large spheres tend to pack in liquid-like order apart from the surface to the center of the cavity as the crowding level is increased. Whereas, for a soft cavity, larger spheres tend to blend with small spheres in the peripheral region at near the boundary of the cavity, and are susceptible to be depleted from the interior of the cavity as the cavity becomes more crowded. These results may help future elucidation of the thermodynamic pathways to stabilize the inhomogeneous structure of mixtures confined in cavities, such as the derepression of genome materials around the interior rim of the nucleus in a cancerous cell.

Description: We present a Nested Markov chain Monte Carlo (NMC) scheme for building equilibrium averages based on accurate potentials such as density functional theory. Metropolis sampling of a reference system, defined by an inexpensive but approximate potential, was used to substantially decorrelate configurations at which the potential of interest was evaluated, thereby dramatically reducing the number needed to build ensemble averages at a given level of precision. The efficiency of this procedure was maximized on-the-fly through variation of the reference system thermodynamic state (characterized here by its inverse temperature β 0 ), which was otherwise unconstrained. Local density approximation results are presented for shocked states of argon at pressures from 4 to 60 GPa, where—depending on the quality of the reference system potential—acceptance probabilities were enhanced by factors of 1.2–28 relative to unoptimized NMC. The optimization procedure compensated strongly for reference potential shortcomings, as evidenced by significantly higher speedups when using a reference potential of lower quality. The efficiency of optimized NMC is shown to be competitive with that of standard ab initio molecular dynamics in the canonical ensemble.

Description: A stochastic algorithm is proposed that can compute the basis-set-incompleteness correction to the second-order many-body perturbation (MP2) energy of a polyatomic molecule. It evaluates the sum of two-, three-, and four-electron integrals over an explicit function of electron-electron distances by a Monte Carlo (MC) integration at an operation cost per MC step increasing only quadratically with size. The method can reproduce the corrections to the MP2/cc-pVTZ energies of H 2 O, CH 4 , and C 6 H 6 within a few m E h after several million MC steps. It circumvents the resolution-of-the-identity approximation to the nonfactorable three-electron integrals usually necessary in the conventional explicitly correlated (R12 or F12) methods.

Description: For classical Brownian systems driven out of equilibrium, we derive inhomogeneous two-time correlation functions from functional differentiation of the one-body density and current with respect to external fields. In order to allow for appropriate freedom upon building the derivatives, we formally supplement the Smoluchowski dynamics by a source term, which vanishes at the physical solution. These techniques are applied to obtain a complete set of dynamic Ornstein-Zernike equations, which serve for the development of approximation schemes. The rules of functional calculus lead naturally to non-Markovian equations of motion for the two-time correlators. Memory functions are identified as functional derivatives of a unique space- and time-nonlocal dissipation power functional.

Description: The static and dynamic properties of striped colloidal particles are obtained using molecular dynamics computer simulations. Striped particles with n = 2 to n = 7 stripes of alternating electric charge are modeled at a high level of detail through a pointwise (PW) representation of the particle surface. We also consider the extent to which striped particles are similar to comparable isotropically attractive particles—such as depletion attracting colloids—by modeling striped particles with an isotropic pair interaction computed by coarse-graining (CG) over orientations at a pair level. Surprisingly, the CG models reproduce the static structure of the PW models for a range of volume fractions and interaction strengths consistent with the fluid region of the phase diagram for all n . As a corollary, different n -striped particle systems with comparable pair affinities (e.g., dimer equilibrium constant) have similar static structure. Stronger pair interactions lead to a collapsed structure in simulation as consistent with a glass-like phase. Different n -striped particle systems are found to have different phase boundaries and for certain n 's no glass-like state is observed in any of our simulations. The CG model is found to have accelerated dynamics relative to the PW model for the same range of fluid conditions for which the models have identical static structure. This suggests striped electrostatic particles have slower dynamics than comparable isotropically attractive colloids. The slower dynamics result from a larger number of long-duration reversible bonds between pairs of striped particles than seen in isotropically attractive systems. We also found that higher n -striped particles systems generally have slower dynamics than lower n -striped systems with comparable pair affinities.

Description: We present a novel implementation of Kohn–Sham density-functional theory utilizing London atomic orbitals as basis functions. External magnetic fields are treated non-perturbatively, which enable the study of both magnetic response properties and the effects of strong fields, using either standard density functionals or current-density functionals—the implementation is the first fully self-consistent implementation of the latter for molecules. Pilot applications are presented for the finite-field calculation of molecular magnetizabilities, hypermagnetizabilities, and nuclear magnetic resonance shielding constants, focusing on the impact of current-density functionals on the accuracy of the results. Existing current-density functionals based on the gauge-invariant vorticity are tested and found to be sensitive to numerical details of their implementation. Furthermore, when appropriately regularized, the resulting magnetic properties show no improvement over standard density-functional results. An advantage of the present implementation is the ability to apply density-functional theory to molecules in very strong magnetic fields, where the perturbative approach breaks down. Comparison with high accuracy full-configuration-interaction results show that the inadequacies of current-density approximations are exacerbated with increasing magnetic field strength. Standard density-functionals remain well behaved but fail to deliver high accuracy. The need for improved current-dependent density-functionals, and how they may be tested using the presented implementation, is discussed in light of our findings.

Description: We present the first analytic implementation of cubic and quartic force constants at the level of Kohn–Sham density-functional theory. The implementation is based on an open-ended formalism for the evaluation of energy derivatives in an atomic-orbital basis. The implementation relies on the availability of open-ended codes for evaluation of one- and two-electron integrals differentiated with respect to nuclear displacements as well as automatic differentiation of the exchange–correlation kernels. We use generalized second-order vibrational perturbation theory to calculate the fundamental frequencies of methane, ethane, benzene, and aniline, comparing B3LYP, BLYP, and Hartree–Fock results. The Hartree–Fock anharmonic corrections agree well with the B3LYP corrections when calculated at the B3LYP geometry and from B3LYP normal coordinates, suggesting that the inclusion of electron correlation is not essential for the reliable calculation of cubic and quartic force constants.

Description: The interactions between proteins/peptides and materials are crucial to research and development in many biomedical engineering fields. The energetics of such interactions are key in the evaluation of new proteins/peptides and materials. Much research has recently focused on the quality of free energy profiles by Jarzynski's equality, a widely used equation in biosystems. In the present work, considerable discrepancies were observed between the results obtained by Jarzynski's equality and those derived by umbrella sampling in biomaterial-water model systems. Detailed analyses confirm that such discrepancies turn up only when the target molecule moves in the high-density water layer on a material surface. Then a hybrid scheme was adopted based on this observation. The agreement between the results of the hybrid scheme and umbrella sampling confirms the former observation, which indicates an approach to a fast and accurate estimation of adsorption free energy for large biomaterial interfacial systems.

Description: Following the arrival of Cassini at Titan in 2004, the Titan atmosphere has been shown to contain large complex polycyclic-aromatic hydrocarbons. Since Cassini has provided a great deal of data, there exists a need for kinetic rate data to help with modeling this atmosphere. One type of kinetic data needed is electron-ion dissociative recombination (e-IDR) rate constants. These data are not readily available for larger compounds, such as naphthalene, or oxygen containing compounds, such as 1,4 dioxane or furan. Here, the rate constants for naphthalene, 1,4 dioxane, and furan have been measured and their temperature dependencies are determined when possible, using the University of Georgia's Variable Temperature Flowing Afterglow. The rate constants are compared with those previously published for other compounds; these show trends which illustrate the effects which multi-rings and oxygen heteroatoms substitutions have upon e-IDR rate constants.

Description: Long known as a fully polarized band with a near vanishing depolarization ratio [η s = 0.05, W. Holzer and R. Ouillon, Chem. Phys. Lett.24, 589 (1974)], the 2ν 5 Raman overtone of SF 6 has so far been considered as of having a prohibitively weak anisotropic spectrum [D. P. Shelton and L. Ulivi, J. Chem. Phys.89, 149 (1988)]. Here, we report the first anisotropic spectrum of this overtone, at room temperature and for 13 gas densities ranging between 2 and 27 amagat. This spectrum is 10 times broader and 50 times weaker than the isotropic counterpart of the overtone [D. Kremer, F. Rachet, and M. Chrysos, J. Chem. Phys.138, 174308 (2013)] and its profile much more sensitive to pressure effects than the profile of the isotropic spectrum. From our measurements an accurate value for the anisotropy matrix-element |⟨000020|Δα|000000⟩| was derived and this value was found to be comparable to that of the mean-polarizability 000020 α ¯ 000000 . Among other conclusions our study offers compelling evidence that, in Raman spectroscopy, highly polarized bands or tiny depolarization ratios are not necessarily incompatible with large polarizability anisotropy transition matrix-elements. Our findings and the way to analyze them suggest that new strategies should be developed on the basis of the complementarity inherent in independent incoherent Raman experiments that run with two different incident-beam polarizations, and on concerted efforts to ab initio calculate accurate data for first and second polarizability derivatives. Values for these derivatives are still rarities in the literature of SF 6 .

Description: We observed homogeneous nucleation process of supercooled liquid Fe by molecular dynamics simulations. Using bond-orientational order parameters together with Voronoi polyhedron method, we characterized local structure, calculated the volume of Voronoi polyhedra of atoms and identified the structure and density fluctuations. We monitored the formation of nucleus and analyzed its inner structure. The birth and growth of the pre-nucleus and nucleus are accompanied with aggregating and disaggregating processes in the time scale of femtosecond. Only the initial solid-like clusters (ISLC), ranging from 1 to 7 atoms, pop up directly from liquid. The relation between the logarithm of number of clusters and the cluster size was found to be linear for ISLCs and was observed to be parabolic for all solid-like clusters (SLC) due to aggregating and disaggregating effects. The nucleus and pre-nuclei mainly consist of body centered cubic (BCC) and hexagonal close packed atoms, while the BCC atoms tend to be located at the surface. Medium-range structure fluctuations induce the birth of ISLCs, benefit the aggregation of embryos and remarkably promote the nucleation. But density fluctuations contribute little to nucleation. The lifetime of most icosahedral-like atoms (ICO) is shorter than 0.7 ps. No obvious relationship was found between structure/density fluctuations and the appearance of ICO atoms.

Description: Velocity-mapped imaging and theoretical calculations have been used to study the angular distribution of the products of NO predissociation following its excitation to the 11s , 10p , 11p , and 9f Rydberg levels based on the NO + ( X 2 Σ + ) core. The Rydberg states were reached from the NO ( A 2 Σ + , v = 0, N = 2, J = 1.5) level prepared with strong alignment by excitation with linear polarization from NO ( X 2 Π, v = 0, N = 1, J = 0.5). Ion dip spectra of the Rydberg states were recorded along with velocity-mapped images at the major peaks. The results are compared to calculations based on a previous theoretical approach modified to include transitions to states of Hund's case (d) coupling. The reasonable agreement shows the predictive value of the theory. The theory has also been used to reassess and explain previous results and to understand variations in the rate of photodissociation with components of the 10p and 11p Rydberg states.

Description: We study the conformational properties of heteropolymers containing two types of monomers A and B , modeled as self-attracting self-avoiding random walks on a regular lattice. Such a model can describe in particular the sequences of hydrophobic and hydrophilic residues in proteins [K. F. Lau and K. A. Dill, Macromolecules 22, 3986 (1989)] and polyampholytes with oppositely charged groups [Y. Kantor and M. Kardar, Europhys. Lett.28, 169 (1994)]. Treating the sequences of the two types of monomers as quenched random variables, we provide a systematic analysis of possible generalizations of this model. To this end we apply the pruned-enriched Rosenbluth chain-growth algorithm, which allows us to obtain the phase diagrams of extended and compact states coexistence as function of both the temperature and fraction of A and B monomers along the heteropolymer chain.

Description: We report a comparative study of the features in dissociative double ionization by high energy electron impact of N 2 and CO molecules. The ratio of cross-section of charge symmetric dissociative ionization to non-dissociative ionization (CSD-to-ND ratio) and the kinetic energy release (KER) spectra of dissociation are experimentally measured and carefully corrected for various ion transmission losses and detector inefficiencies. Given that the double ionization cross sections of these iso-electronic diatomics are very similar, the large difference in the CSD-to-ND ratios must be attributable to the differences in the evolution dynamics of the dications. To understand these differences, potential energy curves (PECs) of dications have been computed using multi-reference configuration interaction method. The Franck-Condon factors and tunneling life times of vibrational levels of dications have also been computed. While the KER spectrum of N 2 + + can be readily explained by considering dissociation via repulsive states and tunneling of meta-stable states, indirect dissociation processes such as predissociation and autoionization have to be taken into account to understand the major features of the KER spectrum of CO ++ . Direct and indirect processes identified on the basis of the PECs and experimental KER spectra also provide insights into the differences in the CSD-to-ND ratios.

Description: Pure rotational spectra of the allenyloxy radical (CH 2 =CCHO) were observed by Fourier transform microwave (FTMW) and FTMW-millimeter wave double-resonance spectroscopy. Molecular constants including the hyperfine interaction constants of CH 2 =CCHO in the 2 A ″ ground electronic state were precisely determined. Ab initio calculations indicate that CH 2 =CCHO has a linear C−C−C backbone with C s symmetry, where the formyl group is in the C s plane and perpendicular to the methylene group. The determined rotational constants and the inertial defect agree well with those derived from the calculations, implying that the calculated molecular structure is reasonable. The fine and hyperfine constants also agree with those derived from the calculated spin density, where the unpaired electron is located mainly on the central carbon atom. The ground state CH 2 =CCHO can, thus, be described as taking the formylvinyl (CH 2 =Ċ−CH=O) form rather than as the allenyloxy (CH 2 =C=CH− O ̇ ) form.

Description: By means of lattice-Boltzmann simulations the drag force on a sphere of radius R approaching a superhydrophobic striped wall has been investigated as a function of arbitrary separation h . Superhydrophobic (perfect-slip vs. no-slip) stripes are characterized by a texture period L and a fraction of the gas area ϕ. For very large values of h / R , we recover the macroscopic formulae for a sphere moving towards a hydrophilic no-slip plane. For h / R = O (1), the drag force is smaller than predicted by classical theories for hydrophilic no-slip surfaces, but larger than expected for a sphere interacting with a uniform perfectly slipping wall. At a thinner gap, h ≪ R the force reduction compared to a classical result becomes more pronounced, and is maximized by increasing ϕ. In the limit of very small separations, our simulation data are in quantitative agreement with an asymptotic equation, which relates a correction to a force for superhydrophobic slip to texture parameters. In addition, we examine the flow and pressure field and observe their oscillatory character in the transverse direction in the vicinity of the wall, which reflects the influence of the heterogeneity and anisotropy of the striped texture. Finally, we investigate the lateral force on the sphere, which is detectable in case of very small separations and is maximized by stripes with ϕ = 0.5.

Description: We numerically analyze the influence of hydrodynamic interactions (HI) on scale-free polymer networks under external forces. We envisage that these polymer networks move under external applied fields: mechanical or electrical type. We work in the framework of the bead-spring model and include the HI in a preaveraged Oseen manner. We focus on physical quantities which depend on the eigenvalue spectrum: averaged monomer displacement and the mechanical relaxation modulus with its two components: the storage and the loss modulus. We analyze how these quantities depend on the underlying topology, such as the number of monomers and γ, which is the parameter that measures the connectivity of a network. The influence of topology appears in the intermediate region of time or frequency. For intermediate time/frequency region we do not encounter regions of constant slope by considering the hydrodynamic interactions, differently from the Rouse model, which neglects these HI. As expected we observe a more linear-like topology for high values of γ and a more star-like shape for low γs.

Description: The first band in the electronic spectrum of pyrrole is calculated from wavepacket propagations performed using the MCTDH method. To do so, two model Hamiltonians are constructed to describe seven low-lying excited electronic states of pyrrole. These Hamiltonians are based on the vibronic coupling model, and are parameterised via fitting to extensive CASPT2 and EOM-CCSD calculations. A detailed analysis of the structure of pyrrole's electronic spectrum in the range 5.5 to 6.5 eV is made. The role of intensity borrowing from transitions to ππ * states by lower-lying 3 s and 3 p Rydberg states is assessed, and reassignments of much of the spectrum are subsequently made which indicate that most of the states in the spectrum are predominantly Rydberg in character. The resulting conclusions drawn serve to highlight the limitations of assignments based on the matching of calculated vertical excitation energies and the positions of peak maxima observed in electronic spectra.

Description: The mesoscopic dynamical properties of oil-in-water microemulsions (MEs) bridged with telechelic polymers of different number of arms and with different lengths of hydrophobic stickers were studied with neutron spin-echo (NSE) probing the dynamics in the size range of individual ME droplets. These results then were compared to those of dynamicic light scattering (DLS) which allow to investigate the dynamics on a much larger length scale. Studies were performed as a function of the polymer concentration, number of polymer arms, and length of the hydrophobic end-group. In general it is observed that the polymer bridging has a rather small influence on the local dynamics, despite the fact that the polymer addition leads to an increase of viscosity by several orders of magnitude. In contrast to results from rheology and DLS, where the dynamics on much larger length and time scales are observed, NSE shows that the linear polymer is more efficient in arresting the motion of individual ME droplets. This finding can be explained by a simple simulation, merely by the fact that the interconnection of droplets becomes more efficient with a decreasing number of arms. This means that the dynamics observed on the short and on the longer length scale depend in an opposite way on the number of arms and hydrophobic stickers.

Description: We extend the exact solutions of the Di Marzio-Rubin matrix method for the thermodynamic properties, including chain density, of a linear polymer molecule confined to walk on a lattice of finite size. Our extensions enable (a) the use of higher dimensions (explicit 2D and 3D lattices), (b) lattice boundaries of arbitrary shape, and (c) the flexibility to allow each monomer to have its own energy of attraction for each lattice site. In the case of the large chain limit, we demonstrate how periodic boundary conditions can also be employed to reduce computation time. Advantages to this method include easy definition of chemical and physical structure (or surface roughness) of the lattice and site-specific monomer-specific energetics, and straightforward relatively fast computations. We show the usefulness and ease of implementation of this extension by examining the effect of energy variation along the lattice walls of an infinite rectangular cylinder with the idea of studying the changes in properties caused by chemical inhomogeneities on the surface of the box. Herein, we look particularly at the polymer density profile as a function of temperature in the confined region for very long polymers. One particularly striking result is the shift in the critical condition for adsorption due to surface energy inhomogeneities and the length scale of the inhomogeneities; an observation that could have important implications for polymer chromatography. Our method should have applications to both copolymers and biopolymers of arbitrary molar mass.

Description: Ab initio simulations that account for nuclear quantum effects have been used to examine the order-disorder transition in squaric acid, a prototypical H-bonded antiferroelectric crystal. Our simulations reproduce the 〉100 K difference in transition temperature observed upon deuteration as well as the strong geometrical isotope effect observed on intermolecular separations within the crystal. We find that collective transfer of protons along the H-bonding chains – facilitated by quantum mechanical tunneling – is critical to the order-disorder transition and the geometrical isotope effect. This sheds light on the origin of isotope effects and the importance of tunneling in squaric acid which likely extends to other H-bonded ferroelectrics.

Description: We present a combined experimental and theoretical study of NO( v = 3 → 3, 2, 1) scattering from a Au(111) surface at incidence translational energies ranging from 0.1 to 1.2 eV. Experimentally, molecular beam–surface scattering is combined with vibrational overtone pumping and quantum-state selective detection of the recoiling molecules. Theoretically, we employ a recently developed first-principles approach, which employs an Independent Electron Surface Hopping (IESH) algorithm to model the nonadiabatic dynamics on a Newns-Anderson Hamiltonian derived from density functional theory. This approach has been successful when compared to previously reported NO/Au scattering data. The experiments presented here show that vibrational relaxation probabilities increase with incidence energy of translation. The theoretical simulations incorrectly predict high relaxation probabilities at low incidence translational energy. We show that this behavior originates from trajectories exhibiting multiple bounces at the surface, associated with deeper penetration and favored (N-down) molecular orientation, resulting in a higher average number of electronic hops and thus stronger vibrational relaxation. The experimentally observed narrow angular distributions suggest that mainly single-bounce collisions are important. Restricting the simulations by selecting only single-bounce trajectories improves agreement with experiment. The multiple bounce artifacts discovered in this work are also present in simulations employing electronic friction and even for electronically adiabatic simulations, meaning they are not a direct result of the IESH algorithm. This work demonstrates how even subtle errors in the adiabatic interaction potential, especially those that influence the interaction time of the molecule with the surface, can lead to an incorrect description of electronically nonadiabatic vibrational energy transfer in molecule-surface collisions.

Description: We present the first implementation of the relativistic quantum chemical two- and four-component density matrix renormalization group algorithm that includes a variational description of scalar-relativistic effects and spin–orbit coupling. Numerical results based on the four-component Dirac–Coulomb Hamiltonian are presented for the standard reference molecule for correlated relativistic benchmarks: thallium hydride.

Description: We elucidate the origin of large differences (two-fold or more) in the fixed-node errors between the first- vs second-row systems for single-configuration trial wave functions in quantum Monte Carlo calculations. This significant difference in the valence fixed-node biases is studied across a set of atoms, molecules, and also Si, C solid crystals. We show that the key features which affect the fixed-node errors are the differences in electron density and the degree of node nonlinearity. The findings reveal how the accuracy of the quantum Monte Carlo varies across a variety of systems, provide new perspectives on the origins of the fixed-node biases in calculations of molecular and condensed systems, and carry implications for pseudopotential constructions for heavy elements.

Description: Crystal structure prediction at high pressures unbiased by any prior known structure information has recently become a topic of considerable interest. We here present a short overview of recently developed structure prediction methods and propose current challenges for crystal structure prediction. We focus on first-principles crystal structure prediction at high pressures, paying particular attention to novel high pressure structures uncovered by efficient structure prediction methods. Finally, a brief perspective on the outstanding issues that remain to be solved and some directions for future structure prediction researches at high pressure are presented and discussed.

Description: We present quantum calculations of the relaxation matrix for the Q branch of N 2 at room temperature using a recently proposed N 2 -N 2 rigid rotor potential. Close coupling calculations were complemented by coupled states studies at high energies and provide about 10 200 two-body state-to state cross sections from which the needed one-body cross-sections may be obtained. For such temperatures, convergence has to be thoroughly analyzed since such conditions are close to the limit of current computational feasibility. This has been done using complementary calculations based on the energy corrected sudden formalism. Agreement of these quantum predictions with experimental data is good, but the main goal of this work is to provide a benchmark relaxation matrix for testing more approximate methods which remain of a great utility for complex molecular systems at room (and higher) temperatures.

Description: We develop a method that combines data mining and first principles calculation to guide the designing of distorted cubane Mn 4 + Mn 3 3 + single molecule magnets. The essential idea of the method is a process consisting of sparse regressions and cross-validation for analyzing calculated data of the materials. The method allows us to demonstrate that the exchange coupling between Mn 4 + and Mn 3 + ions can be predicted from the electronegativities of constituent ligands and the structural features of the molecule by a linear regression model with high accuracy. The relations between the structural features and magnetic properties of the materials are quantitatively and consistently evaluated and presented by a graph. We also discuss the properties of the materials and guide the material design basing on the obtained results.

Description: In this paper, based on the effective intermolecular potential with well separated density and configuration contributions and the definition of the isothermal bulk modulus, we derive two similar equations of state dedicated to describe volumetric data of supercooled liquids studied in the extremely wide pressure range related to the density range, which is extremely wide in comparison with the experimental range reached so far in pressure-volume-temperature measurements of glass-forming liquids. Both the equations comply with the generalized density scaling law of molecular dynamics versus h (ρ)/ T at different densities ρ and temperatures T , where the scaling exponent can be in general only a density function γ(ρ) = d ln  h / d ln ρ as recently argued by the theory of isomorphs. We successfully verify these equations of state by using data obtained from molecular dynamics simulations of the Kob-Andersen binary Lennard-Jones liquid. As a very important result, we find that the one-parameter density function h ( ρ ) analytically formulated in the case of this prototypical model of supercooled liquid, which implies the one-parameter density function γ ( ρ ), is able to scale the structural relaxation times with the value of this function parameter determined by fitting the volumetric simulation data to the equations of state. We also show that these equations of state properly describe the pressure dependences of the isothermal bulk modulus and the configurational isothermal bulk modulus in the extremely wide pressure range investigated by the computer simulations. Moreover, we discuss the possible forms of the density functions h ( ρ ) and γ ( ρ ) for real glass formers, which are suggested to be different from those valid for the model of supercooled liquid based on the Lennard-Jones intermolecular potential.

Description: We present a fully general derivation of the Laplace–Young formula and discuss the interplay between the intrinsic surface geometry and the extrinsic one ensuing from the immersion of the surface in the ordinary Euclidean three-dimensional space. We prove that the (reversible) work done in a general surface deformation can be expressed in terms of the surface stress tensor and the variation of the intrinsic surface metric.

Description: The behavior of Liquid N,N-dimethylformamide subjected to a wide range of externally applied electric fields (from 0.001 V/nm to 1 V/nm) has been investigated through molecular dynamics simulation. To approach the objective the AMOEBA polarizable force field was extended to include the interaction of the external electric field with atomic partial charges and the contribution to the atomic polarization. The simulation results were evaluated with quantum mechanical calculations. The results from the present force field for the liquid at normal conditions were compared with the experimental and molecular dynamics results with non-polarizable and other polarizable force fields. The uniform external electric fields of higher than 0.01 V/nm have a significant effect on the structure of the liquid, which exhibits a variation in numerous properties, including molecular polarization, local cluster structure, rotation, alignment, energetics, and bulk thermodynamic and structural properties.

Description: A systematic ab initio treatment of the nitryl halides (XNO 2 ) and the cis - and trans - conformers of the halide nitrites (XONO), where X = Cl, Br, and I, have been carried out using highly correlated methods with sequences of correlation consistent basis sets. Equilibrium geometries and harmonic frequencies have been accurately calculated in all cases at the explicitly correlated CCSD(T)-F12b level of theory, including the effects of core-valence correlation for the former. Where experimental values are available for the equilibrium structures (ClNO 2 and BrNO 2 ), the present calculations are in excellent agreement; however, the X-O distances are slightly too long by about 0.01 Å due to missing multireference effects. Accurate predictions for the iodine species are made for the first time. The vertical electronic excitation spectra have been calculated using equation-of-motion coupled cluster methods for the low-lying singlet states and multireference configuration interaction for both singlet and triplet states. The latter also included the effects of spin-orbit coupling to provide oscillator strengths for the ground state singlet to excited triplet transitions. While for ClNO 2 the transitions to excited singlet states all occur at wavelengths shorter than 310 nm, there is one longer wavelength singlet transition in BrNO 2 and two in the case of INO 2 . The long wavelength tail in the XNO 2 species is predicted to be dominated by transitions to triplet states. In addition to red-shifting from X = Cl to I, the triplet transitions also increase in oscillator strength, becoming comparable to many of the singlet transitions in the case of INO 2 . Hence in particular, the latter species should be very photolabile. Similar trends are observed and reported for the halogen nitrites, many of which for the first time.

Description: A general method to study classical scattering in n -dimension is developed. Through classical trajectory calculations, the three-body recombination is computed as a function of the collision energy for helium atoms, as an example. Quantum calculations are also performed for the J Π = 0 + symmetry of the three-body recombination rate in order to compare with the classical results, yielding good agreement for E ≳ 1 K. The classical threshold law is derived and numerically confirmed for the Newtonian three-body recombination rate. Finally, a relationship is found between the quantum and classical three-body hard hypersphere elastic cross sections which is analogous to the well-known shadow scattering in two-body collisions.

Description: F1s and C1s photoelectron angular distributions are considered for CH 3 F, a molecule which does not support any shape resonance. In spite of the absence of features in the photoionization cross section profile, the recoil frame photoelectron angular distributions (RFPADs) exhibits dramatic changes depending on both the photoelectron energy and polarization geometry. Time-dependent density functional theory calculations are also given to rationalize the photoionization dynamics. The RFPADs have been compared with the theoretical calculations, in order to assess the accuracy of the theoretical method and rationalize the experimental findings. The effect of finite acceptance angles for both ionic fragments and photoelectrons has been included in the calculations, as well as the effect of rotational averaging around the fragmentation axis. Excellent agreement between theory and experiment is obtained, confirming the good quality of the calculated dynamical quantities (dipole moments and phase shifts).

Description: The nucleation of particles from trace gases in the atmosphere is an important source of cloud condensation nuclei, and these are vital for the formation of clouds in view of the high supersaturations required for homogeneous water droplet nucleation. The methods of quantum chemistry have increasingly been employed to model nucleation due to their high accuracy and efficiency in calculating configurational energies; and nucleation rates can be obtained from the associated free energies of particle formation. However, even in such advanced approaches, it is typically assumed that the nuclei have a classical nature, which is questionable for some systems. The importance of zero-point motion (also known as quantum nuclear dynamics) in modelling small clusters of sulphuric acid and water is tested here using the path integral molecular dynamics method at the density functional level of theory. The general effect of zero-point motion is to distort the mean structure slightly, and to promote the extent of proton transfer with respect to classical behaviour. In a particular configuration of one sulphuric acid molecule with three waters, the range of positions explored by a proton between a sulphuric acid and a water molecule at 300 K (a broad range in contrast to the confinement suggested by geometry optimisation at 0 K) is clearly affected by the inclusion of zero point motion, and similar effects are observed for other configurations.

Description: An extension of the well established Reverse Monte Carlo (RMC) method for modeling systems under close confinement has been developed. The method overcomes limitations induced by close confinement in systems such as fluids adsorbed in microporous materials. As a test of the method, we investigate a model system of 36 Ar adsorbed into two zeolites with significantly different pore sizes: Silicalite-I (a pure silica form of ZSM-5 zeolite, characterized by relatively narrow channels forming a 3D network) at partial and full loadings and siliceous Faujasite (which exhibits relatively wide channels and large cavities). The model systems are simulated using grand canonical Monte Carlo and, in each case, its structure factor is used as input for the proposed method, which shows a rapid convergence and yields an adsorbate microscopic structure in good agreement with that of the model system, even to the level of three body correlations, when these are induced by the confining media. The application to experimental systems is straightforward incorporating factors such as the experimental resolution and appropriate q-sampling, along the lines of previous experiences of RMC modeling of powder diffraction data including Bragg and diffuse scattering.

Description: We use molecular dynamics simulations to investigate the arrangement of polymer chains when absorbed onto a long, single-wall carbon nano-tube (SWCNT). We study the conformation and organization of the polymer chains on the SWCNT and their dependence on the tube's diameter and the rate of cooling. We use two types of cooling processes: direct quenching and gradual cooling. The radial density distribution function and bond orientational order parameter are used to characterize the polymer chain structure near the surface. In the direct cooling process, the beads of the polymer chain organize in lamella-like patterns on the surface of the SWCNT with the long axis of the lamella parallel to the axis of the SWCNT. In a stepwise, gradual cooling process, the polymer beads form a helical pattern on the surface of a relatively thick SWCNT, but form a lamella-like pattern on the surface of a very thin SWCNT. We develop a theoretical (free energy) model to explain this difference in pattern structures for the gradual cooling process and also provide a qualitative explanation for the pattern that forms from the direct cooling process.

Description: RNA folding kinetics is directly tied to RNA biological functions. We introduce here a new approach for predicting the folding kinetics of RNA secondary structure with pseudoknots. This approach is based on our previous established helix-based method for predicting the folding kinetics of RNA secondary structure. In this approach, the transition rates for an elementary step: (1) formation, (2) disruption of a helix stem, and (3) helix formation with concomitant partial melting of an incompatible helix, are calculated with the free energy landscape. The folding kinetics of the Hepatitis delta virus (HDV) ribozyme and the mutated sequences are studied with this method. The folding pathways are identified by recursive searching the states with high net flux-in(out) population starting from the native state. The theory results are in good agreement with that of the experiments. The results indicate that the bi-phasic folding kinetics for the wt HDV sequence is ascribed to the kinetic partitioning mechanism: Part of the population will quickly fold to the native state along the fast pathway, while another part of the population will fold along the slow pathway, in which the population is trapped in a non-native state. Single mutation not only changes the folding rate but also the folding pathway.

Description: We introduce a stable and efficient complex Langevin (CL) scheme to enable the first direct numerical simulations of the coherent-states (CS) formulation of polymer field theory. In contrast with Edwards’ well-known auxiliary-field (AF) framework, the CS formulation does not contain an embedded nonlinear, non-local, implicit functional of the auxiliary fields, and the action of the field theory has a fully explicit, semi-local, and finite-order polynomial character. In the context of a polymer solution model, we demonstrate that the new CS-CL dynamical scheme for sampling fluctuations in the space of coherent states yields results in good agreement with now-standard AF-CL simulations. The formalism is potentially applicable to a broad range of polymer architectures and may facilitate systematic generation of trial actions for use in coarse-graining and numerical renormalization-group studies.

Description: Microgels are soft particles that can be deformed and compressed, which would induce intriguing phase behaviors at high packing fractions. Poly( N -isopropylacrylamide) (PNIPAM) microgels, with a lower critical solution temperature (LCST) of 33 °C, have attracted considerable interests as model colloids, since the volume of them and the interaction between the microgels can be tuned precisely by temperature. In this work, the linear viscoelastic properties of PNIPAM microgel suspensions have been investigated using mechanical spectroscopy. A particular attention is focused on the phase behaviors at high concentrations. With increasing concentration the system undergoes a repulsive glass-to-gel transition below the LCST, while, as temperature is raised across the LCST, the system undergoes a gel-to-attractive glass transition. A mechanism of these transitions for the microgels is proposed based on the directional interaction between the particles. In moderate concentration or de-swelling microgels the interaction is isotropic leading to the glass phase, while in concentrated and deformed microgels the interaction is directional leading to the gel phase. Our results enrich the current understanding of the phase transition in microgel systems and shed new light on the phase diagram of colloidal suspensions in general.

Description: We report on the shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil, as obtained from fixed-nuclei elastic scattering calculations performed with the Schwinger multichannel method with pseudopotentials. Our results are in good agreement with the available electron transmission spectroscopy data, and support the existence of three π* resonances in uracil and 5-fluorouracil. As expected, the anion states are more stable in the substituted molecules than in uracil. Since the stabilization is stronger in 5-chlorouracil, the lowest π* resonance in this system becomes a bound anion state. The present results also support the existence of a low-lying σ CCl * shape resonance in 5-chlorouracil. Exploratory calculations performed at selected C–Cl bond lengths suggest that the σ CCl * resonance could couple to the two lowest π* states, giving rise to a very rich dissociation dynamics. These facts would be compatible with the complex branching of the dissociative electron attachment cross sections, even though we cannot discuss any details of the vibration dynamics based only on the present fixed-nuclei results.

Description: The entropic forces on the self-retracting granular chains, which are confined in channels with different widths, are determined. The time dependence of the length of chain remaining in the channel L in (t) is measured. The entropic force is treated as the only parameter in fitting the solution of the nonlinear equation of motion of L in (t) to the experimental data. The dependence of the entropic force on the width of the confining channel can be expressed as a power-law with an exponent of 1.3, which is consistent with the previous theoretical predictions for the entropy loss due to confinement.

Description: A second order classical perturbation theory is developed and applied to elastic atom corrugated surface scattering. The resulting theory accounts for experimentally observed asymmetry in the final angular distributions. These include qualitative features, such as reduction of the asymmetry in the intensity of the rainbow peaks with increased incidence energy as well as the asymmetry in the location of the rainbow peaks with respect to the specular scattering angle. The theory is especially applicable to “soft” corrugated potentials. Expressions for the angular distribution are derived for the exponential repulsive and Morse potential models. The theory is implemented numerically to a simplified model of the scattering of an Ar atom from a LiF(100) surface.

Description: A kinetic model is developed to describe the heterogeneous condensation of vapor on an insoluble spherical particle. This new model considers two mechanisms of cluster growth: direct addition of water molecules from the vapor and surface diffusion of adsorbed water molecules on the particle. The effect of line tension is also included in the model. For the first time, the exact expression of evaporation coefficient is derived for heterogeneous condensation of vapor on an insoluble spherical particle by using the detailed balance. The obtained expression of evaporation coefficient is proved to be also correct in the homogeneous condensation and the heterogeneous condensation on a planar solid surface. The contributions of the two mechanisms to heterogeneous condensation including the effect of line tension are evaluated and analysed. It is found that the cluster growth via surface diffusion of adsorbed water molecules on the particle is more important than the direct addition from the vapor. As an example of our model applications, the growth rate of the cap shaped droplet on the insoluble spherical particle is derived. Our evaluation shows that the growth rate of droplet in heterogeneous condensation is larger than that in homogeneous condensation. These results indicate that an explicit kinetic model is benefit to the study of heterogeneous condensation on an insoluble spherical particle.

Description: We propose in the present paper a SCC-DFTB/FF (Self-Consistent-Charge Density Functional based Tight Binding/Force-Field) scheme adapted to the investigation of molecules trapped in rare gas environments. With respect to usual FF descriptions, the model involves the interaction of quantum electrons in a molecule with rare gas atoms in an anisotropic scheme. It includes polarization and dispersion contributions and can be used for both neutral and charged species. Parameters for this model are determined for hydrocarbon-argon complexes and the model is validated for small hydrocarbons. With the future aim of studying polycyclic aromatic hydrocarbons in Ar matrices, extensive benchmark calculations are performed on (C 6 H 6 ) +/0 Ar n clusters against DFT and CCSD(T) calculations for the smaller sizes, and more generally against other experimental and theoretical data. Results on the structures and energetics (isomer ordering and energy separation, cohesion energy per Ar atom) are presented in detail for n = 1–8, 13, 20, 27, and 30, for both neutrals and cations. We confirm that the clustering of Ar atoms leads to a monotonous decrease of the ionization potential of benzene for n ⩽ 20, in line with previous experimental and FF data.

Description: In our previous paper [B. T. Sutcliffe and R. G. Woolley, J. Chem. Phys.137, 22A544 (2012)] we argued that the Born-Oppenheimer approximation could not be based on an exact transformation of the molecular Schrödinger equation. In this Comment we suggest that the fundamental reason for the approximate nature of the Born-Oppenheimer model is the lack of a complete set of functions for the electronic space, and the need to describe the continuous spectrum using spectral projection.

Description: Ionic solutions exhibiting multiple association are described within the binding mean spherical approximation (BiMSA). This model is based on the Wertheim formalism, in the framework of the primitive model at the McMillan-Mayer level. The cation and the anion form the various complexes according to stepwise complexation-equilibria. Analytic expressions for the Helmholtz energy, the internal energy, the speciation, and for the osmotic and activity coefficients are given considering a binary solution with an arbitrary number of association sites on one type of ion (polyion) and one site on the ions of opposite sign (counterions). As an alternative, mean field expressions, as developed in SAFT-type theories, are also presented. The result obtained from the latter approximate method exhibits a reasonable agreement with those from BiMSA for the speciation, and a remarkable one for the osmotic coefficient.

Description: A theory of nonlinear response of chemical kinetics, in which multiple perturbations are used to probe the time evolution of nonlinear chemical systems, is developed. Expressions for nonlinear chemical response functions and susceptibilities, which can serve as multidimensional measures of the kinetic pathways and rates, are derived. A new class of multidimensional measures that combine multiple perturbations and measurements is also introduced. Nonlinear fluctuation-dissipation relations for steady-state chemical systems, which replace operations of concentration measurement and perturbations, are proposed. Several applications to the analysis of complex reaction mechanisms are provided.

Description: The development of efficient theoretical methods for describing electron transfer (ET) reactions in condensed phases is important for a variety of chemical and biological applications. Previously, dynamical dielectric continuum theory was used to derive Langevin equations for a single collective solvent coordinate describing ET in a polar solvent. In this theory, the parameters are directly related to the physical properties of the system and can be determined from experimental data or explicit molecular dynamics simulations. Herein, we combine these Langevin equations with surface hopping nonadiabatic dynamics methods to calculate the rate constants for thermal ET reactions in polar solvents for a wide range of electronic couplings and reaction free energies. Comparison of explicit and implicit solvent calculations illustrates that the mapping from explicit to implicit solvent models is valid even for solvents exhibiting complex relaxation behavior with multiple relaxation time scales and a short-time inertial response. The rate constants calculated for implicit solvent models with a single solvent relaxation time scale corresponding to water, acetonitrile, and methanol agree well with analytical theories in the Golden rule and solvent-controlled regimes, as well as in the intermediate regime. The implicit solvent models with two relaxation time scales are in qualitative agreement with the analytical theories but quantitatively overestimate the rate constants compared to these theories. Analysis of these simulations elucidates the importance of multiple relaxation time scales and the inertial component of the solvent response, as well as potential shortcomings of the analytical theories based on single time scale solvent relaxation models. This implicit solvent approach will enable the simulation of a wide range of ET reactions via the stochastic dynamics of a single collective solvent coordinate with parameters that are relevant to experimentally accessible systems.

Description: Underlying physical principles for the high efficiency of excitation energy transfer in light-harvesting complexes are not fully understood. Notably, the degree of robustness of these systems for transporting energy is not known considering their realistic interactions with vibrational and radiative environments within the surrounding solvent and scaffold proteins. In this work, we employ an efficient technique to estimate energy transfer efficiency of such complex excitonic systems. We observe that the dynamics of the Fenna-Matthews-Olson (FMO) complex leads to optimal and robust energy transport due to a convergence of energy scales among all important internal and external parameters. In particular, we show that the FMO energy transfer efficiency is optimum and stable with respect to important parameters of environmental interactions including reorganization energy λ, bath frequency cutoff γ, temperature T , and bath spatial correlations. We identify the ratio of k B λ T / ℏ γ ⁢ g as a single key parameter governing quantum transport efficiency, where g is the average excitonic energy gap.

Description: Decomposition of energetic material 3,4-dinitropyrazole (DNP) and two model molecules 4-nitropyrazole and 1-nitropyrazole is investigated both theoretically and experimentally. The initial decomposition mechanisms for these three nitropyrazoles are explored with complete active space self-consistent field (CASSCF) level. The NO molecule is observed as an initial decomposition product from all three materials subsequent to UV excitation. Observed NO products are rotationally cold (

Description: By modifying the tetrahedrality (the strength of the three body interactions) in the well-known Stillinger-Weber model for silicon, we study the diffusivity of a series of model liquids as a function of tetrahedrality and temperature at fixed pressure. Previous work has shown that at constant temperature, the diffusivity exhibits a maximum as a function of tetrahedrality, which we refer to as the diffusivity anomaly, in analogy with the well-known anomaly in water upon variation of pressure at constant temperature. We explore to what extent the structural and thermodynamic changes accompanying changes in the interaction potential can help rationalize the diffusivity anomaly, by employing the Rosenfeld relation between diffusivity and the excess entropy (over the ideal gas reference value), and the pair correlation entropy, which provides an approximation to the excess entropy in terms of the pair correlation function. We find that in the modified Stillinger-Weber liquids, the Rosenfeld relation works well above the melting temperatures but exhibits deviations below, with the deviations becoming smaller for smaller tetrahedrality. Further we find that both the excess entropy and the pair correlation entropy at constant temperature go through maxima as a function of the tetrahedrality, thus demonstrating the close relationship between structural, thermodynamic, and dynamical anomalies in the modified Stillinger-Weber liquids.

Description: The electronic coupling that mediates energy transfer in molecular aggregates is theoretically investigated using the principles of quantum electrodynamics (QED). In this context, both the electromagnetic tensor and rate equation relating to these couplings are re-examined with a focus on the role of the relative distance and orientation of transition dipole moment pairs, considering near-, intermediate-, and far-zone contributions to the coupling. The QED based coupling terms are investigated both analytically and numerically, and they are physically interpreted in terms of the character of the mediating (virtual) photons. The spatial dependence of the couplings for a two-dimensional molecular aggregate of ordered and isotropic transition dipole moments is numerically calculated. Further, Pauli Master Equations are employed for a one-dimensional chain of molecules and donor-acceptor pairs, to investigate the importance of intermediate- and far-zone contributions to the electronic coupling on electronic energy transfer dynamics. The results indicate that although Förster theory is often qualitatively and quantitatively correct for describing electronic energy transfer (EET) processes, intermediate- and far-zone coupling terms could sometimes be non-negligible for correctly describing EET in natural and artificial, mesoscopic, solar energy harvesting systems. In particular, the results indicate that these terms are non-negligible when using Förster resonance energy transfer spectroscopic ruler techniques for distances 〉10 nm.

Description: We present results of molecular dynamics simulations for coarse-grained polymer brushes in a wide temperature range to investigate the factors that affect the glass transition in these systems. We focus on the influences of free surface, polymer-substrate interaction strength, grafting density, and chain length not only on the change of glass transition temperature T g , but also the fragility D of the glass former. It is found that the confinement can enhance the dependence of the T g on the cooling rate as compared to the bulk melt. Our layer-resolved analysis demonstrates that it is possible to control the glass transition temperature T g of polymer brushes by tuning the polymer-substrate interaction strength, the grafting density, and the chain length. Moreover, we find quantitative differences in the influence range of the substrate and the free surface on the density and dynamics. This stresses the importance of long range cooperative motion in glass formers near the glass transition temperature. Furthermore, the string-like cooperative motion analysis demonstrates that there exists a close relation among glass transition temperature T g , fragility D , and string length ⟨ S ⟩. The polymer brushes that possess larger string length ⟨ S ⟩ tend to have relatively higher T g and smaller D . Our results suggest that confining a fragile glass former through forming polymer brushes changes not only the glass transition temperature T g , but also the very nature of relaxation process.

Description: Thermodynamic and diffusion properties of water + methane mixtures in a single liquid phase are studied using NVT molecular dynamics. An extensive comparison is reported for the thermal pressure coefficient, compressibilities, expansion coefficients, heat capacities, Joule-Thomson coefficient, zero frequency speed of sound, and diffusion coefficient at methane concentrations up to 15% in the temperature range of 298–650 K. The simulations reveal a complex concentration dependence of the thermodynamic properties of water + methane mixtures. The compressibilities, heat capacities, and diffusion coefficients decrease with increasing methane concentration, whereas values of the thermal expansion coefficients and speed of sound increase. Increasing methane concentration considerably retards the self-diffusion of both water and methane in the mixture. These effects are caused by changes in hydrogen bond network, solvation shell structure, and dynamics of water molecules induced by the solvation of methane at constant volume conditions.

Description: Quasi-two-dimensional fluids can be generated by confining a fluid between two parallel walls with narrow separation. Such fluids exhibit an inhomogeneous structure perpendicular to the walls due to the loss of translational symmetry. Taking the transversal degrees of freedom as a perturbation to an appropriate 2D reference fluid we provide a systematic expansion of the m -particle density for arbitrary m . To leading order in the slit width this density factorizes into the densities of the transversal and lateral degrees of freedom. Explicit expressions for the next-to-leading order terms are elaborated analytically quantifying the onset of inhomogeneity. The case m = 1 yields the density profile with a curvature given by an integral over the pair-distribution function of the corresponding 2D reference fluid, which reduces to its 2D contact value in the case of pure excluded-volume interactions. Interestingly, we find that the 2D limit is subtle and requires stringent conditions on the fluid-wall interactions. We quantify the rapidity of convergence for various structural quantities to their 2D counterparts.

Description: We propose an adiabatic reweighting algorithm for computing the free energy along an external parameter from adaptive molecular dynamics simulations. The adaptive bias is estimated using Bayes identity and information from all the sampled configurations. We apply the algorithm to a structural transition in a cluster and to the migration of a crystalline defect along a reaction coordinate. Compared to standard adaptive molecular dynamics, we observe an acceleration of convergence. With the aid of the algorithm, it is also possible to iteratively construct the free energy along the reaction coordinate without having to differentiate the gradient of the reaction coordinate or any biasing potential.

Description: A fluid constituted of hard spherocylinders is studied using a density functional theory for non-spherical hard particles, which can be written as a function of weighted densities. This is based on an extended deconvolution of the Mayer f -function for arbitrarily shaped convex hard bodies in tensorial weight functions, which depend each only on the shape and orientation of a single particle. In the course of an examination of the isotropic-nematic interface at coexistence the functional is applied to anisotropic and inhomogeneous problems for the first time. We find good qualitative agreement with other theoretical predictions and also with Monte Carlo simulations.

Description: We introduce a database (HAB11) of electronic coupling matrix elements ( H ab ) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n -electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute H ab values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the economic FODFT and highly cost-effective FODFTB methods underestimate couplings by 37.6% and 42.4%, respectively, due to neglect of interaction between donor and acceptor. The errors are systematic, however, and can be significantly reduced by applying a uniform scaling factor for each method. Applications to dimers outside the database, specifically rotated thiophene dimers and larger acenes up to pentacene, suggests that the same scaling procedure significantly improves the FODFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.

Description: The complex quantum Hamilton-Jacobi equation-Bohmian trajectories (CQHJE-BT) method is introduced as a synthetic trajectory method for integrating the complex quantum Hamilton-Jacobi equation for the complex action function by propagating an ensemble of real-valued correlated Bohmian trajectories. Substituting the wave function expressed in exponential form in terms of the complex action into the time-dependent Schrödinger equation yields the complex quantum Hamilton-Jacobi equation. We transform this equation into the arbitrary Lagrangian-Eulerian version with the grid velocity matching the flow velocity of the probability fluid. The resulting equation describing the rate of change in the complex action transported along Bohmian trajectories is simultaneously integrated with the guidance equation for Bohmian trajectories, and the time-dependent wave function is readily synthesized. The spatial derivatives of the complex action required for the integration scheme are obtained by solving one moving least squares matrix equation. In addition, the method is applied to the photodissociation of NOCl. The photodissociation dynamics of NOCl can be accurately described by propagating a small ensemble of trajectories. This study demonstrates that the CQHJE-BT method combines the considerable advantages of both the real and the complex quantum trajectory methods previously developed for wave packet dynamics.

Description: Multichannel quantum-defect theory (MQDT) is used to calculate the electron binding energies of np Rydberg states of H 2 , HD, and D 2 around n = 60 at an accuracy of better than 0.5 MHz. The theory includes the effects of rovibronic channel interactions and the hyperfine structure, and has been extended to the calculation of the asymmetric hyperfine structure of Rydberg states of a heteronuclear diatomic molecule (HD). Starting values for the eigenquantum-defect parameters of MQDT were extracted from ab initio potential-energy functions for the low-lying p Rydberg states of molecular hydrogen and subsequently refined in a global weighted fit to available experimental data on the singlet and triplet Rydberg states of H 2 and D 2 . The electron binding energies of high- np Rydberg states derived in this work represent important quantities for future determinations of the adiabatic ionization energies of H 2 , HD, and D 2 at sub-MHz accuracy.

Description: Interaction forces between carboxylate colloidal latex particles of about 2 μ m in diameter immersed in aqueous solutions of monovalent salts were measured with the colloidal probe technique, which is based on the atomic force microscope. We have systematically varied the ionic strength, the type of salt, and also the surface charge densities of the particles through changes in the solution pH. Based on these measurements, we have accurately measured the dispersion forces acting between the particles and estimated the apparent Hamaker constant to be (2.0 ± 0.5) × 10 −21 J at a separation distance of about 10 nm. This value is basically independent of the salt concentration and the type of salt. Good agreement with Lifshitz theory is found when roughness effects are taken into account. The combination of retardation and roughness effects reduces the value of the apparent Hamaker constant and its ionic strength dependence with respect to the case of ideally smooth surfaces.

Description: We propose a method for computing the activation barrier for chemical reactions involving molecules subjected to mechanical stress. The method avoids reactant and transition-state saddle optimizations at every force by, instead, solving the differential equations governing the force dependence of the critical points (i.e., minima and saddles) on the system's potential energy surface (PES). As a result, only zero-force geometry optimization (or, more generally, optimization performed at a single force value) is required by the method. In many cases, minima and transition-state saddles only exist within a range of forces and disappear beyond a certain critical point. Our method identifies such force-induced instabilities as points at which one of the Hessian eigenvalues vanishes. We elucidate the nature of those instabilities as fold and cusp catastrophes, where two or three critical points on the force-modified PES coalesce, and provide a classification of various physically distinct instability scenarios, each illustrated with a concrete chemical example.

Description: The equilibrium vacancy concentration in nickel was determined from ab initio calculations performed with both generalized gradient approximation and local density approximation up to the melting point. We focus the study on the vacancy formation entropy expressed as a sum of a vibration and an electronic contribution, which were determined from the vibration modes and the electronic densities of states. Applying a method based on the quasi-harmonic approximation, the temperature dependence of the defect formation energy and entropy were calculated. We show that the vibrations of the first shell of atoms around the defect are predominant to the vibration formation entropy. On the other hand, the electronic formation entropy is very sensitive to the exchange-correlation potential used for the calculations. Finally, the vacancy concentration is computed at finite temperature with the calculated values for the defect formation energy and entropy. In order to reconcile point-defects concentration obtained with our calculations and experimental data, we conducted complementary calorimetric measurements of the vacancy concentration in the 1073–1273 K temperature range. Close agreement between theory and experiments at high temperature is achieved if the calculations are performed with the generalized gradient approximation and both vibration and electronic contributions to the formation entropy are taken into account.

Description: Combining the reduced hierarchy equations of motion (HEOM) approach with the Wigner-function formalism, we investigate nonperturbatively exciton dissociation under the influence of a phonon bath in an organic heterojunction. The exciton is modeled by an electron-hole pair with the electron moving in the presence of both an external electric field and the Coulomb attraction potential from the hole. In the absence of a phonon bath, calculated HEOM results reproduce those from the Onsager-Braun theory in weak electric fields. In the presence of a phonon bath, substantial deviations from the Onsager-Braun theory are found, signaling phonon-induced quantum effects. Furthermore, time evolution of the spatial current distribution is examined, and an initial spike followed by a polarity change of the transient photocurrent have been recovered.

Description: We report a hybrid mesophase consisting of magnetic nanorods confined between the non-ionic surfactant bilayers of a lamellar phase. The magnetic field-induced ordering of the nanorods was measured experimentally and modeled by a two-dimensional Onsager theory including the third virial coefficient. The nanorods are strongly confined in layers, with no orientational coupling from one layer to the next. At high volume concentration they exhibit spontaneous in-plane orientational ordering and form a stack of independent two-dimensional nematic systems. This isotropic-nematic transition is first-order.

Description: Literature values for the energy gap of long one-dimensional carbon chains vary from as little as 0.2 eV to more than 4 eV. To resolve this discrepancy, we use the GW many-body approach to calculate the band gap E g of an infinite carbon chain. We also compute the energy dependence of the attenuation coefficient β governing the decay with chain length of the electrical conductance of long chains and compare this with recent experimental measurements of the single-molecule conductance of end-capped carbon chains. For long chains, we find E g = 2.16 eV and an upper bound for β of 0.21 Å −1 .

Description: The effect of nano-sized fillers on the equilibrium and dynamical properties of a linear polymer is studied by using off-lattice Monte Carlo simulation. Fillers are arranged periodically in the system with period d and Lennard-Jones interaction between polymer and fillers is considered. Results show that the statistical dimension and dynamical diffusion of the polymer are dependent on the polymer-filler interaction strength ɛ pf and the relative size between R G0 and d , here R G0 is the radius of gyration of polymer in dilute solution. Normal diffusion of polymer is always observed in the regime 2 R G0 〉 d . And the diffusion coefficient D is scaled with chain length N as D ∼ N −α , where the exponent α increases with ɛ pf . Whereas in the regime 2 R G0 〈 d ≪ Nl 0 with l 0 the mean bond length of polymer, normal diffusion is observed only at ɛ pf 〈 2, but the polymer will be adsorbed on the fillers and cannot diffuse at ɛ pf 〉 2. In addition, we find that there is a critical interaction strength ɛ * pf = 2 in our model system.

Description: A new implementation of the approximate coupled cluster singles and doubles CC2 linear response model is reported. It employs a Cholesky decomposition of the two-electron integrals that significantly reduces the computational cost and the storage requirements of the method compared to standard implementations. Our algorithm also exploits a partitioning form of the CC2 equations which reduces the dimension of the problem and avoids the storage of doubles amplitudes. We present calculation of excitation energies of benzene using a hierarchy of basis sets and compare the results with conventional CC2 calculations. The reduction of the scaling is evaluated as well as the effect of the Cholesky decomposition parameter on the quality of the results. The new algorithm is used to perform an extrapolation to complete basis set investigation on the spectroscopically interesting benzylallene conformers. A set of calculations on medium-sized molecules is carried out to check the dependence of the accuracy of the results on the decomposition thresholds. Moreover, CC2 singlet excitation energies of the free base porphin are also presented.

Description: We combine explicit correlation via the canonical transcorrelation approach with the density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods to compute a near-exact beryllium dimer curve, without the use of composite methods. In particular, our direct density matrix renormalization group calculations produce a well-depth of D e = 931.2 cm −1 which agrees very well with recent experimentally derived estimates D e = 929.7±2 cm −1 [J. M. Merritt, V. E. Bondybey, and M. C. Heaven, Science324, 1548 (2009)] and D e = 934.6 cm −1 [K. Patkowski, V. Špirko, and K. Szalewicz, Science326, 1382 (2009)], as well the best composite theoretical estimates, D e = 938±15 cm −1 [K. Patkowski, R. Podeszwa, and K. Szalewicz, J. Phys. Chem. A111, 12822 (2007)] and D e =935.1±10 cm −1 [J. Koput, Phys. Chem. Chem. Phys.13, 20311 (2011)]. Our results suggest possible inaccuracies in the functional form of the potential used at shorter bond lengths to fit the experimental data [J. M. Merritt, V. E. Bondybey, and M. C. Heaven, Science324, 1548 (2009)]. With the density matrix renormalization group we also compute near-exact vertical excitation energies at the equilibrium geometry. These provide non-trivial benchmarks for quantum chemical methods for excited states, and illustrate the surprisingly large error that remains for 1 Σ g − 1 state with approximate multi-reference configuration interaction and equation-of-motion coupled cluster methods. Overall, we demonstrate that explicitly correlated density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods allow us to fully converge to the basis set and correlation limit of the non-relativistic Schrödinger equation in small molecules.

Description: We report high-resolution anion photoelectron spectra of the cryogenically cooled cyanomethide anion, CH 2 CN − , and its isotopologue, CD 2 CN − , using slow photoelectron velocity-map imaging (SEVI) spectroscopy. Electron affinities of 12 468(2) cm −1 for CH 2 CN and 12 402(2) cm −1 for CD 2 CN are obtained, demonstrating greater precision than previous experiments. New vibrational structure is resolved for both neutral species, especially activity of the ν 5 hydrogen umbrella modes. The ν 6 out-of-plane bending mode fundamental frequency is measured for the first time in both systems and found to be 420(10) cm −1 for CH 2 CN and 389(8) cm −1 for CD 2 CN. Some rotational structure is resolved, allowing for accurate extraction of vibrational frequencies. Temperature-dependent SEVI spectra show marked effects ascribed to controlled population of low-lying anion vibrational levels. We directly measure the inversion splitting between the first two vibrational levels of the anion ν 5 umbrella mode in both species, finding a splitting of 130(20) cm −1 for CH 2 CN − and 81(20) cm −1 for CD 2 CN − . Franck-Condon forbidden activity is observed and attributed to mode-specific vibrational autodetachment from the CH 2 CN − and CD 2 CN − dipole bound excited states. We also refine the binding energy of the anion dipole bound states to 39 and 42 cm −1 , respectively, for CH 2 CN − and CD 2 CN − .

Description: A series of nine commonly used density functional methods were assessed to accurately predict the oxidation potential of the (C 2 H 2 S 2 −2 /C 2 H 2 S 2 •− ) redox couple. It was found that due to their greater tendency for charge delocalization the GGA functionals predict a structure where the radical electron is delocalized within the alkene backbone of C 2 H 2 S 2 •− , whereas the hybrid functionals and the reference QCISD/cc-pVTZ predict that the radical electron remains localized on the sulfurs. However, chemical intuition suggests that the results obtained with the GGA functionals should be correct. Indeed, with the use of the geometries obtained at the HCTH/6-311++G(3df,3pd) level of theory both the QCISD and hybrid DFT methods yield a molecule with a delocalized electron. Notably, this new molecule lies at least 53 kJ mol −1 lower in energy than the previously optimized one that had a localized radical. Using these new structures the calculated oxidation potential was found to be 2.71–2.97 V for the nine DFT functionals tested. The M06-L functional provided the best agreement with the QCISD/cc-pVTZ reference oxidation potential of 3.28 V.

Description: We develop a classical mapping approach suitable to describe vibrationally coupled charge transport in molecular junctions based on the Cartesian mapping for many-electron systems [B. Li and W. H. Miller, J. Chem. Phys.137, 154107 (2012)]. To properly describe vibrational quantum effects in the transport characteristics, we introduce a simple transformation rewriting the Hamiltonian in terms of occupation numbers and use a binning function to facilitate quantization. The approach provides accurate results for the nonequilibrium Holstein model for a range of bias voltages, vibrational frequencies, and temperatures. It also captures the hallmarks of vibrational quantum effects apparent in step-like structure in the current-voltage characteristics at low temperatures as well as the phenomenon of Franck-Condon blockade.

Description: We have developed a method for calculating the cubic anharmonic couplings in molecular crystals for normal modes with the zero wave vector in the framework of classical mechanics, and have applied it, combined with perturbation theory, to obtain the linewidths of all infrared absorption lines of crystalline pentaerythritol tetranitrate in the terahertz region (

Description: The refinement of the Robert-Bonamy (RB) formalism by considering the line coupling for isotropic Raman Q lines of linear molecules developed in our previous study [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys.139, 034305 (2013)] has been extended to infrared P and R lines. In these calculations, the main task is to derive diagonal and off-diagonal matrix elements of the Liouville operator iS 1 − S 2 introduced in the formalism. When one considers the line coupling for isotropic Raman Q lines where their initial and final rotational quantum numbers are identical, the derivations of off-diagonal elements do not require extra correlation functions of the S ̂ operator and their Fourier transforms except for those used in deriving diagonal elements. In contrast, the derivations for infrared P and R lines become more difficult because they require a lot of new correlation functions and their Fourier transforms. By introducing two dimensional correlation functions labeled by two tensor ranks and making variable changes to become even functions, the derivations only require the latters’ two dimensional Fourier transforms evaluated at two modulation frequencies characterizing the averaged energy gap and the frequency detuning between the two coupled transitions. With the coordinate representation, it is easy to accurately derive these two dimensional correlation functions. Meanwhile, by using the sampling theory one is able to effectively evaluate their two dimensional Fourier transforms. Thus, the obstacles in considering the line coupling for P and R lines have been overcome. Numerical calculations have been carried out for the half-widths of both the isotropic Raman Q lines and the infrared P and R lines of C 2 H 2 broadened by N 2 . In comparison with values derived from the RB formalism, new calculated values are significantly reduced and become closer to measurements.

Description: We investigated the effect of the exchange-correlation kernels of Dobson and Wang (DW) [Phys. Rev. B62, 10038 (2000)] and Corradini, Del Sole, Onida, and Palummo (CDOP) [Phys. Rev. B57, 14569 (1998)] in the framework of the adiabatic connection fluctuation-dissipation theorem. The original CDOP kernel was generalized to treat inhomogeneous systems, and an efficient numerical implementation was developed. We found that both kernels improve the correlation energy in bulk silicon as compared to that evaluated from the random phase approximation (RPA). In particular, the correlation energy from the CDOP kernel is in excellent agreement with the diffusion Monte Carlo result. In the case of the Kr dimer, while the DW kernel leads to stronger binding than RPA, the CDOP kernel does the opposite. The cause of this quite different behavior of the two kernels is discussed. Our study suggests that special attention needs to be paid to describe the effective interaction at the low density regions when developing model exchange-correlation kernels.

Description: We simulate the linear and nonlinear rheology of two different viscoelastic polymer solutions, a polyisobutylene solution in pristane and an aqueous solution of hydroxypropylcellulose, using a highly coarse-grained approach known as Responsive Particle Dynamics (RaPiD) model. In RaPiD, each polymer has originally been depicted as a spherical particle with the effects of the eliminated degrees of freedom accounted for by an appropriate free energy and transient pairwise forces. Motivated by the inability of this spherical particle representation to entirely capture the nonlinear rheology of both fluids, we extended the RaPiD model by introducing a deformable particle capable of elongation. A Finite-Extensible Non-Linear Elastic potential provides a free energy penalty for particle elongation. Upon disentangling, this deformability allows more time for particles to re-entangle with neighbouring particles. We show this process to be integral towards recovering the experimental nonlinear rheology, obtaining excellent agreement. We show that the nonlinear rheology is crucially dependent upon the maximum elongation and less so on the elasticity of the particles. In addition, the description of the linear rheology has been retained in the process.

Description: Using forced Rayleigh scattering, the tracer diffusion of methyl red through water and eight other solvents at different temperatures is investigated and the results are compared to the Cohen-Turnbull theory of free-volume diffusion. In seven solvents the effective non-Arrhenius activation energy measured experimentally agrees with the Cohen-Turnbull energy. In water, however, the diffusion can be described mathematically by the free volume model but there is a disagreement of more than an order of magnitude between these energies. We propose that the unique “zero point” free volume forced onto water by the strong hydrogen bonding requires a different mechanism for tracer motion.

Description: The concentration induced shape transitions of linear model surfactants, H x T y , on a lattice have been studied using Monte Carlo simulation. It has been found that a sphere to cylinder shape transition is generally found on shortening the hydrophilic part of the surfactant and anticipates an eventual phase transition. Asymmetric surfactants with longer heads than tails ( x 〉 y ) prefer to form only spherical micelles independent of total surfactant concentration while asymmetric surfactants with longer tails than heads ( x 〈 y ) form spherical micelles at lower concentration and undergo a shape transition to cylindrical micelles on increasing the total concentration. Finally, in the case of symmetric surfactants with x = y , only the shortest surfactants H 1 T 1 and H 2 T 2 undergo a sphere to cylinder shape transition on increasing surfactant concentration. Longer symmetric surfactants are always found to prefer to form spherical micelles.