Nanostructuration of CoSi by mechanical milling and mechanical alloying

CoSi is an inexpensive thermoelectric material for medium temperature (200–500 °C). Its power factor is as large as the state of the art materials; however, its thermal conductivity is too large. Then, improving its thermoelectric performances implies increasing the scattering of phonons, which can be performed by nanostructuring the material. In this paper we investigate the effect of nanostructuration on the structure, microstructure, lattice dynamics and stability of CoSi. We obtained powders of about 13 nm by mechanical milling bulk CoSi for only four hours or by mechanical alloying pure elements for twelve hours. Nanostructuration induces a 0.1% expansion of the lattice parameter. Raman spectroscopy, associated to ab initio calculations, highlights the effectiveness of nanostructuration on phonon scattering, showing a reduction of the phonon relaxation time by as much as 80%. Powders are stable up to 450 °C; then grains coarsen and a partial degradation of the material occurs, probably due to silicon sublimation. Our results indicate that nanostructuration should be considered when interested to reduce CoSi thermal conductivity.     

Dynamic mechanical analysis

Using dynamic mechanical analysis (DMA) we have studied thermal degradation for a system containing a diglycidyl ether of bisphenol A (DGEBA) and 1,3-bisaminomethylcylohexane (1,3-BAC). The changes of dynamic mechanical properties during thermal degradation indicated a shift of the glass transition temperature (T _g) to higher temperatures and a decrease in the peak value of the dynamic loss factor (tan δ) with an increasing of aging time. The value of dynamic storage modulus (E′) at the rubbery state showed an increase with aging time, whiteE′ at the glassy state only underwent a moderate change with increased thermal degradation. From these results it can be argued that thermal degradation during the stage prior to the onset of the severe degradation involves structural changes in the epoxy system, as further crosslinking and loss of dangling chains in the crosslinked network.     

Specific quantum mechanical/molecular mechanical capping-potentials for biomolecular functional groups

We present a series of capping-potentials designed as link atoms to saturate dangling bonds at the quantum/classical interface within density functional theory-based hybrid QM/MM calculations. We aim at imitating the properties of different carbon-carbon bonds by means of monovalent analytic pseudopotentials. These effective potentials are optimized such that the perturbations of the quantum electronic density are minimized. This optimization is based on a stochastic scheme, which helps to avoid local minima trapping. For a series of common biomolecular groups, we find capping-potentials that outperform the more common hydrogen-capping in view of structural and spectroscopic properties. To demonstrate the transferability to complex systems, we also benchmark our potentials with a hydrogen-bonded dimer, yielding systematic improvements in structural and spectroscopic parameters.     

Quantum mechanical/molecular mechanical study of anthrax lethal factor catalysis

We report a hybrid quantum mechanical and molecular mechanical study of the catalysis of anthrax lethal factor. The calculations suggest that the zinc peptidase uses the same general base-general acid mechanism as in thermolysin and carboxypeptidase A, in which a zinc-bound water is activated by Glu687 to nucleophilically attack the scissile carbonyl carbon in the substrate. The catalysis is aided by an oxyanion hole formed by the zinc ion and the side chain of Tyr728, which provide stabilization for the fractionally charged carbonyl oxygen. The assigned role of Tyr728 differs from previous suggestions but is consistent with the established mechanism of other zinc proteases.     

Improved pseudobonds for combined ab initio quantum mechanical/molecular mechanical methods

The pseudobond approach offers a smooth connection at the quantum mechanical/molecular mechanical interface which passes through covalent bonds. It replaces the boundary atom of the environment part with a seven-valence-electron atom to form a pseudobond with the boundary atom of the active part [Y. Zhang, T. S. Lee, and W. Yang, J. Chem. Phys. 110, 46 (1999)]. In its original formulation, the seven-valence-electron boundary atom has the basis set of fluorine and a parametrized effective core potential. Up to now, only the Cps(sp3)-C(sp3) pseudobond has been successfully developed; thus in the case of proteins, it can only be used to cut the protein side chains. Here we employ a different formulation to construct this seven-valence-electron boundary atom, which has its own basis set as well as the effective core potential. We have not only further improved Cps(sp3)-C(sp3) pseudobond, but also developed Cps(sp3)-C(sp2,carbonyl) and Cps(sp3)-N(sp3) pseudobonds for the cutting of protein backbones and nucleic acid bases. The basis set and effective core potential for the seven-valence-electron boundary atom are independent of the molecular mechanical force field. Although the parametrization is performed with density functional calculations using hybrid B3LYP exchange-correlation functional, it is found that the same set of parameters is also applicable to Hartree-Fock and MP2 methods, as well as DFT calculations with other exchange-correlation functionals. Tests on a series of molecules yield very good structural, electronic, and energetic results in comparison with the corresponding full ab initio quantum mechanical calculations.     

