Quantum chemical investigations on molecular clusters

Towards a converged strategy for including microsolvation in reaction mechanism calculations

A major part of chemical conversions is carried out in the fluid phase, where an accurate modeling of the involved reactions requires to also take into account solvation effects. Implicit solvation models often cover these effects with sufficient accuracy but can fail drastically when specific solvent-solute interactions are important. In those cases, microsolvation, i.e., the explicit inclusion of one or more solvent molecules, is a commonly used strategy. Nevertheless, microsolvation also introduces new challenges-a consistent workflow as well as strategies how to systematically improve prediction performance are not evident. For the COSMO and COSMO-RS solvation models, this work proposes a simple protocol to decide if microsolvation is needed and how the corresponding molecular model has to look like. To demonstrate the improved accuracy of the approach, specific application examples are presented and discussed, i.e., the computation of aqueous pK_a values and a mechanistic study of the methanol mediated Morita-Baylis-Hillman reaction.

Infrared spectra of PEHA molecule and its resistance to oxidation in water and methanol media at 298.15 K: solvent cluster size dependency

The present work investigates the infrared spectra and solvation free energies (SE) of PEHA ((E)-2-(Pyridin-2-yl) ethyl 3-(3,4-dihydroxyphenyl) acrylate) and their impact on the oxidation. The latter was examined through the ionization potential parameter (IP). These investigations were carried out by the DFT method at B3LYP/6-31G(d) for optimization and frequency calculations and corrected for BSSE. X3LYP/6-311++G(2d,2p) was employed for single-point energy calculations. Water and methanol cluster sizes were used for solvation through the explicit solvent model. Thus, the infrared spectra show that the overview frequencies of PEHA compare well with the experimental results. The intense infrared absorptions of complexes are due to the stretching of O-H bonds of solvent clusters in the range 2600-3850 cm^-1. The binding energy per solvent molecule of complexes was calculated and shows that water and methanol clusters mimic the liquid state as from 5 to 10 solvent molecules. The SE of PEHA increases with the increase of the cluster size of water and methanol in the direction of the limit. The latter was censured by the solvation done using the combined explicit-implicit solvent model. As for IP parameter, the results are largely above the IP limit and lower than the IP from gas phase. Thus, water and methanol media have an effect of lowering the IP of PEHA compound. Consequently, both media favour the oxidation of PEHA molecule, which facilitates its metabolism in human organism.

Direct and Reliable Method for Estimating the Hydrogen Bond Energies and Cooperativity in Water Clusters, Wn, n = 3 to 8

No direct method for estimating the individual O-H···O hydrogen bond (H-bond) energies in water clusters (W_n) exists in the literature. In this work, we propose such a direct method based on the molecular tailoring approach, which also enables the estimation of the cooperativity contributions. The calculated H-bond energies at MP2(full)/aug-cc-pVTZ and CCSD(T)/aug-cc-pVDZ levels for W_n, n = 3 to 8, agree well with one another and fall between 0.3 and 11.6 kcal mol^-1 with the cooperativity contributions in the range of -1.2 and 7.0 kcal mol^-1. For gauging the accuracy of our H-bond energies for a cluster, the H-bond energy sum is added to the sum of monomer energies, and the results are compared with the respective total energy. These two values agree with each other to within 8.3 mH (∼5 kcal mol^-1), testifying the accuracy of our estimated H-bond energies. Further, these H-bond strengths show a good correlation with the respective O-H stretching frequencies and the molecular electron density values at the (3, -1) O-H···O H-bond critical point.

Fragment-based quantum mechanical approach to biomolecules, molecular clusters, molecular crystals and liquids

To study large molecular systems beyond the system size that the current state-of-the-art ab initio electronic structure methods could handle, fragment-based quantum mechanical (QM) approaches have been developed over the past years, and proved to be efficient in dealing with large molecular systems at various ab initio levels. According to the fragmentation approach, a large molecular system can be divided into subsystems (fragments), and subsequently the property of the whole system can be approximately obtained by taking a proper combination of the corresponding terms of individual fragments. Therefore, the standard QM calculation of a large system could be circumvented by carrying out a series of calculations on small fragments, which significantly promotes computational efficiency. The electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method is one of the fragment-based QM approaches which has been developed by our research group in recent years. This Perspective presents the theoretical framework of this fragmentation method and its applications in biomolecules, molecular clusters, molecular crystals and liquids, including total energy calculation, protein-ligand/protein binding affinity prediction, geometry optimization, vibrational spectrum simulation, ab initio molecular dynamics simulation, and prediction of excited-state properties.

