Design of electroceramic materials using artificial neural networks and multiobjective evolutionary algorithms

We describe the computational design of electroceramic materials with optimal permittivity for application as electronic components. Given the difficulty of large-scale manufacture and characterization of these materials, including the theoretical prediction of their materials properties by conventional means, our approach is based on a recently established database containing composition and property information for a wide range of ceramic compounds. The electroceramic materials composition-function relationship is encapsulated by an artificial neural network which is used as one of the objectives in a multiobjective evolutionary algorithm. Evolutionary algorithms are stochastic optimization techniques which we employ to search for optimal materials based on chemical composition. The other objectives optimized include the reliability of the neural network prediction and the overall electrostatic charge of the material. The evolutionary algorithm searches for materials which simultaneously have high relative permittivity, minimum overall charge, and good prediction reliability. We find that we are able to predict a range of new electroceramic materials with varying degrees of reliability. In some cases the materials are similar to those contained in the database; in others, completely new materials are predicted.     

Transfer measurements for the Ti+Ni systems at near barrier energies

Large enhancements have been observed in the sub-barrier fusion cross sections for Ti+Ni systems in our previous studies. Coupled channel calculations incorporating couplings to 2⁺ and 3⁻ states failed to explain these enhancements completely. A possibilty of transfer channels contributing to the residual enhancements had been suggested. In order to investigate the role of relevant transfer channels, measurements of one- and two-nucleon transfer were carried out for ^46,48Ti+⁶¹Ni systems. The present paper gives the results of these studies.     

Energy controlled insertion of polar molecules in dense fluids

We present a method to search low energy configurations of polar molecules in the complex potential energy surfaces associated with dense fluids. The search is done in the configurational space of the translational and rotational degrees of freedom of the molecule, combining steepest-descent and Newton-Raphson steps which embed information on the average sizes of the potential energy wells obtained from prior inspection of the liquid structure. We perform a molecular dynamics simulation of a liquid water shell which demonstrates that the method enables fast and energy-controlled water molecule insertion in aqueous environments. The algorithm finds low energy configurations of incoming water molecules around three orders of magnitude faster than direct random insertion. This method represents an important step towards dynamic simulations of open systems and it may also prove useful for energy-biased ensemble average calculations of the chemical potential.     

Rapid and accurate prediction of binding free energies for saquinavir-bound HIV-1 proteases

To explain drug resistance by computer simulations at the molecular level, we first have to assess the accuracy of theoretical predictions. Herein we report an application of the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) technique to the ranking of binding affinities of the inhibitor saquinavir with the wild type (WT) and three resistant mutants of HIV-1 protease: L90M, G48V, and G48V/L90M. For each ligand-protein complex we report 10 ns of fully unrestrained molecular dynamics (MD) simulations with explicit solvent. We investigate convergence, internal consistency, and model dependency of MM/PBSA ligand binding energies. Converged enthalpy and entropy estimates produce ligand binding affinities within 1.5 kcal/mol of experimental values, with a remarkable level of correlation to the experimentally observed ranking of resistance levels. A detailed analysis of the enthalpic/entropic balance of drug-protease interactions explains resistance in L90M in terms of a higher vibrational entropy than in the WT complex, while G48V disrupts critical hydrogen bonds at the inhibitor's binding site and produces an altered, more unfavorable balance of Coulomb and polar desolvation energies.     

Force field validation for nucleic acid simulations: comparing energies and dynamics of a DNA dodecamer

Important questions exist regarding the quality of force fields used in molecular dynamics (MD) simulations and their interoperable use with other available MD implementations. NAMD is one of the most efficient and scalable parallel molecular dynamics codes for large-scale biomolecular simulations in the open source domain. It is the aim of this article to analyze and compare the dynamics of a benchmark DNA dodecamer d(CTTTTGCAAAAG)2 system, including its binding to a specific drug molecule arising from the use of various simulation protocols in NAMD using Amber98, with the dynamics arising from simulations of the same dodecamer using Amber98 in the AMBER package, one of the most well-established simulation codes for nucleic acids. Based upon a set of validation benchmarks, the details of which are discussed, we find that nucleic acid simulations using NAMD give meaningful results and that the essential features of the resulting dynamics are similar to those arising from the AMBER package. This sets the stage for reliable large-scale simulations of nucleic acids using NAMD.     

Theory, modelling and simulation in origins of life studies

Origins of life studies represent an exciting and highly multidisciplinary research field. In this review we focus on the contributions made by theory, modelling and simulation to addressing fundamental issues in the domain and the advances these approaches have helped to make in the field. Theoretical approaches will continue to make a major impact at the "systems chemistry" level based on the analysis of the remarkable properties of nonlinear catalytic chemical reaction networks, which arise due to the auto-catalytic and cross-catalytic nature of so many of the putative processes associated with self-replication and self-reproduction. In this way, we describe inter alia nonlinear kinetic models of RNA replication within a primordial Darwinian soup, the origins of homochirality and homochiral polymerization. We then discuss state-of-the-art computationally-based molecular modelling techniques that are currently being deployed to investigate various scenarios relevant to the origins of life.     

Role of host layer flexibility in DNA guest intercalation revealed by computer simulation of layered nanomaterials

Layered double hydroxides (LDHs) have been shown to form staged intermediate structures in experimental studies of intercalation. However, the mechanism by which staged structures are produced remains undetermined. Using molecular dynamics simulations, we show that LDHs are flexible enough to deform around bulky intercalants such as deoxyribonucleic acid (DNA). The flexibility of layered materials has previously been shown to affect the pathway by which staging occurs. We explore three possible intermediate structures which may form during intercalation of DNA into Mg2Al LDHs and study how the models differ energetically. When DNA strands are stacked directly on top of each other, the LDH system has a higher potential energy than when they are stacked in a staggered or interstratified structure. It is generally thought that staged intercalation occurs through a Daumas-Herold or a Rudorff model. We find, on average, greater diffusion coefficients for DNA strands in a Daumas-Herold configuration compared to a Rudorff model and a stage-1 structure. Our simulations provide evidence for the presence of peristaltic modes of motion within Daumas-Herold configurations. This is confirmed by spectral analysis of the thickness variation of the basal spacing. Peristaltic modes are more prominent in the Daumas-Herold structure compared to the Rudorff and stage-1 structures and support a mechanism by means of which bulky intercalated molecules such as DNA rapidly diffuse within an LDH interlayer.