Evaluation of an ab initio quantum mechanical/molecular mechanical hybrid-potential link-atom method

 Hybrid potentials have become a common tool in the study of many condensed-phase processes and are the subject of much active research. An important aspect of the formulation of a hybrid potential concerns how to handle covalent bonds between atoms that are described with different potentials and, most notably, those at the interface of the quantum mechanical (QM) and molecular mechanical (MM) regions. Several methods have been proposed to deal with this problem, ranging from the simple link-atom method to more sophisticated hybrid-orbital techniques. Although it has been heavily criticized, the link-atom method has probably been the most widely used in applications, especially with hybrid potentials that use semiempirical QM methods. Our aim in this paper has been to evaluate the link-atom method for ab initio QM/MM hybrid potentials and to compare the results it gives with those of previously published studies. Given its simplicity and robustness, we find that the link-atom method can produce results of comparable accuracy to other methods as long as the charge distribution on the MM atoms at the interface is treated appropriately. Received: 27 September 2002 / Accepted: 21 October 2002 / Published online: 8 January 2003 Correspondence to: M. J. Field e-mail: mjfield@ibs.fr Acknowledgements. The authors thank the Institut de Biologie Structurale – Jean-Pierre Ebel, the Commissariat à l'Energie Atomique and the Centre National de la Recherche Scientifique for support of this work.     

YinYang atom: a simple combined ab initio quantum mechanical molecular mechanical model

A simple interface is proposed for combined quantum mechanical (QM) molecular mechanical (MM) calculations for the systems where the QM and MM regions are connected through covalent bonds. Within this model, the atom that connects the two regions, called YinYang atom here, serves as an ordinary MM atom to other MM atoms and as a hydrogen-like atom to other QM atoms. Only one new empirical parameter is introduced to adjust the length of the connecting bond and is calibrated with the molecule propanol. This model is tested with the computation of equilibrium geometries and protonation energies for dozens of molecules. Special attention is paid on the influence of MM point charges on optimized geometry and protonation energy, and it is found that it is important to maintain local charge-neutrality in the MM region in order for the accurate calculation of the protonation and deprotonation energies. Overall the simple YinYang atom model yields comparable results to some other QM/MM models.     

Accelerating quantum mechanical/molecular mechanical sampling using pure molecular mechanical potential as an importance function: the case of effective fragment potential

Acceleration of sampling from a quantum mechanical/effective fragment mechanical (QM/EFP) potential is explored with effective fragment potential (EFP) as an importance function. EFP, generated on the fly, is found to be an excellent choice for an importance function for a QM/EFP potential. This technique is used to find nine stationary points of a blocked amino acid with twelve waters in a semi-automated way.     

Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP(+), while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme-reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP(+) cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP(+) cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.     

Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis

The molecular mechanism of hairpin ribozyme catalysis is studied with molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential with a recently developed semiempirical AM1/d-PhoT model for phosphoryl transfer reactions. Simulations are used to derive one- and two-dimensional potentials of mean force to examine specific reaction paths and assess the feasibility of proposed general acid and base mechanisms. Density-functional calculations of truncated active site models provide complementary insight to the simulation results. Key factors utilized by the hairpin ribozyme to enhance the rate of transphosphorylation are presented, and the roles of A38 and G8 as general acid and base catalysts are discussed. The computational results are consistent with available experimental data, provide support for a general acid/base mechanism played by functional groups on the nucleobases, and offer important insight into the ability of RNA to act as a catalyst without explicit participation by divalent metal ions.     

The catalytic activity of proline racemase: a quantum mechanical/molecular mechanical study

The enzyme proline racemase from the eukaryotic parasite Trypanosoma cruzi (responsible for endemic Chagas disease) catalyzes the reversible stereoinversion of chiral Calpha in proline. We employed a new combined quantum mechanical and molecular mechanical (QM/MM) potential to study the reaction mechanism of the enzyme. Three critical points were found: two almost isoenergetic minima (M1a and M2a), in which the enzyme is bound to L- and D-Pro, respectively, and a transition state (TSCa), unveiling a highly asynchronous concerted process. A systematic analysis was performed on the optimized geometries to point out the key role played by some residues in stabilizing the transition state.     