Microhydration of a Carbonyl Group: How does the Molecular Electrostatic Potential (MESP) Impact the Formation of (H2O)n:(R2C═O)Complexes?

The presence of a carbonyl group in a molecule usually leads to the identification of a π-hole on the molecular electrostatic potential (MESP) of the species. How does this electrophilic site influence the formation of microhydrated complexes? To address this point, a panel of R₂CO solutes with various MESPs was selected, and we identified the structures and properties of several complexes containing one, two, three and six water molecules. The following solutes were considered in the present study: H₂CO, F₂CO, Cl₂CO,(NC)₂CO and H₂C═CO. Geometry optimizations and frequency calculations were carried out at the LC-ωPBE/6-311++G(d,p) level, with the GD3BJ empirical correction for dispersion. For a number of n water molecules around the R₂CO solute, the structure and the features of the most stable (H₂O)_n:(R₂CO) complexes are highly dependent on the MESP of the isolated R₂CO solute. The formation of pi-hole bondings appears to play a decisive role in the initiation of a three-dimensional organization of water molecules around the solute.

Local Acid Strength of Solutions and Its Quantitative Evaluation Using Excess Infrared Nitrile Probes

We propose the concept of local acidity in condensed-phase chemistry in this work. The feature is demonstrated in trifluoroethanol (TFE) by employing two Fourier-transform infrared spectroscopy (FTIR) nitrile probes, acetonitrile (CH₃CN) and benzonitrile (PhCN). Specifically, three positive excess peaks were found in the binary systems composed of TFE and a probe using excess spectroscopy. To characterize the local acidity quantitatively, we have tried to correlate the wavenumbers of the positive excess peaks of the probes and the pK_a values in water of a series of XH-containing compounds (X = O, N, and C). Good linear relationships were discovered. Accordingly, three different pK_a values of TFE were determined based on the three positive excess infrared peaks, which are attributed to the monomer, dimer, and trimer of TFE with the help of quantum-chemical calculations. The concept of local acidity and its quantitative evaluation enrich our knowledge of acid-base chemistry and will shed light on a better understanding of microstructures of solutions.

Large-Sized Ammonia Clusters and Solvation Energies of the Proton in Ammonia

The absolute solvation energies (free energies and enthalpies) of the proton in ammonia are used to compute the pK_a of species embedded in ammonia. They are also used to compute the solvation energies of other ions in ammonia. Despite their importance, it is not possible to determine experimentally the solvation energies of the proton in a given solvent. We propose in this work a direct approach to compute the solvation energies of the proton in ammonia from large-sized neutral and protonated ammonia clusters. To undertake this investigation, we performed a geometry optimization of neutral and protonated ammonia 30-mer, 40-mer, and 50 mer to locate stable structures. These structures have been fully optimized at both APFD/6-31++g(d,p) and M06-2X/6-31++g(d,p) levels of theory. An infrared spectroscopic study of these structures has been provided to assess the reliability of our investigation. Using these structures, we have computed the absolute solvation free energy and the absolute solvation enthalpy of the proton in ammonia. It comes out that the absolute solvation free energy of the proton in ammonia is calculated to be -1192 kJ mol^-1 , whereas the absolute solvation enthalpy is evaluated to be -1214 kJ mol^-1 . © 2019 Wiley Periodicals, Inc.

Exploration of the potential energy surfaces of small ethanol clusters

The potential energy surfaces of small ethanol clusters, from dimer to pentamer, have been thoroughly explored using two different levels of theory. There is a clear relative energy gap between cyclic, linear and branched cyclic structures.

Integrating Machine Learning with the Multilayer Energy-Based Fragment Method for Excited States of Large Systems

