Summary of 17 chemicals evaluated by OECD TG229 using Japanese Medaka, Oryzias latipes in EXTEND 2016.

In June 2016, the Ministry of the Environment of Japan announced a program "EXTEND2016" on the implementation of testing and assessment for endocrine active chemicals, consisting of a two-tiered strategy. The aim of the Tier 1 screening and the Tier 2 testing is to identify the impacts on the endocrine system and to characterize the adverse effects to aquatic animals by endocrine disrupting chemicals detected in the aquatic environment in Japan. For the consistent assessment of the effects on reproduction associated with estrogenic, anti-estrogenic, androgenic, and/or anti-androgenic activities of chemicals throughout Tier 1 screening to Tier 2 testing, a unified test species, Japanese medaka (Oryzias latipes), has been used. For Tier 1 screening, the in vivo Fish Short-Term Reproduction Assay (OECD test guideline No. 229) was conducted for 17 chemicals that were nominated based on the results of environmental monitoring, existing knowledge obtained from a literature survey, and positive results in reporter gene assays using the estrogen receptor of Japanese medaka. In the 17 assays using Japanese medaka, adverse effects on reproduction (i.e., reduction in fecundity and/or fertility) were suggested for 10 chemicals, and a significant increase of hepatic vitellogenin in males, indicating estrogenic (estrogen receptor agonistic) potency, was found for eight chemicals at the concentrations in which no overt toxicity was observed. Based on these results, and the frequency and the concentrations detected in the Japanese environment, estrone, 4-nonylphenol (branched isomers), 4-tert-octylphenol, triphenyl phosphate, and bisphenol A were considered as high priority candidate substances for the Tier 2 testing.

Scaling up electronic structure calculations on quantum computers: The frozen natural orbital based method of increments.

The method of increments and frozen natural orbital (MI-FNO) framework is introduced to help expedite the application of noisy, intermediate-scale quantum (NISQ) devices for quantum chemistry simulations. The MI-FNO framework provides a systematic reduction of the occupied and virtual orbital spaces for quantum chemistry simulations. The correlation energies of the resulting increments from the MI-FNO reduction can then be solved by various algorithms, including quantum algorithms such as the phase estimation algorithm and the variational quantum eigensolver (VQE). The unitary coupled-cluster singles and doubles VQE framework is used to obtain correlation energies for the case of small molecules (i.e., BeH2, CH4, NH3, H2O, and HF) using the cc-pVDZ basis set. The quantum resource requirements are estimated for a constrained geometry complex catalyst that is utilized in industrial settings for the polymerization of α-olefins. We show that the MI-FNO approach provides a significant reduction in the quantum bit (qubit) requirements relative to the full system simulations. We propose that the MI-FNO framework can create scalable examples of quantum chemistry problems that are appropriate for assessing the progress of NISQ devices.

Summary of reference chemicals evaluated by the fish short-term reproduction assay, OECD TG229, using Japanese Medaka, Oryzias latipes.

Under the Organisation for Economic Co-operation and Development (OECD), the Ministry of the Environment of Japan (MOE) added Japanese medaka (Oryzias latipes) to the test guideline fish short-term reproduction assay (FSTRA) developed by the United States Environmental Protection Agency (US EPA) using fathead minnow (Pimephales promelas). The FSTRA was designed to detect endocrine disrupting effects of chemicals interacting with the hypothalamic-pituitary-gonadal axis (HPG axis) such as agonists or antagonists on the estrogen receptor (Esr) and/or the androgen receptor (AR) and steroidogenesis inhibitors. We conducted the FSTRA with Japanese medaka, in accordance with OECD test guideline number 229 (TG229), for 16 chemicals including four Esr agonists, two Esr antagonists, three AR agonists, two AR antagonists, two steroidogenesis inhibitors, two progesterone receptor agonists, and a negative substance, and evaluated the usability and the validity of the FSTRA (TG229) protocol. In addition, in vitro reporter gene assays (RGAs) using Esr1 and ARß of Japanese medaka were performed for the 16 chemicals, to support the interpretation of the in vivo effects observed in the FSTRA. In the present study, all the test chemicals, except an antiandrogenic chemical and a weak Esr agonist, significantly reduced the reproductive status of the test fish, that is, fecundity or fertility, at concentrations where no overt toxicity was observed. Moreover, vitellogenin (VTG) induction in males and formation of secondary sex characteristics (SSC), papillary processes on the anal fin, in females was sensitive endpoints to Esr and AR agonistic effects, respectively, and might be indicators of the effect concentrations in long-term exposure. Overall, it is suggested that the in vivo FSTRA supported by in vitro RGA data can adequately detect effects on the test fish, O. latipes, and probably identify the mode of action (MOA) of the chemicals tested.