The dynamo library for molecular simulations using hybrid quantum mechanical and molecular mechanical potentials

The Dynamo module library has been developed for the simulation of molecular systems using hybrid quantum mechanical (QM) and molecular mechanical (MM) potentials. Dynamo is not a program package but is a library of Fortran 90 modules that can be employed by those interested in writing their own programs for performing molecular simulations. The library supports a range of different types of molecular calculation including geometry optimizations, reaction‐path determinations and molecular dynamics and Monte Carlo simulations. This article outlines the general structure and capabilities of the library and describes in detail Dynamo's semiempirical QM/MM hybrid potential. Results are presented to indicate three particular aspects of this implementation—the handling of long‐range nonbonding interactions, the nature of the boundary between the quantum mechanical and molecular mechanical atoms and how to perform path‐integral hybrid‐potential molecular dynamics simulations. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1088–1100, 2000     

The search for stationary points on a quantum mechanical/molecular mechanical potential-energy surface

 We propose a methodology to locate stationary points on a quantum mechanical/molecular mechanical potential-energy surface. This algorithm is based on a suitable approximation of an initial full Hessian matrix, either a modified Broyden–Fletcher–Goldfarg–Shanno or a Powell update formula for the location of, respectively, a minimum or a transition state, and the so-called rational function optimization. The latter avoids the Hessian matrix inversion required by a quasi-Newton–Raphson method. Some examples are presented and analyzed. Received: 16 July 2001 / Accepted: 9 October 2001 / Published online: 9 January 2002     

Study of the mechanical damping behavior of SBR-modified cement pastes by dynamic mechanical analyzer

Improvement in mechanical damping of SBR-modified cement pastes had been evaluated by dynamic mechanical analyzer. Specimens were fabricated and tested under 3-point flexure clamp at frequencies of 0.5–50 Hz or temperatures of ?30 to 70 °C. Significant improvement in damping was observed in cementitious-SBR composite specimens when SBR latex to cement ratio was 0.12, which is hypothesized to occur due to improvements in viscosity of cement paste. Furthermore, the SBR-modified cement pastes showed a decreased damping variation tendency with an increase of frequency. They also showed a peak damping variation tendency under the effect of the glass transition temperature. Based on the three element model, mechanical parameters are calculated by fitting the dynamic modulus of SBR-modified cement pastes.     

Generalized quantum mechanical two-centre problems

For the one-electron Schrödinger equation among the solutions of which the Slater-Zener-type functions can be found, it is shown, that it can be generalized to the two-centre case only in one way, if one demands separability in prolate spheroidal coordinates, and if in addition to the Coulomb term of the potential energy there shall be an additional function of the product r ₁ · r ₂ only. The generalized problem with a potential energy of the form V(r) = ? Z₁/r₁ ? Z₂/r₂ ? Q(R)/r₁r₂ is studied for the case of two equal centres Z ₁=Z₂=Z≧0 with regard to the existence and number of bound states. The results are extended as far as possible also to the case with unequal centres. For some examples with equal centres wave functions and correlation diagrams have been computed exactly for the lowest electronic states.     

A multifunctional mechanical Fourier spectrometer

A multifunctional instrument for measuring viscoelastic characteristics of polymeric materials is described. The operational principle of this instrument is based on the method of mechanical Fourier spectrometry. An incoming signal represents white noise that, like an outgoing signal, is expanded into a Fourier series in frequency. The use of a highly sensitive displacement transducer and a high-precision electromagnetic actuator makes it possible to collect spectral functions in a broad frequency range with a fair level of reliability. For the tests of flowable media, deformation is performed so that the test sample is compressed between two plane-parallel plates. In the tests of rubbery and solid samples, loading is accomplished via a three-point bending procedure or deformation of a membrane. The instrument makes it possible to perform tests on small samples (a volume less than 0.25 μl) and solid films with a thickness of 10–300 μm. The interval of complex elastic-modulus measurements ranges from 10 Pa to 50 MPa in the regime of a rheometer and from 10 kPa to 100 GPa in the regime of a solid-body spectrometer. The frequency interval is 0–1 Hz, and the temperature interval ranges from −100 to +300°C under the isothermic regime or under scanning at a heating rate from 0.1 to 10 K/min. Examples of the experimental results are presented.