In this work we have combined machine learning techniques with our recently developed multilayer energy-based fragment method for studying excited states of large systems. The photochemically active and inert regions are separately treated with the complete active space self-consistent field method and the trained models. This method is demonstrated to provide accurate energies and gradients leading to essentially the same excited-state potential energy surfaces and nonadiabatic dynamics compared with full ab initio results. Furthermore, in conjunction with the use of machine learning models, this method is highly parallel and exhibits low-scaling computational cost. Finally, the present work could encourage the marriage of machine learning with fragment-based electronic structure methods to explore photochemistry of large systems.

Relationships between NMR shifts and interaction energies in biphenyls, alkanes, aza-alkanes, and oxa-alkanes with X─H… H─Y and X─H… Z (X, Y = C or N; Z = N or O) hydrogen bonding

Hydrogen-hydrogen C─H^… H─C bonding between the bay-area hydrogens in biphenyls, and more generally in congested alkanes, very strained polycyclic alkanes, and cis-2-butene, has been investigated by calculation of proton nuclear magnetic resonance (NMR) shifts and atom-atom interaction energies. Computed NMR shifts for all protons in the biphenyl derivatives correlate very well with experimental data, with zero intercept, unit slope, and a root mean square deviation of 0.06 ppm. For some congested alkanes, there is generally good agreement between computed values for a selected conformer and the experimental data, when it is available. In both cases, the shift of a given proton or pair of protons tends to increase with the corresponding interaction energy. Computed NMR shift differences for methylene protons in polycyclic alkanes, where one is involved in a very short contact ("in") and the other is not ("out"), show a rough correlation with the corresponding C─H^… H─C exchange energies. The "in" and "in,in" isomers of selected aza- and diaza-cycloalkanes, respectively, are X─H^… H─N hydrogen bonded, whereas the "out" and "in,out" isomers display X─H^… N hydrogen bonds (X = C or N). Oxa-alkanes and the "in" isomers of aza-oxa-alkanes are X─H^… O hydrogen bonded. There is a very good general correlation, including both N─H^… H─Y (Y = C or N) and N─H^… Z (Z = N or O) interactions, for NH proton shifts against the exchange energy. For "in" CH protons, the data for the different C─H^… H─Y and C─H^… Z interactions are much more dispersed and the overall shift/exchange energy correlation is less satisfactory.

A multi-layer energy-based fragment method for excited states and nonadiabatic dynamics

We developed a multi-layer energy-based fragment (MLEBF) method within the many-body energy expansion framework. It supplies accurate energies and gradients, and accurately reproduces excited-state topological structures. Moreover, MLEBF-based nonadiabatic dynamics simulations give nearly the same results compared with full ab initio ones. The present work could stimulate developing energy-based fragment methods for photochemistry of large systems.

Electrostatics-Assisted Building-Up Procedure for Capturing Energy Minima of Metal Clusters: Test Case of Agn Clusters

Global geometry optimization of metal clusters is an important problem in nanophysics. The starting geometries of the clusters generated with empirical or other model potentials are generally optimized further by density functional theory (DFT)-based energy minimization. For this purpose, several algorithms such as simulated annealing, genetic algorithms, basin hopping, etc. are used. Our building-up procedure generates putative lower-energy structures of metal (M) clusters, M_n+1, M_n+2, etc., by anchoring one or more metal atoms in the vicinity of the minima of the molecular electrostatic potential (MESP) of M_n. Here, we report an application of this method to Ag_n clusters, for 5 ≤ n ≤ 20, followed up by DFT-based geometry optimization, generating several lower-energy structures than those reported in the literature. New low-energy isomers are obtained by applying the same procedure to the test case of mixed-metal clusters, Ni_nAg_m, for n + m = 4 and 5. In conclusion, our MESP-based building-up procedure offers a new general methodology for generating lower-energy geometries of metal clusters.

Molecular Hydration of Carbonic Acid: Ab Initio Quantum Chemical and Density Functional Theory Investigation