A Simple Model for Relative Energies of All Fullerenes Reveals the Interplay between Intrinsic Resonance and Structural Deformation Effects in Medium-Sized Fullerenes.

Fullerenes are sheets of sp2 carbon atoms wrapped around to form spheres. With this simple consideration, we have in the present study devised and (with over 3600 DFT data points) successfully validated a simple model, termed R+D, for estimating the relative energies of fullerenes. This model contains a resonance component to account for the intrinsic differences between the π-energies of different fullerenes, and a deformation component for treating the distortions from planarity. Notably, we find that both terms (and they alone) are required to obtain good relative energies, which lends support to the formulation of the R+D model. An interesting finding is that for some medium-sized IPR fullerenes, their isomers show similar variations in the two components. We deduce that these fullerenes may represent a good opportunity for tuning molecular properties for practical applications. We hope that the promising results of the present study will encourage further investigations into fullerenes from a fundamental perspective.

A path integral molecular dynamics study on intermolecular hydrogen bond of acetic acid-arsenic acid anion and acetic acid-phosphoric acid anion clusters.

We apply ab initio path integral molecular dynamics simulation employing ωB97XD as the quantum chemical calculation method to acetic acid-arsenic acid anion and acetic acid-phosphoric acid anion clusters to investigate the difference of the hydrogen bond structure and its fluctuation such as proton transfer. We found that the nuclear quantum effect enhanced the fluctuation of the hydrogen bond structure and proton transfer, which shows treatment of the nuclear quantum effect was essential to investigate these systems. The hydrogen bond in acetic acid-arsenic acid anion cluster showed characters related to low-barrier hydrogen bonds, while acetic acid-phosphoric acid anion cluster did not. We found non-negligible distinction between these two systems, which could not be found in conventional calculations. We suggest that the difference in amount of atomic charge of the atoms consisting the hydrogen bond is the origin of the difference between acetic acid-arsenic acid and acetic acid-phosphoric acid anion cluster. © 2018 Wiley Periodicals, Inc.

The reHISS Three-Range Exchange Functional with an Optimal Variation of Hartree-Fock and Its Use in the reHISSB-D Density Functional Theory Method.

In the present study, we have reparametrized the HISS exchange functional. The new "reHISS" exchange provides a balance between short- and mid-range Hartree-Fock exchange (HFX) and a large total HFX coverage, with a fast convergence to zero HFX in the long range. The five parameters in this functional (according to equations 3 and 4 in the main text) are cSR = 0.15, cMR = 2.5279, cLR = 0, ωSR = 0.27, and ωLR = 0.2192. The combination of reHISS exchange with a reparametrized B97c-type correlation functional (Chan et al., J. Comput. Chem. 2017, 38, 2307) and a D2 dispersion term (s6 = 0.6) gives the reHISSB-D method. We find it to be more accurate than related screened-exchange methods and, importantly, its accuracy is more uniform across different properties. Fundamentally, our analysis suggests that the good performance of the reHISS exchange is related to it capturing a near-optimal proportion of HFX in the range of interelectronic distance that is important for many chemical properties, and we propose this range to be approximately 1-4Å. © 2018 Wiley Periodicals, Inc.

Solvent Dependence of Double Proton Transfer in the Formic Acid-Formamidine Complex: Path Integral Molecular Dynamics Investigation.

