Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “quantum computing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,603 records · Page 89Linked to original sources

Theoretical study of anion binding to calix[4]pyrrole: the effects of solvent, fluorine substitution, cosolute, and water traces.

The binding of different anions to calix[4]pyrrole has been studied by means of molecular dynamics coupled to thermodynamic integration calculations. The effect of different apolar solvents, octafluoro substitution, and the change in binding free energy derived from the presence of cosolute and water traces (the hydrated salt used to introduce the anion in the solution) were examined. Calculations allow us to rationalize the differential binding of ions to calix[4]pyrrole and octafluorocalix[4]pyrrole as well as to predict the behavior in new solvents for which experimental data are not available yet. It is found that both calix[4]pyrrole and octafluorocalix[4]pyrrole have a dramatic preference for F- in the gas phase and pure aprotic solvents, but the situation can change dramatically in protic solvents or in the presence of the hydrated cation which is used as cosolute of the anion. Overall, our results provide interesting clues for a better understanding of the process detected experimentally as "binding".

Anions↗

Dynamics of an enzymatic substitution reaction in haloalkane dehalogenase.

Reactive flux molecular dynamics simulations have been carried out using a combined QM/MM potential to study the dynamics of the nucleophilic substitution reaction of dichloroethane by a carboxylate group in haloalkane dehalogenase and in water. We found that protein dynamics accelerates the reaction rate by a factor of 2 over the uncatalyzed reaction. Compared to the thermodynamic effect in barrier reduction, protein dynamic contribution is relatively small. However, analyses of the friction kernel reveal that the origins of the reaction dynamics in water and in the enzyme are different. In aqueous solution, there is significant electrostatic solvation effect, which is reflected by the slow reorganization relaxation of the solvent. On the other hand, there is no strong electrostatic coupling in the enzyme and the major effect on reaction coordinate motion is intramolecular energy relaxation.

Computer Simulation↗

Understanding the template preorganization step of an artificial arginine receptor.

A biomimetic complex which mimics the arginine-phosphonate diester interaction of the arginine fork is investigated with respect to structure and energetics of stable configurations. Within this work, we provide knowledge on local minima of the isolated system obtained from first-principles calculations. Non-negligible solvation effects are studied in a microsolvation approach. The interactions which govern the structural patterns of molecular recognition in this tweezer-guest complex can be significantly modulated by the action of hydrogen bond accepting and donating solvent molecules, such as dimethyl sulfoxide or water, which were present in experimental investigations on this system. Different tweezer-guest structures are evaluated with respect to their temperature-dependent thermodynamical properties as products of the first association reaction step of the bisphosphonate tweezer template and the guanidinium moiety.

Arginine↗

Cope elimination: elucidation of solvent effects from QM/MM simulations.

The Cope elimination reactions for threo- and erythro-N,N-dimethyl-3-phenyl-2-butylamine oxide have been investigated using QM/MM calculations in water, THF, and DMSO. The aprotic solvents provide up to million-fold rate accelerations. The effects of solvation on the reactants, transition structures, and rates of reaction are elucidated here using two-dimensional potentials of mean force (PMF) derived from free-energy perturbation calculations in Monte Carlo simulations (MC/FEP). The resultant free energies of activation in solution are in close agreement with experiment. Ab initio calculations at the MP2/6-311+G-(2d,p) level using the PCM continuum solvent model were also carried out; however, only the QM/MM methodology was able to reproduce the large rate increases in proceeding from water to the dipolar aprotic solvents. Solute-solvent interaction energies and radial distribution functions are also analyzed and show that poorer solvation of the reactant in the aprotic solvents is primarily responsible for the observed rate enhancements. It is found that the amine oxide oxygen is the acceptor of three hydrogen bonds from water molecules for the reactant but only one to two weaker ones at the transition state. The overall quantitative success of the computations supports the present QM/MM/MC approach, featuring PDDG/PM3 as the QM method.

Alkenes↗

Classification of drugs by discriminant analysis using fragment molecular connectivity values.

An investigation was made into the use of linear and quadratic discriminant analysis, along with K nearest-neighbor analysis, in the classification of a set of 51 compounds which were divided into five therapeutic categories. By superimposing each compound on a pattern structure, as first proposed by Cammarata, eight positions were assigned on the molecule. Each position was coded with the numerical value of a descriptor index. Relative molar refraction, which was the index used by Cammarata, was compared with a number of molecular connective indices. For each of the indices studied, it was found that only four of the eight positions contributed significantly to between-class differences. It was also found that first-order molecular connectivity, calculated as the sum of the contributions of each of the bonds joining a given position, resulted in consistently fewer misclassifications as compared with the other indices. Using first-order molecular connectivity, validation procedures were performed on the original set of compounds, on random samples drawn from this set, and on a set of ten compounds not included in the analysis. The results obtained were highly data dependent, but they, nevertheless, suggest that molecular connectivity indices should prove useful in structural classification procedures.