Molecular hydration of carbonic acid (H₂CO₃) is investigated in terms of bonding patterns in H₂CO₃···(H₂O) _n [ n = 1-4] hydrogen-bonded clusters within ab initio quantum chemical and density functional theory (DFT) frameworks. Successive addition of water molecules to H₂CO₃···H₂O entails elongation of O-H (hydroxyl) bond as well as contraction of specific intermolecular hydrogen bonds signifying hydration of carbonic acid; these structural features get markedly enhanced under the continuum solvation framework. A comparison between the structurally similar clusters H₂CO₃···(H₂O) _n and HCOOH···(H₂O) _n [ n = 1-3] brings out the structural stability of the former. The present investigation in conjunction with the binding energy behavior of approaching water molecule(s) should serve as a precursor for pathways exploring aqueous dissociation of H₂CO₃ for larger clusters, as well as development of force-field potentials for acid dissociation process.

An integrated instrument of DUV-IR photoionization mass spectrometry and spectroscopy for neutral clusters

We have developed an integrated instrument combining deep ultraviolet laser ionization mass spectrometry (DUV-LIMS) and infrared multiphoton dissociation (IR-MPD) spectroscopy, abbreviated as DUV-IR. The 177.3 nm DUV laser (7 eV single-photon energy) has short pulse duration (15 ps) and appropriate pulse energy (∼20 µJ), which is found to be highly efficient for low-fragment photoionization of neutral metal clusters and molecules. A home-made cluster source is designed with an adjustable formation channel suitable for the generation of different cluster series. The well-aligned components of the reflection time-of-flight mass spectrometer, as well as the coaxial design of DUV laser and molecular beam, bring forth high sensitivity and high resolution of the DUV-LIMS. Taking these advantages, well-resolved neutral V_n (n = 1-43) and (Benzene)_n (n = 1-25) clusters have been generated free of fragmentation. In addition to the generation and detection of neutral clusters, a fast-flow reaction tube is also designed downstream of the cluster source allowing to study their reactivity. In particular, a broad-range tunable IR laser (1.3-16 µm) is coupled with the DUV laser to attain IR-MPD spectroscopic analysis. This integrated system offers a general protocol to prepare various clusters to study their gas-phase reactivity and to determine their structures.

Exploring Reaction Energy Profiles Using the Molecules-in-Molecules Fragmentation-Based Approach

The Molecules-in-Molecules (MIM) fragmentation-based approach has been successfully used in previous studies to obtain the energies, optimized geometries, and spectroscopic properties of large molecular systems. The present work delineates a protocol to study the potential energy profiles for multistep chemical reactions using the MIM methodology. In a complex multistep chemical reaction, the fragmentation scheme needs to be changed as the reacting species transition into a new reaction step, resulting in a discontinuity in the potential energy curve of the reaction. In our approach, the fragmentation scheme for a particular step in a reaction is chosen on the basis of the nature of the bonding changes associated with that step. Thus, the reactant, transition state, and product are treated consistently throughout the reaction step, leading to an accurate energy barrier for that step. The discontinuity now occurs in describing the energies of reaction intermediates at the transition point between two reaction steps that are treated by two different fragmentation schemes. To address this issue, we propose a systematic procedure for obtaining continuous potential energy curves that are least shifted from their initial positions. The corrected MIM potential energy curves are continuous with activation energies preserved. Following this approach, energy profiles of complex reactions involving large molecular species can be obtained at high levels of theory with a reasonable computational cost.

Pragmatic Many-Body Approach for Economic MP2 Energy Estimation of Molecular Clusters

We propose a procedure, within the many-body analysis (MBA) framework, for an economic yet accurate estimation of the correlated method-based energies of large molecular clusters employing Dunning's augmented basis sets. The basis of the procedure is to segregate the Hartree-Fock ( E_HF) and correlation energy ( E_C) estimations. E_HF is found to differ by tens of millihartrees (mH) from its full-calculation (FC) counterpart on truncating the MBA expansion at the two-body (MBA-2) level. On the contrary, E_C is estimated with smaller error on modest hardware with limited computation time at the (MBA-2) level. In view of this, we adopt a pragmatic method wherein the E_HF (accurate to five decimal places) is taken from the FC, whereas E_C is estimated at the MBA-2 level. This method is applied to a variety of medium to large molecular clusters at the MP2 level. Preliminary results at the CCSD(T) level for (H₂O)₁₆ and (H₂O)₁₇ are also reported with tremendous savings in wall-clock time and resources. The typical errors in MP2 and CCSD(T) energies per monomer are up to 0.1 and 0.2 mH, respectively. Thus the present method, balancing accuracy and computational economy, opens a way for estimating energies of large molecular clusters using correlated theories with large basis sets.