Solvent dependence of double proton transfer in the formic acid-formamidine (FA-FN) complex at room temperature was investigated by means of ab initio path integral molecular dynamics (AIPIMD) simulation with taking nuclear quantum and thermal effects into account. The conductor-like screening model (COSMO) was applied for solvent effect. In comparison with gas phase, double proton delocalization between two heavy atoms (O and N) in FA-FN were observed with reduced proton transfer barrier height in low dielectric constant medium (<4.8). For dielectric constant medium at 4.8, the chance of finding these two protons are more pronounced due to the solvent effect which completely washes out the proton transfer barrier. In the case of higher dielectric constant medium (>4.8), the ionic species becomes more stable than the neutral ones and the formate anion and formamidium cation are thermodynamically stable. For ab initio molecular dynamics simulation, in low dielectric constant medium (<4.8) a reduction of proton transfer barrier with solvent effect is found to be less pronounced than the AIPIMD due to the absence of nuclear quantum effect. Moreover, the motions of FA-FN complex are significantly different with increasing dielectric constant medium. Such a difference is revealed in detail by the principal component analysis.

Solvatochromism and preferential solvation of Brooker's merocyanine in water-methanol mixtures.

The excitation energy of Brooker's merocyanine in water-methanol mixtures shows nonlinear behavior with respect to the mole fraction of methanol, and it was suggested that this behavior is related to preferential solvation by methanol. We investigated the origin of this behavior and its relation to preferential solvation using the three-dimensional reference interaction site model self-consistent field method and time-dependent density functional theory. The calculated excitation energies were in good agreement with the experimental behavior. Analysis of the coordination numbers revealed preferential solvation by methanol. The free energy component analysis implied that solvent reorganization and solvation entropy drive the preferential solvation by methanol, while the direct solute-solvent interaction promotes solvation by water. The difference in the preferential solvation effect on the ground and excited states causes the nonlinear excitation energy shift. © 2017 Wiley Periodicals, Inc.

Correlation functional in screened-exchange density functional theory procedures.

In the present study, we have explored several prospects for the further development of screened-exchange density functional theory (SX-DFT) procedures. Using the performance of HSE06 as our measure, we find that the use of alternative correlation functionals (as oppose to PBEc in HSE06) also yields adequate results for a diverse set of thermochemical properties. We have further examined the performance of new SX-DFT procedures (termed HSEB-type methods) that comprise the HSEx exchange and a (near-optimal) reparametrized B97c (cOS,0 = cSS,0 = 1, cOS,1 = -1.5, cOS,2 = -0.644, cSS,1 = -0.5, and cSS,2 = 1.10) correlation functionals. The different variants of HSEB all perform comparably to or slightly better than the original HSE-type procedures. These results, together with our fundamental analysis of correlation functionals, point toward various directions for advancing SX-DFT methods. © 2017 Wiley Periodicals, Inc.

Assessment of range-separated functionals in the presence of implicit solvent: Computation of oxidation energy, reduction energy, and orbital energy.

Recently, we have investigated the ionization potential (IP) theorem for some small molecules in the presence of external electric field [M. P. Borpuzari et al., J. Chem. Phys. 144, 164113 (2016)]. In this article, we assess the performance of some density functionals, local density approximation, generalized-gradient approximation (GGA), hybrid, meta-GGA hybrid, and range-separated functionals in the presence of two different solvent dielectrics, water and cyclohexane, in reproducing the vertical oxidation energy, reduction energy, and the frontier orbital energies. We also study the accessibility of different computational solvent models like the polarizable continuum model (PCM) and non-equilibrium PCM (NEPCM) in reproducing the desired properties. In general, the range-separated functionals do not perform well in reproducing orbital energies in the PCM. Range separation with the NEPCM is better. It is found that CAM-B3LYP, M06-2X, and ωB97XD functionals reproduce highest occupied molecular orbital energy in solvents, which may be due to the cancellation of PCM and density functional theory errors. Finally, we have tested the validity of the IP theorem in the solvent environment.

Singularity Correction for Long-Range-Corrected Density Functional Theory with Plane-Wave Basis Sets.