Analysis of Variance↗

First-principles calculation of 17O and 25Mg NMR shieldings in MgO at finite temperature: rovibrational effect in solids.

The temperature dependence of (17)O and (25)Mg NMR chemical shifts in solid MgO have been calculated using a first-principles approach. Density functional theory, pseudopotentials, a plane-wave basis set, and periodic boundary conditions were used both to describe the motion of the nuclei and to compute the NMR chemical shifts. The chemical shifts were obtained using the gauge-including projector augmented wave method. In a crystalline solid, the temperature dependence is due to both (i) the variation of the averaged equilibrium structure and (ii) the fluctuation of the atoms around this structure. In MgO, the equilibrium structure at each temperature is uniquely defined by the cubic lattice parameters, which we take from experiment. We evaluate the effect of the fluctuations within a quasiharmonic approximation. In particular, the dynamical matrix, defining the harmonic Hamiltonian, has been computed for each equilibrium volume. This harmonic Hamiltonian was used to generate nuclear configurations that obey quantum statistical mechanics. The chemical shifts were averaged over these nuclear configurations. The results reproduce the previously published experimental NMR data measured on MgO between room temperature and 1000 degrees C. It is shown that the chemical shift behavior with temperature cannot be explained by thermal expansion alone. Vibrational corrections due to the fluctuations of atoms around their equilibrium position are crucial to reproduce the experimental results.

Journal Article↗

Why isn't the N20 dodecahedron ideal for three-coordinate nitrogen?

Nitrogen molecules are the focus of much attention for their potential as high-energy density materials. The usefulness of such molecules as energy sources depends on the stability of the molecules with respect to dissociation. Many such molecules dissociate too easily to be a stable fuel, and the reasons for such instability are related to the details of structure and bonding of the molecule. Such details will be examined for isomers of the molecule N20. N20 has a highly symmetric isomer in which the 20 atoms occupy the vertexes of a dodecahedron. This isomer is a cage molecule in which all of the faces are regular pentagons that have interior angles of 108 degrees . These angles are very close to the known bond angles in very stable nitrogen compounds such as ammonia. Such a structure with only pentagons should provide an ideal bonding environment for three-coordinate nitrogen. However, by use of theoretical calculations including density-functional theory and fourth-order perturbation theory, along with the Dunning correlation-consistent basis sets, it will be shown that dodecahedral N20 is not the most stable cage for N20. The reasons why will be discussed in terms of the structure and bonding of the molecules.

Computer Simulation↗

Computational model of hole transport in DNA.

A computational model based on the molecular dynamics (MD) simulation for the hole transport in DNA has been developed and applied to study hole current in DNA strands consisting of different numbers of GC pairs. The approach is based on the hopping mechanism which is thermally activated. The calculations show that the hole hopping intensifies with the temperature and the transport rate increases in agreement with the experimental evidence. It is also determined that the degree of structural ordering in the DNA strand enhances the hole conductivity and reasons are provided why this may occur.

Base Pairing↗

Electric field fluctuations drive vibrational dephasing in water.

We present a microscopic description of the vibrational spectroscopy of the OH stretch of HOD in liquid D(2)O. Our model predicts that OH frequency correlations decay with a sharp and rapid ( approximately 35 fs) decrease, followed by a beat at approximately 125 fs from intermolecular oxygen vibrations. On a short time scale ( approximately 200 fs), ultrafast infrared spectroscopy of the OH stretch is sensitive to localized intermolecular motions. For times longer than approximately 200 fs, cooperative molecular rearrangements drive dephasing. The interplay of electric field fluctuations, both local and cooperative, dictate vibrational frequency shifts and destroy vibrational coherence in water.

Computer Simulation↗

Vibrational spectroscopy of the G...C base pair: experiment, harmonic and anharmonic calculations, and the nature of the anharmonic couplings.