C-H/π and C-H-O Interactions in Concert: A Study of the Anisole-Methane Complex using Resonant Ionization and Velocity Mapped Ion Imaging

Noncovalent forces such as hydrogen bonding, halogen bonding, π-π stacking, and C-H/π and C-H/O interactions hold the key to such chemical processes as protein folding, molecular self-assembly, and drug-substrate interactions. Invaluable insight into the nature and strength of these forces continues to come from the study of isolated molecular clusters. In this work, we report on a study of the isolated anisole-methane complex, where both C-H/π and C-H/O interactions are possible, using a combination of theory and experiments that include mass-selected two-color resonant two-photon ionization spectroscopy, two-color appearance potential (2CAP) measurements, and velocity mapped ion imaging (VMI). Using 2CAP and VMI, we derive the binding energies of the complex in ground, excited, and cation radical states. The experimental values from the two methods are in excellent agreement, and they are compared with selected theoretical values calculated using density functional theory and ab initio methods. The optimized ground-state cluster geometry, which is consistent with the experimental observations, shows methane sitting above the ring, interacting with anisole via both C-H/π and C-H/O interactions, and this dual mode of interaction is reflected in a larger ground-state binding energy as compared with the prototypical benzene-methane system.

Structures, relative stability and binding energies of neutral water clusters, (H2O)2–30

We have revised the structures of neutral water clusters, (H₂O)_n=2–30, with the affordable M06-2X functional, presenting up to 25 isomers for each cluster size.

Dynamic simulation of liquid molecular nanoclusters: structure, stability and quantification of internal (pseudo)symmetries

In a few hours on a standard laptop, AA-CLP MD correctly reproduces the thermodynamic properties of bulk liquids and provides information on the nanoscale dynamics of liquid nanoclusters.

Experimental and theoretical investigation of intramolecular cooperativity in cyclic benzene trimer motif

A series of new symmetrical tripodal molecules 1a–4b with a central benzene scaffold substituted with methyl/ethyl groups and three benzimidazolyl units having a bithiophene/biphenyl/5-alkylthiophene motif at the 2-position via a –CH₂– unit were synthesized and characterized by elemental analysis, HR-MS, and NMR spectroscopy. NMR spectral data reveal that all molecules adopt a cyclic benzene trimer (CBT) using three benzimidazolyl units. Intramolecular cooperative edge-to-face C–H⋯π interactions stabilize the CBT motif in solution and are strong in ethyl substituted molecules (1b–4b) compared to methyl substituted (1a–4a) ones. However, the strength of the CBT unit in the tripodal molecule is independent of the length of the substituent at the 2-position of the benzimidazolyl unit. The relative ¹H NMR chemical shift calculated at the MPW1PW91/6-311+G(d,p) level of theory corroborates the experimental values, and the calculations predict the distribution of the structures into syn isomers. The relative change in the NMR chemical shift is justified by the relative change in the magnitude of the (3,+3) critical point (CP) in the molecular electrostatic potential (MESP) topography. Also, a linear correlation of the intramolecular C–H⋯π interactions evaluated at M062X/6-311+G(d,p) with the relative NMR chemical shift suggest the latter as a measure of intramolecular cooperativity.

A family of biaryl/alkylthiophene (R–R) benzimidazolyl-based tripodal molecules with cyclic benzene trimer (CBT) motif was synthesized and studied by NMR spectroscopy and MPW1PW91/6-311+G(d,p) theory.

Evidence of Reactivity in the Membrane for the Unstable Monochloramine during MIMS Analysis

Membrane Inlet Mass Spectrometry (MIMS) was used to analyze monochloramine solutions (NH₂Cl) and ammonia solutions in a compact FTICR. Chemical ionization enables identification and quantification of the products present in the permeate. The responses of protonated monochloramine and ammonium increase linearly with the solution concentration. The enrichments were respectively 1.2 and 5.5. Pervaporation is dependent on pH and only the basic form of ammonia NH₃ pervaporates through the membrane. Unexpectedly, the small ammonia molecule permeated very slowly. It could be due to interactions with water molecules inside the membrane that create clusters. Moreover, NH₂Cl solutions, in addition to the NH₃Cl⁺ signal, presented a strong NH₄⁺ signal at m/z 18.034. Ammonia presence in the low-pressure zone before ionization is probable as NH₄⁺ was detected with all the precursors used, particularly CF₃⁺ and trimethylbenzene that presents a proton affinity higher than monochloramine. Ammonia may be formed inside the membrane due to the fact that NH₂Cl is unstable and may react with the water present in the membrane. Those results highlight the need for caution when dealing with chloramines in MIMS and more generally with unstable molecules.