We introduced two methods to correct the singularity in the calculation of long-range Hartree-Fock (HF) exchange for long-range-corrected density functional theory (LC-DFT) calculations in plane-wave basis sets. The first method introduces an auxiliary function to cancel out the singularity. The second method introduces a truncated long-range Coulomb potential, which has no singularity. We assessed the introduced methods using the LC-BLYP functional by applying it to isolated systems of naphthalene and pyridine. We first compared the total energies and the HOMO energies of the singularity-corrected and uncorrected calculations and confirmed that singularity correction is essential for LC-DFT calculations using plane-wave basis sets. The LC-DFT calculation results converged rapidly with respect to the cell size as the other functionals, and their results were in good agreement with the calculated results obtained using Gaussian basis sets. LC-DFT succeeded in obtaining accurate orbital energies and excitation energies. We next applied LC-DFT with singularity correction methods to the electronic structure calculations of the extended systems, Si and SiC. We confirmed that singularity correction is important for calculations of extended systems as well. The calculation results of the valence and conduction bands by LC-BLYP showed good convergence with respect to the number of k points sampled. The introduced methods succeeded in overcoming the singularity problem in HF exchange calculation. We investigated the effect of the singularity correction on the excitation state calculation and found that careful treatment of the singularities is required compared to ground-state calculations. We finally examined the excitonic effect on the band gap of the extended systems. We calculated the excitation energies to the first excited state of the extended systems using a supercell model at the Γ point and found that the excitonic binding energy, supposed to be small for inorganic semiconductors, was quite large. Our findings suggest that more investigation on the effect of the excitonic binding energy on band gaps is necessary.

Recommended approaches to the scientific evaluation of ecotoxicological hazards and risks of endocrine-active substances.

A SETAC Pellston Workshop® "Environmental Hazard and Risk Assessment Approaches for Endocrine-Active Substances (EHRA)" was held in February 2016 in Pensacola, Florida, USA. The primary objective of the workshop was to provide advice, based on current scientific understanding, to regulators and policy makers; the aim being to make considered, informed decisions on whether to select an ecotoxicological hazard- or a risk-based approach for regulating a given endocrine-disrupting substance (EDS) under review. The workshop additionally considered recent developments in the identification of EDS. Case studies were undertaken on 6 endocrine-active substances (EAS-not necessarily proven EDS, but substances known to interact directly with the endocrine system) that are representative of a range of perturbations of the endocrine system and considered to be data rich in relevant information at multiple biological levels of organization for 1 or more ecologically relevant taxa. The substances selected were 17α-ethinylestradiol, perchlorate, propiconazole, 17ß-trenbolone, tributyltin, and vinclozolin. The 6 case studies were not comprehensive safety evaluations but provided foundations for clarifying key issues and procedures that should be considered when assessing the ecotoxicological hazards and risks of EAS and EDS. The workshop also highlighted areas of scientific uncertainty, and made specific recommendations for research and methods-development to resolve some of the identified issues. The present paper provides broad guidance for scientists in regulatory authorities, industry, and academia on issues likely to arise during the ecotoxicological hazard and risk assessment of EAS and EDS. The primary conclusion of this paper, and of the SETAC Pellston Workshop on which it is based, is that if data on environmental exposure, effects on sensitive species and life-stages, delayed effects, and effects at low concentrations are robust, initiating environmental risk assessment of EDS is scientifically sound and sufficiently reliable and protective of the environment. In the absence of such data, assessment on the basis of hazard is scientifically justified until such time as relevant new information is available. Integr Environ Assess Manag 2017;13:267-279. © 2017 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals, Inc. on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Challenges in assigning endocrine-specific modes of action: Recommendations for researchers and regulators.