The results of harmonic and anharmonic frequency calculations on a guanine-cytosine complex with an enolic structure (a tautomeric form with cytosine in the enol form and with a hydrogen at the 7-position on guanine) are presented and compared to gas-phase IR-UV double resonance spectral data. Harmonic frequencies were obtained at the RI-MP2/cc-pVDZ, RI-MP2/TZVPP, and semiempirical PM3 levels of electronic structure theory. Anharmonic frequencies were obtained by the CC-VSCF method with improved PM3 potential surfaces; the improved PM3 potential surfaces are obtained from standard PM3 theory by coordinate scaling such that the improved PM3 harmonic frequencies are the same as those computed at the RI-MP2/cc-pVDZ level. Comparison of the data with experimental results indicates that the average absolute percentage deviation for the methods is 2.6% for harmonic RI-MP2/cc-pVDZ (3.0% with the inclusion of a 0.956 scaling factor that compensates for anharmonicity), 2.5% for harmonic RI-MP2/TZVPP (2.9% with a 0.956 anharmonicity factor included), and 2.3% for adapted PM3 CC-VSCF; the empirical scaling factor for the ab initio harmonic calculations improves the stretching frequencies but decreases the accuracy of the other mode frequencies. The agreement with experiment supports the adequacy of the improved PM3 potentials for describing the anharmonic force field of the G...C base pair in the spectroscopically probed region. These results may be useful for the prediction of the pathways of vibrational energy flow upon excitation of this system. The anharmonic calculations indicate that anharmonicity along single mode coordinates can be significant for simple stretching modes. For several other cases, coupling between different vibrational modes provides the main contribution to anharmonicity. Examples of strongly anharmonically coupled modes are the symmetric stretch and group torsion of the hydrogen-bonded NH2 group on guanine, the OH stretch and torsion of the enol group on cytosine, and the NH stretch and NH out-of-plane bend of the non-hydrogen-bonded NH group on guanine.

Base Pairing↗

Length and substituent-scrambling energies of parent and halogen-substituted conjugated polyynes.

Conjugated polyynes are a class of species of diverse and increasing interest. Length-scrambling and substituent scrambling reaction energies were examined using ab initio quantum chemistry calculations to investigate issues concerning the energetic effects of the molecular ends (substituent communication). Computations were performed for the parent, monohalogenated, and dihalogenated (F, Cl, Br, I) polyynes of up to 60 carbon atoms. A study of resonance effects using natural resonance theory and bond lengths demonstrates lone-pair-donating effects that increase in the series F < Cl < Br < I, but run counter to the halogen inductive effects which decrease in this series and dominate energetic effects.

Journal Article↗

An evaluation of various computational methods for the treatment of organoselenium compounds.

A reliable computational method for the prediction of organoselenium geometries and bond dissociation energies (BDEs) has been determined on the basis of the performance of density functional theory (DFT: B3LYP and B3PW91) and ab initio molecular orbital procedures (Hartree-Fock (HF)) in conjunction with various Pople basis sets including (but not limited to) the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p), 6-311G(2df,p), and 6-311G(3df,3pd) sets. Predicted geometries and BDEs are compared with available experimental data and quadratic configuration interaction including single and double substitutions (QCISD) results. The B3PW91/6-311G(2df,p) level of theory is recommended for the prediction of the geometries and energetics of organoselenium compounds.

Computer Simulation↗

Photodissociation of cyclobutyl bromide at 234 nm studied using velocity map imaging.

This study investigates the 234 nm photodissociation dynamics of cyclobutyl bromide using a two-dimensional photofragment velocity imaging technique. The spin-orbit ground- and excited-state Br(2P) atoms are state-selectively detected via [2+1] resonance enhanced multiphoton ionization (REMPI), whereas the cyclobutyl radicals are ionized using 157 nm laser light. The Br(2P(3/2)) and the Br(2P(1/2)) atoms and their c-C4H7 radical cofragments evidence a single-peaked, Gaussian-shaped translational energy distribution ranging from approximately 14 to approximately 39 kcal/mol and angular distributions with significant parallel character. The Br(2P(1/2))/ Br(2P(3/2)) spin-orbit branching ratio is determined to be 0.11 +/- 0.07 by momentum match between the Br(2P) photofragments and the recoiling c-C4H7 fragments, assuming a uniform photoionization probability of the c-C4H7 radicals with an internal energy range of 10-35 kcal/mol. The REMPI line strength ratio for the detection of Br(2P(3/2)) and Br(2P(1/2)) atoms at 233.681 and 234.021 nm, respectively, is therefore derived to be 0.10 +/- 0.07. The measured recoil kinetic energies of the c-C4H7 radicals, and the resulting distribution of internal energies, indicates some of the radicals are formed with total internal energies above the barrier to isomerization and subsequent dissociation, but our analysis indicates they may be stable due to the substantial fraction of the internal energy which is partitioned to rotational energy of the radicals.

Bromides↗

Interpretation of hydrogen bonding in the weak and strong regions using conceptual DFT descriptors.