UV-Vis, Fluorescence, and Resonance Raman Spectroscopic and Density Functional Theoretical Studies on 3-Amino-1,2,4-triazole: Microsolvation and Solvent-Dependent Nonadiabatic Excited State Decay in Solution

The microsolvation and photophysics of 3-amino-1,2,4-triazole (3AT) after excitation to the light-absorbing S₂(nπ*) state were studied by using resonance Raman spectroscopy and single component artificial force-induced reaction (SC-AFIR) in a global reaction route mapping (GRRM) strategy. The vibrational spectra were assigned on the basis of experimental data and density functional theory (DFT) calculations. The resonance Raman spectra of 3AT were measured to probe the excited state structural dynamics in the Franck-Condon region. The conformations of 3AT(CH₃CN)₁, 3AT(CH₃OH)₂, and 3AT(H₂O)₂ clusters were determined by combining vibrational spectrum experiments and B3LYP/6-311++G(d,p) computations. DFT calculations were carried out to obtain the minimal excitation energies of the lower-lying singlet excited states, and the curve-crossing points. It was revealed that the short-time structural dynamics of 3AT were dominated by the N-N stretching coordinates. An excited state decay mechanism is proposed: 3AT is initially excited to the S₂(nπ*) state, then the conical intersection (CI) of the S₂(nπ*)/S₁(ππ*) potential energy surfaces is crossed, and 3AT then decays to the lower solvent-dependent excited state S₁(ππ*). It subsequently returns to the S₀ state, accompanied by a large Stokes fluorescence shift, which was interpreted as the stabilized S₁(ππ*) excited state bonding to several water molecules via intermolecular hydrogen bonding.

Harmonizing accuracy and efficiency: A pragmatic approach to fragmentation of large molecules

Fragmentation methods offer an attractive alternative for ab initio treatment of large molecules and molecular clusters. However, balancing the accuracy and efficiency of these methods is a tight-rope-act. With this in view, we present an algorithm for automatic molecular fragmentation within Molecular Tailoring Approach (MTA) achieving this delicate balance. The automated code is tested out on a variety of molecules and clusters at the Hartree-Fock (HF)- and Møller-Plesset second order perturbation theory as well as density functional theory employing augmented Dunning basis sets. The results show remarkable accuracy and efficiency vis-à-vis the respective full calculations. Thus the present work forms an important step toward the development of an MTA-based black box code for implementation of HF as well as correlated quantum chemical calculations on large molecular systems.

Intramolecular O-H⋯O and C-H⋯O hydrogen bond cooperativity in D-glucopyranose and D-galactopyranose-A DFT/GIAO, QTAIM/IQA, and NCI approach

Density functional theory calculations are used to compute proton nuclear magnetic resonance (NMR) chemical shifts, interatomic distances, atom-atom interaction energies, and atomic charges for partial structures and conformers of α-D-glucopyranose, β-D-glucopyranose, and α-D-galactopyranose built up by introducing OH groups into 2-methyltetrahydropyran stepwisely. For the counterclockwise conformers, the most marked effects on the NMR shift and the charge on the OH1 proton are produced by OH2, those of OH3 and OH4 being somewhat smaller. This argues for a diminishing cooperative effect. The effect of OH6 depends on the configuration of the hydroxymethyl group and the position, axial or equatorial, of OH4, which controls hydrogen bonding in the 1,3-diol motif. Variations in the interaction energies reveal that a "new" hydrogen bond is sometimes formed at the expense of a preexisting one, probably due to geometrical constraints. Whereas previous work showed that complexing a conformer with pyridine affects only the nearest neighbour, successive OH groups increase the interaction energy of the N⋯H1 hydrogen bond and reduce its length. Analogous results are obtained for the clockwise conformers. The interaction energies for C-H⋯OH hydrogen bonding between axial CH protons and OH groups in certain conformers are much smaller than for O-H⋯OH bonds but they are largely covalent, whereas those of the latter are predominantly coulombic. These interactions are modified by complexation with pyridine in the same way as O-H⋯OH interactions: the computed NMR shifts of the CH protons increase, the atom-atom distances are shorter, and interaction energies are enhanced.