As regulatory programs evaluate substances for their endocrine-disrupting properties, careful study design and data interpretation are needed to distinguish between responses that are truly endocrine specific and those that are not. This is particularly important in regulatory environments where criteria are under development to identify endocrine-disrupting properties to enable hazard-based regulation. Irrespective of these processes, most jurisdictions use the World Health Organization/International Programme on Chemical Safety definition of an endocrine disruptor, requiring that a substance is demonstrated to cause a change in endocrine function that consequently leads to an adverse effect in an intact organism. Such a definition is broad, and at its most cautious might capture many general mechanisms that would not specifically denote an endocrine disruptor. In addition, endocrine responses may be adaptive in nature, designed to maintain homeostasis rather than induce an irreversible adverse effect. The likelihood of indirect effects is increased in (eco)toxicological studies that require the use of maximum tolerated concentrations or doses, which must produce some adverse effect. The misidentification of indirect effects as truly endocrine mediated has serious consequences for prompting animal- and resource-intensive testing and regulatory consequences. To minimize the risk for misidentification, an objective and transparent weight-of-evidence procedure based on biological plausibility, essentiality, and empirical evidence of key events in an adverse outcome pathway is recommended to describe the modes of action that may be involved in toxic responses in nontarget organisms. Confounding factors such as systemic toxicity, general stress, and infection can add complexity to such an evaluation and should be considered in the weight of evidence. A recommended set of questions is proffered to help guide researchers and regulators in discerning endocrine and nonendocrine responses. Although many examples provided in this study are based on ecotoxicology, the majority of the concepts and processes are applicable to both environmental and human health assessments. Integr Environ Assess Manag 2017;13:280-292. © 2016 SETAC.

Performance of the OP correlation functional in relation to its formulation: Influence of the exchange component and the effect of incorporating same-spin correlations.

In the present study, we have investigated two significant features of the OP correlation functional, namely the incorporation of the exchange functional into itself, and the inclusion of only opposite-spin (OS) effects. To explore the latter feature, we have compared OP with B95 and a new functional introduced in the present study - the OPB method that combines OP with the same-spin (SS) component of B95. In general, we find that B95 and OPB perform comparably. Our comparisons of the various DFT procedures suggest that the incorporation of a meta-GGA (e.g., TPSS) into OP and OPB does not necessarily lead to a chemically more accurate procedure than the use of a related GGA (e.g., PBE). An important finding is the more notable (and somewhat more consistent) improvement in performance with the incorporation of SS correlation, particularly for longer-range chemical properties. Nonetheless, on average across our test sets of over 800 systems, the difference between the performances of OP versus B95 or OPB is not exceedingly large. By drawing a parallel between these DFT methods and the wavefunction scaled-MP2-type methods, we reason that one can further develop the OP functional, and perhaps a wider range of correlation functionals by combining it with the technique of range separation.

From C60 to Infinity: Large-Scale Quantum Chemistry Calculations of the Heats of Formation of Higher Fullerenes.

We have carried out large-scale computational quantum chemistry calculations on the K computer to obtain heats of formation for C60 and some higher fullerenes with the DSD-PBE-PBE/cc-pVQZ double-hybrid density functional theory method. Our best estimated values are 2520.0 ± 20.7 (C60), 2683.4 ± 17.7 (C70), 2862.0 ± 18.5 (C76), 2878.8 ± 13.3 (C78), 2946.4 ± 14.5 (C84), 3067.3 ± 15.4 (C90), 3156.6 ± 16.2 (C96), 3967.7 ± 33.4 (C180), 4364 (C240) and 5415 (C320) kJ mol(-1). In our assessment, we also find that the B3-PW91-D3BJ and BMK-D3(BJ) functionals perform reasonably well. Using the convergence behavior for the calculated per-atom heats of formation, we obtained the formula ΔfH per carbon = 722n(-0.72) + 5.2 kJ mol(-1) (n = the number of carbon atoms), which enables an estimation of ΔfH for higher fullerenes more generally. A slow convergence to the graphene limit is observed, which we attribute to the relatively small proportion of fullerene carbons that are in "low-strain" regions. We further propose that it would take tens, if not hundreds, of thousands of carbons for a fullerene to roughly approach the limit. Such a distinction may be a contributing factor to the discrete properties between the two types of nanomaterials. During the course of our study, we also observe a fairly reliable means for the theoretical calculation of heats of formation for medium-sized fullerenes. This involves the use of isodesmic-type reactions with fullerenes of similar sizes to provide a good balance of the chemistry and to minimize the use of accompanying species.