Hydrogen bonding is among the most fundamental interactions in biology and chemistry, providing an extra stabilization of 1-40 kcal/mol to the molecular systems involved. This wide range of stabilization energy underlines the need for a general and comprehensive theory that will explain the formation of hydrogen bonds. While a simple electrostatic model is adequate to describe the bonding patterns in the weak and moderate hydrogen bond regimes, strong hydrogen bonds, on the other hand, require a more complete theory due to the appearance of covalent interactions. In this study, conceptual DFT tools such as local hardness, eta(r) and local softness, s(r), have been used in order to get an alternative view on solving this hydrogen-bonding puzzle as described by Gilli et al. [J. Mol. Struct. 2000, 552, 1]. A series of both homonuclear and heteronuclear resonance-assisted hydrogen bonds of the types O-H...N, N-H...O, N-H...N, and O-H...O with strength varying from weak to very strong have been studied. First of all, DeltaPA and DeltapK(a) values were calculated and correlated to the hydrogen bond energy. Then the electrostatic effects were examined as hard-hard interactions accessible through molecular electrostatic potential, natural population analysis (NPA) charge, and local hardness calculations. Finally, secondary soft-soft interaction effects were entered into the picture described by the local softness values, providing insight into the covalent character of the strong hydrogen bonds.

Acrylates↗

High-resolution infrared spectroscopy of HCN-Znn (n = 1-4) clusters: structure determination and comparisons with theory.

High-resolution infrared laser spectroscopy has been used to obtain rotationally resolved spectra of HCN-Zn(n) (n = 1-4) complexes formed in helium nanodroplets. In the present study the droplets passed through a metal oven, where the zinc vapor pressure was adjusted until one or more atoms were captured by the droplets. A second pickup cell was then used to dope the droplets with a single HCN molecule. Rotationally resolved infrared spectra are obtained for all of these complexes, providing valuable information concerning their structures. Stark spectra are reported and used to determine the corresponding permanent electric dipole moments. Ab initio calculations are also reported for these complexes for comparison with the experimental results.

Computer Simulation↗

Orbital signatures of methyl in L-alanine.

Molecular orbital signatures of the methyl substituent in L-alanine have been identified with respect to those of glycine from information obtained in coordinate and momentum space, using dual space analysis. Electronic structural information in coordinate space is obtained using ab initio (MP2/TZVP) and density functional theory (B3LYP/TZVP) methods, from which the Dyson orbitals are simulated based on the plane wave impulse approximation into momentum space. In comparison to glycine, relaxation in geometry and valence orbitals in L-alanine is found as a result of the attachment of the methyl group. Five orbitals rather than four orbitals are identified as methyl signatures. That is, orbital 6a in the core shell, orbitals 11a and 12a in the inner valence shell, and orbitals 19a and 20a in the outer valence shell. In the inner valence shell, the attachment of methyl to glycine causes a splitting of its orbital 10a' into orbitals 11a and 12a of L-alanine, whereas in the outer valence shell the methyl group results in an insertion of an additional orbital pair of 19a and 20a. The frontier molecular orbitals, 24a and 23a, are found without any significant role in the methylation of glycine.

Alanine↗

Geometry and solvent dependence of the electronic spectra of the amide group and consequences for peptide circular dichroism.

The influence of geometry variations and solvent environment of N-methylacetamide on its energies and absorption intensities was systematically analyzed with the aid of the time-dependent density functional theory (TD DFT). Selective and often complicated reactions of individual electronic levels on the perturbations were found important for the resultant spectral profile. For example, the n-pi band position varied by tens of nanometers due to the C=O bond length oscillations, while it was rather unaffected by surrounding water. On the contrary, pi-pi type transition energies and intensities were broadly dispersed by the aqueous environment but exhibited a modest coordinate dependence. A simple electrostatic model used previously for absorption in the IR region (J. Chem. Phys. 2005, 122, 144501) explained these changes only partially. Additionally, electronic transfer between the solute and the solvent had to be considered for faithful modeling of the ultraviolet light absorption. The inclusion of the environment and dynamics in the modeling then provided more accurate positions, intensities, and realistic inhomogeneous widths of spectral lines. These factors were found important for absorption and circular dichroism spectra of larger peptides and proteins. This was demonstrated with a combined DFT/coupled oscillator model providing principal features observed in electronic circular dichroism spectra of standard peptide conformations.

Acetamides↗

Computational study of pharmacophores: beta-lactams.

The thermochemistry of bicyclic lactams (penams, penems, cephams, cephems), which are key pharmacophores in beta-lactam antibiotics, has been investigated by high-level ab initio methods. Particular attention has been paid to estimating magnitudes of amide resonance (ARE) and ring strain (RSE) in the four-member lactam ring because these quantities are difficult to measure and distinguish experimentally. The ring strain destabilization effect is greater than the stabilization arising from amide resonance. However, in cephemes the amide resonance stabilization slightly exceeds destabilization due to the beta-lactam ring strain.

Computer Simulation↗