Anhydrous octyl-glucoside phase transition from lamellar to isotropic induced by electric and magnetic fields

A static deuterium nuclear magnetic resonance (²HNMR) technique (magnetic field, B = 7.05 T) was employed to monitor the thermotropic lamellar phase of the anhydrous 1:1 mixture sample of octyl-b-D-glucoside (βOG) and that of partially deuterium labelled at the alpha position on the chain, i.e.,βOG-d₂ In the absence of an electric field, the ²H NMR spectrum of the mixture gives a typical quadrupolar doublet representing the aligned lamellar phase. Upon heating to beyond the clearing temperature at 112 °C, this splitting converts to a single line expected for an isotropic phase. Simultaneous application of magnetic and electric fields (E = 0.4 MV/m) at 85 °C in the lamellar phase, whose direction was set to be parallel or perpendicular to the magnetic field, resulted in the change of the doublet into a single line and this recovers to the initial doublet with time for both experimental geometries. This implies E- and B-field-induced phase transitions from the lamellar to an isotropic phase and a recovery to the lamellar phase again with time. Moreover, these phase transformations are accompanied by a transient current. A similar observation was made in a computational study when an electric field was applied to a water cluster system. Increasing the field strength distorts the water cluster and weakens its hydrogen bonds leading to a structural breakdown beyond a threshold field-strength. Therefore, we suggest the observed field-induced transition is likely due to a structure change of the βOG lamellar assembly caused by the field effect and not due to Joule heating.

Accurate prediction of the structure and vibrational spectra of ionic liquid clusters with the generalized energy-based fragmentation approach: critical role of ion-pair-based fragmentation

The GEBF method with the ion-pair-based fragmentation has been developed to facilitate ab initio calculations of general ionic liquid clusters.

Microsolvation of F−

The effective charge in F⁻ is strong enough to induce water dissociation even for small molecularities in the gas phase.

Development of a Flexible-Monomer Two-Body Carbon Dioxide Potential and Its Application to Clusters up to (CO2 )13

A flexible-monomer two-body potential energy function was developed that approaches the high level CCSD(T)/CBS potential energy surface (PES) of carbon dioxide (CO₂ ) systems. This function was generated by fitting the electronic energies of unique CO₂ monomers and dimers to permutationally invariant polynomials. More than 200,000 CO₂ configurations were used to train the potential function. Comparisons of the PESs of six orientations of flexible CO₂ dimers were evaluated to demonstrate the accuracy of the potential. Furthermore, the potential function was used to determine the minimum energy structures of CO₂ clusters containing as many as 13 molecules. For isomers of (CO₂ )₃ , the potential demonstrated energetic agreement with the M06-2X functional and structural agreement of the B2PLYP-D functional at substantially reduced computational costs. A separate function, fit to MP2/aug-cc-pVDZ reference energies, was developed to directly compare the two-body potential to the ab initio MP2 level of theory. © 2017 Wiley Periodicals, Inc.

Hydrogen-Bonding Interactions in Luminescent Quinoline-Triazoles with Dominant 1D Crystals

Quinoline-triazoles 2-((4-(diethoxymethyl)-1H-1,2,3-triazol-1-yl)methyl)quinoline (1), 2-((4-(m-tolyl)-1H-1,2,3-triazol-1-yl)methyl)quinoline (2) and 2-((4-(p-tolyl)-1H-1,2,3-triazol-1-yl)methyl)quinoline (3) have been prepared with CuAAC click reactions and used as a model series to probe the relationship between lattice H-bonding interaction and crystal direction of growth. Crystals of 1-3 are 1D tape and prism shapes that correlate with their intermolecular and solvent 1D lattice H-bonding interactions. All compounds were thermally stable up to about 200 C and blue-green emissive in solution.