Theoretical vibrational spectra of OH(-)(H2O)2: the effect of quantum distribution and vibrational coupling.

We performed ab initio path integral molecular dynamics simulations for the hydroxide-water cluster, OH(-)(H2O)2, at 50 K, 100 K, and 150 K to investigate its flexible structure. From our simulations, we found that nuclear quantum effects enhance hydroxide hydrogen atom inversion and the conformational change between isomers occurs by simultaneous rotation of the free hydrogen atom. We propose the importance of including the transition state conformer with C2 symmetry, for the description of this system at temperatures realized in predissociation experiments. Temperature dependence of relative populations of each conformer along with multidimensional vibrational calculations were used to simulate the vibrational spectra and compare with the experimental spectra of Johnson and coworkers. We assign the doublet peaks seen in the experiment at 2500 to 3000 cm(-1), as the mixture of the ionic hydrogen bonded OH stretching overtone, ionic hydrogen bonded OH bending overtone, and the combination band of the ionic hydrogen bonded OH stretch and bend, which are modulated by the van der Waals OO vibrations. We concluded that for OH(-)(H2O)2, the vibrational couplings between the ionic hydrogen bonded motion and floppy modes contribute to the broadening of peaks observed in the 2500 to 3000 cm(-1) region.

Theoretical study on the stability of double-decker type metal phthalocyanines, M(Pc)2 and M(Pc)2(+) (M = Ti, Sn and Sc): a critical assessment on the performance of density functionals.

We report the results of theoretical calculations on the optimized structures and relative energies between the D4d and D2 symmetry structures for double-decker type phthalocyanine compounds, Ti(Pc)2, Ti(Pc)2(+), Sn(Pc)2, Sn(Pc)2(+), Sc(Pc)2 and Sc(Pc)2(+), using eighteen types of functionals: B3LYP, B3PW91, B3P86, PBE1PBE, BHandHLYP, BPW91, BP86, M06, M06-2x, M06-HF, M06L, LC-BPW91, LC-ωPBE, CAM-B3LYP, B97D, ωB97, ωB97X and ωB97XD. Two phthalocyanine moieties are stacked in a face-to-face configuration in the D4d structure, but they are stapled by two σ-bonds in the D2 one. We found that the molecular symmetry of M(Pc)2 and M(Pc)2(+) depends on the balance between stabilization due to electron delocalization and exchange repulsion of π-electrons in the phthalocyanine moieties. We assessed the performance of the well-established functionals to select the appropriate functional for calculations on M(Pc)2 and M(Pc)2(+), and several important aspects came out. Generally, the hybrid GGA and hybrid meta-GGA functionals with 20-27% of the HF exchange term would give the molecular structures consistent with the experimental expectations for the double-decker type phthalocyanine compounds. Pure GGA and pure meta-GGA functionals (BPW91, BP86, M06L and B97D) have the tendency to overestimate the stability of the D4d structure. On the other hand, functionals including HF exchange for 50% and over or including long-range corrections (BHandHLYP, M06, M06-2x, M06-HF, LC-BPW91, LC-ωPBE, CAM-B3LYP, ωB97, ωB97X and ωB97XD) tend to overestimate the stability of the D2 structure. It should be emphasized that the B3LYP functional, one of the most commonly used hybrid GGA functionals with 20% HF exchange, cannot estimate the relative stability between the two molecular structures of Ti(Pc)2 appropriately. The calculation for the systems considered in this article required well-balanced treatment of the HF exchange with the accompanied exchange-correlation functional. Thus, as has been pointed out rigorously and frequently, the selection of the functional is a crucial point for reliability of the calculations.

Toward the complete range separation of non-hybrid exchange-correlation functional.

In this study, we use a very simple scheme to achieve range separation of a total exchange-correlation functional. We have utilized this methodology to combine a short-range pure density functional theory (DFT) functional with a corresponding long-range pure DFT, leading to a "Range-separated eXchange-Correlation" (RXC) scheme. By examining the performance of a range of standard exchange-correlation functionals for prototypical short- and long-range properties, we have chosen B-LYP as the short-range functional and PBE-B95 as the long-range counterpart. The results of our testing using a more diverse range of data sets show that, for properties that we deem to be short-range in nature, the performance of this prescribed RXC-DFT protocol does resemble that of B-LYP in most cases, and vice versa. Thus, this RXC-DFT protocol already provides meaningful numerical results. Furthermore, we envisage that the general RXC scheme can be easily implemented in computational chemistry software packages. This study paves a way for further refinement of such a range-separation technique for the development of better performing DFT procedures.

Ab Initio Path Integral Molecular Dynamics Study of the Nuclear Quantum Effect on Out-of-Plane Ring Deformation of Hydrogen Maleate Anion.

Ab initio path integral molecular dynamics (PIMD) simulation was performed to understand the nuclear quantum effect on the out-of-plane ring deformation of hydrogen maleate anion and investigate the existence of a stable structure with ring deformation, which was suggested in experimental observation (Fillaux et al., Chem. Phys. 1999, 120, 387-403). The isotope effect and the temperature effect are studied as well. We first investigated the nuclear quantum effect on the proton transfer. In static calculation and classical ab initio molecular dynamics simulations, the proton in the hydrogen bond is localized to either oxygen atom. On the other hand, the proton is located at the center of two oxygen atoms in quantum ab initio PIMD simulations. The nuclear quantum effect washes out the barrier of proton transfer. We next examined the nuclear quantum effect on the motion of hydrogen maleate anion. Principal component analysis revealed that the out-of-plane ring bending modes have dominant contribution to the entire molecular motion. In quantum ab initio PIMD simulations, structures with ring deformation were the global minimum for the deuterated isotope at 300 K. We analyzed the out-of-plane ring bending mode further and found that there are three minima along a ring distortion mode. We successfully found a stable structure with ring deformation of hydrogen maleate for the first time, to our knowledge, using theoretical calculation. The structures with ring deformation found in quantum simulation of the deuterated isotope allowed the proton transfer to occur more frequently than the planar structure. Static ab initio electronic structure calculation found that the structures with ring deformation have very small proton transfer barrier compared to the planar structure. We suggest that the "proton transfer driven" mechanism is the origin of stabilization for the structure with out-of-plane ring deformation.

Accurate Prediction of Hyperfine Coupling Constants in Muoniated and Hydrogenated Ethyl Radicals: Ab Initio Path Integral Simulation Study with Density Functional Theory Method.

We performed ab initio path integral molecular dynamics (PIMD) simulations with a density functional theory (DFT) method to accurately predict hyperfine coupling constants (HFCCs) in the ethyl radical (CßH3-CαH2) and its Mu-substituted (muoniated) compound (CßH2Mu-CαH2). The substitution of a Mu atom, an ultralight isotope of the H atom, with larger nuclear quantum effect is expected to strongly affect the nature of the ethyl radical. The static conventional DFT calculations of CßH3-CαH2 find that the elongation of one Cß-H bond causes a change in the shape of potential energy curve along the rotational angle via the imbalance of attractive and repulsive interactions between the methyl and methylene groups. Investigation of the methyl-group behavior including the nuclear quantum and thermal effects shows that an unbalanced CßH2Mu group with the elongated Cß-Mu bond rotates around the Cß-Cα bond in a muoniated ethyl radical, quite differently from the CßH3 group with the three equivalent Cß-H bonds in the ethyl radical. These rotations couple with other molecular motions such as the methylene-group rocking motion (inversion), leading to difficulties in reproducing the corresponding barrier heights. Our PIMD simulations successfully predict the barrier heights to be close to the experimental values and provide a significant improvement in muon and proton HFCCs given by the static conventional DFT method. Further investigation reveals that the Cß-Mu/H stretching motion, methyl-group rotation, methylene-group rocking motion, and HFCC values deeply intertwine with each other. Because these motions are different between the radicals, a proper description of the structural fluctuations reflecting the nuclear quantum and thermal effects is vital to evaluate HFCC values in theory to be comparable to the experimental ones. Accordingly, a fundamental difference in HFCC between the radicals arises from their intrinsic molecular motions at a finite temperature, in particular the methyl-group behavior.