Direct computation of parameters for accurate polarizable force fields
Verstraelen, Toon Vandenbrande, Steven; Ayers, Paul W.
2014-11-21
We present an improved electronic linear response model to incorporate polarization and charge-transfer effects in polarizable force fields. This model is a generalization of the Atom-Condensed Kohn-Sham Density Functional Theory (DFT), approximated to second order (ACKS2): it can now be defined with any underlying variational theory (next to KS-DFT) and it can include atomic multipoles and off-center basis functions. Parameters in this model are computed efficiently as expectation values of an electronic wavefunction, obviating the need for their calibration, regularization, and manual tuning. In the limit of a complete density and potential basis set in the ACKS2 model, the linear response properties of the underlying theory for a given molecular geometry are reproduced exactly. A numerical validation with a test set of 110 molecules shows that very accurate models can already be obtained with fluctuating charges and dipoles. These features greatly facilitate the development of polarizable force fields.
Developing accurate molecular mechanics force fields for conjugated molecular systems.
Do, Hainam; Troisi, Alessandro
2015-10-14
A rapid method to parameterize the intramolecular component of classical force fields for complex conjugated molecules is proposed. The method is based on a procedure of force matching with a reference electronic structure calculation. It is particularly suitable for those applications where molecular dynamics simulations are used to generate structures that are therefore analysed by electronic structure methods, because it is possible to build force fields that are consistent with electronic structure calculations that follow classical simulations. Such applications are commonly encountered in organic electronics, spectroscopy of complex systems and photobiology (e.g. photosynthetic systems). We illustrate the method by parameterizing the force fields of a molecule used in molecular semiconductors (2,2-dicyanovinyl-capped S,N-heteropentacene or DCV-SN5), a polymeric semiconductor (thieno[3,2-b]thiophene-diketopyrrolopyrrole TT-DPP) and a chromophore embedded in a protein environment (15,16-dihydrobiliverdin or DBV) where several hundreds of parameters need to be optimized in parallel. PMID:26349916
Accurate force fields and methods for modelling organic molecular crystals at finite temperatures.
Nyman, Jonas; Pundyke, Orla Sheehan; Day, Graeme M
2016-06-21
We present an assessment of the performance of several force fields for modelling intermolecular interactions in organic molecular crystals using the X23 benchmark set. The performance of the force fields is compared to several popular dispersion corrected density functional methods. In addition, we present our implementation of lattice vibrational free energy calculations in the quasi-harmonic approximation, using several methods to account for phonon dispersion. This allows us to also benchmark the force fields' reproduction of finite temperature crystal structures. The results demonstrate that anisotropic atom-atom multipole-based force fields can be as accurate as several popular DFT-D methods, but have errors 2-3 times larger than the current best DFT-D methods. The largest error in the examined force fields is a systematic underestimation of the (absolute) lattice energy. PMID:27230942
Lee, M.W.; Meuwly, M.
2013-01-01
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories.
Lee, Myung Won; Meuwly, Markus
2013-12-14
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories. PMID:24170171
Geng, Hao; Jiang, Fan; Wu, Yun-Dong
2016-05-19
Cyclic peptides (CPs) are promising candidates for drugs, chemical biology tools, and self-assembling nanomaterials. However, the development of reliable and accurate computational methods for their structure prediction has been challenging. Here, 20 all-trans CPs of 5-12 residues selected from Cambridge Structure Database have been simulated using replica-exchange molecular dynamics with four different force fields. Our recently developed residue-specific force fields RSFF1 and RSFF2 can correctly identify the crystal-like conformations of more than half CPs as the most populated conformation. The RSFF2 performs the best, which consistently predicts the crystal structures of 17 out of 20 CPs with rmsd < 1.1 Å. We also compared the backbone (ϕ, ψ) sampling of residues in CPs with those in short linear peptides and in globular proteins. In general, unlike linear peptides, CPs have local conformational free energies and entropies quite similar to globular proteins. PMID:27128113
An Accurate Quartic Force Field and Vibrational Frequencies for HNO and DNO
NASA Technical Reports Server (NTRS)
Dateo, Christopher E.; Lee, Timothy J.; Schwenke, David W.
1994-01-01
An accurate ab initio quartic force field for HNO has been determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T), in conjunction with the correlation consistent polarized valence triple zeta (cc-pVTZ) basis set. Improved harmonic frequencies were determined with the cc-pVQZ basis set. Fundamental vibrational frequencies were determined using a second-order perturbation theory analysis and also using variational calculations. The N-0 stretch and bending fundamentals are determined well from both vibrational analyses. The H-N stretch, however, is shown to have an unusually large anharmonic correction, and is not well determined using second-order perturbation theory. The H-N fundamental is well determined from the variational calculations, demonstrating the quality of the ab initio quartic force field. The zero-point energy of HNO that should be used in isodesmic reactions is also discussed.
Accurate ab initio quartic force fields for the ions HCO(+) and HOC(+)
NASA Astrophysics Data System (ADS)
Martin, J. M. L.; Taylor, Peter R.; Lee, Timothy J.
1993-07-01
The quartic force fields of HCO(+) and HOC(+) have been computed using augmented coupled cluster methods and basis sets of spdf and spdfg quality. Calculations on HCN, CO, and N2 have been performed to assist in calibrating the computed results. Going from an spdf to an spdfg basis shortens triple bonds by about 0.004 A, and increases the corresponding harmonic frequency by 10-20/cm, leaving bond distances about 0.003 A too long and triple bond stretching frequencies about 5/cm too low. Accurate estimates for the bond distances, fundamental frequencies, and thermochemical quantities are given. HOC(+) lies 37.8 +/- 0.5 kcal/mol (0 K) above HCO(+); the classical barrier height for proton exchange is 76.7 +/- 1.0 kcal/mol.
Accurate ab initio Quartic Force Fields of Cyclic and Bent HC2N Isomers
NASA Technical Reports Server (NTRS)
Inostroza, Natalia; Huang, Xinchuan; Lee, Timothy J.
2012-01-01
Highly correlated ab initio quartic force field (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC2N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T). Dunning s correlation-consistent basis sets cc-pVXZ, X=3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schr dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by perturbation theory). On the other hand, this procedure (a QFF together with either perturbation theory or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC2N. All three isomers possess significant dipole moments, 3.05D, 3.06D, and 1.71D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.
Are current atomistic force fields accurate enough to study proteins in crowded environments?
Petrov, Drazen; Zagrovic, Bojan
2014-05-01
The high concentration of macromolecules in the crowded cellular interior influences different thermodynamic and kinetic properties of proteins, including their structural stabilities, intermolecular binding affinities and enzymatic rates. Moreover, various structural biology methods, such as NMR or different spectroscopies, typically involve samples with relatively high protein concentration. Due to large sampling requirements, however, the accuracy of classical molecular dynamics (MD) simulations in capturing protein behavior at high concentration still remains largely untested. Here, we use explicit-solvent MD simulations and a total of 6.4 µs of simulated time to study wild-type (folded) and oxidatively damaged (unfolded) forms of villin headpiece at 6 mM and 9.2 mM protein concentration. We first perform an exhaustive set of simulations with multiple protein molecules in the simulation box using GROMOS 45a3 and 54a7 force fields together with different types of electrostatics treatment and solution ionic strengths. Surprisingly, the two villin headpiece variants exhibit similar aggregation behavior, despite the fact that their estimated aggregation propensities markedly differ. Importantly, regardless of the simulation protocol applied, wild-type villin headpiece consistently aggregates even under conditions at which it is experimentally known to be soluble. We demonstrate that aggregation is accompanied by a large decrease in the total potential energy, with not only hydrophobic, but also polar residues and backbone contributing substantially. The same effect is directly observed for two other major atomistic force fields (AMBER99SB-ILDN and CHARMM22-CMAP) as well as indirectly shown for additional two (AMBER94, OPLS-AAL), and is possibly due to a general overestimation of the potential energy of protein-protein interactions at the expense of water-water and water-protein interactions. Overall, our results suggest that current MD force fields may distort the
Are Current Atomistic Force Fields Accurate Enough to Study Proteins in Crowded Environments?
Petrov, Drazen; Zagrovic, Bojan
2014-01-01
The high concentration of macromolecules in the crowded cellular interior influences different thermodynamic and kinetic properties of proteins, including their structural stabilities, intermolecular binding affinities and enzymatic rates. Moreover, various structural biology methods, such as NMR or different spectroscopies, typically involve samples with relatively high protein concentration. Due to large sampling requirements, however, the accuracy of classical molecular dynamics (MD) simulations in capturing protein behavior at high concentration still remains largely untested. Here, we use explicit-solvent MD simulations and a total of 6.4 µs of simulated time to study wild-type (folded) and oxidatively damaged (unfolded) forms of villin headpiece at 6 mM and 9.2 mM protein concentration. We first perform an exhaustive set of simulations with multiple protein molecules in the simulation box using GROMOS 45a3 and 54a7 force fields together with different types of electrostatics treatment and solution ionic strengths. Surprisingly, the two villin headpiece variants exhibit similar aggregation behavior, despite the fact that their estimated aggregation propensities markedly differ. Importantly, regardless of the simulation protocol applied, wild-type villin headpiece consistently aggregates even under conditions at which it is experimentally known to be soluble. We demonstrate that aggregation is accompanied by a large decrease in the total potential energy, with not only hydrophobic, but also polar residues and backbone contributing substantially. The same effect is directly observed for two other major atomistic force fields (AMBER99SB-ILDN and CHARMM22-CMAP) as well as indirectly shown for additional two (AMBER94, OPLS-AAL), and is possibly due to a general overestimation of the potential energy of protein-protein interactions at the expense of water-water and water-protein interactions. Overall, our results suggest that current MD force fields may distort the
Sprenger, K G; Jaeger, Vance W; Pfaendtner, Jim
2015-05-01
We have applied molecular dynamics to calculate thermodynamic and transport properties of a set of 19 room-temperature ionic liquids. Since accurately simulating the thermophysical properties of solvents strongly depends upon the force field of choice, we tested the accuracy of the general AMBER force field, without refinement, for the case of ionic liquids. Electrostatic point charges were developed using ab initio calculations and a charge scaling factor of 0.8 to more accurately predict dynamic properties. The density, heat capacity, molar enthalpy of vaporization, self-diffusivity, and shear viscosity of the ionic liquids were computed and compared to experimentally available data, and good agreement across a wide range of cation and anion types was observed. Results show that, for a wide range of ionic liquids, the general AMBER force field, with no tuning of parameters, can reproduce a variety of thermodynamic and transport properties with similar accuracy to that of other published, often IL-specific, force fields. PMID:25853313
NASA Astrophysics Data System (ADS)
Liu, Qianlong; Reifsnider, Kenneth
2012-11-01
The basis of dielectrophoresis (DEP) is the prediction of the force and torque on particles. The classical approach to the prediction is based on the effective moment method, which, however, is an approximate approach, assumes infinitesimal particles. Therefore, it is well-known that for finite-sized particles, the DEP approximation is inaccurate as the mutual field, particle, wall interactions become strong, a situation presently attracting extensive research for practical significant applications. In the present talk, we provide accurate calculations of the force and torque on the particles from first principles, by directly resolving the local geometry and properties and accurately accounting for the mutual interactions for finite-sized particles with both dielectric polarization and conduction in a sinusoidally steady-state electric field. Since the approach has a significant advantage, compared to other numerical methods, to efficiently simulate many closely packed particles, it provides an important, unique, and accurate technique to investigate complex DEP phenomena, for example heterogeneous mixtures containing particle chains, nanoparticle assembly, biological cells, non-spherical effects, etc. This study was supported by the Department of Energy under funding for an EFRC (the HeteroFoaM Center), grant no. DE-SC0001061.
Romanov, V N; Cygan, R T; Myshakin, E M
2012-06-21
Naturally occurring clay minerals provide a distinctive material for carbon capture and carbon dioxide sequestration. Swelling clay minerals, such as the smectite variety, possess an aluminosilicate structure that is controlled by low-charge layers that readily expand to accommodate water molecules and, potentially, CO2. Recent experimental studies have demonstrated the efficacy of intercalating CO2 in the interlayer of layered clays, but little is known about the molecular mechanisms of the process and the extent of carbon capture as a function of clay charge and structure. A series of molecular dynamics simulations and vibrational analyses have been completed to assess the molecular interactions associated with incorporation of CO2 and H2O in the interlayer of montmorillonite clay and to help validate the models with experimental observation. An accurate and fully flexible set of interatomic potentials for CO2 is developed and combined with Clayff potentials to help evaluate the intercalation mechanism and examine the effect of molecular flexibility onthe diffusion rate of CO2 in water.
A procedure for computing accurate ab initio quartic force fields: Application to HO2+ and H2O
NASA Astrophysics Data System (ADS)
Huang, Xinchuan; Lee, Timothy J.
2008-07-01
A procedure for the calculation of molecular quartic force fields (QFFs) is proposed and investigated. The goal is to generate highly accurate ab initio QFFs that include many of the so-called ``small'' effects that are necessary to achieve high accuracy. The small effects investigated in the present study include correlation of the core electrons (core correlation), extrapolation to the one-particle basis set limit, correction for scalar relativistic contributions, correction for higher-order correlation effects, and inclusion of diffuse functions in the one-particle basis set. The procedure is flexible enough to allow for some effects to be computed directly, while others may be added as corrections. A single grid of points is used and is centered about an initial reference geometry that is designed to be as close as possible to the final ab initio equilibrium structure (with all effects included). It is shown that the least-squares fit of the QFF is not compromised by the added corrections, and the balance between elimination of contamination from higher-order force constants while retaining energy differences large enough to yield meaningful quartic force constants is essentially unchanged from the standard procedures we have used for many years. The initial QFF determined from the least-squares fit is transformed to the exact minimum in order to eliminate gradient terms and allow for the use of second-order perturbation theory for evaluation of spectroscopic constants. It is shown that this step has essentially no effect on the quality of the QFF largely because the initial reference structure is, by design, very close to the final ab initio equilibrium structure. The procedure is used to compute an accurate, purely ab initio QFF for the H2O molecule, which is used as a benchmark test case. The procedure is then applied to the ground and first excited electronic states of the HO2+ molecular cation. Fundamental vibrational frequencies and spectroscopic
NASA Astrophysics Data System (ADS)
Sagui, Celeste; Pedersen, Lee G.; Darden, Thomas A.
2004-01-01
The accurate simulation of biologically active macromolecules faces serious limitations that originate in the treatment of electrostatics in the empirical force fields. The current use of "partial charges" is a significant source of errors, since these vary widely with different conformations. By contrast, the molecular electrostatic potential (MEP) obtained through the use of a distributed multipole moment description, has been shown to converge to the quantum MEP outside the van der Waals surface, when higher order multipoles are used. However, in spite of the considerable improvement to the representation of the electronic cloud, higher order multipoles are not part of current classical biomolecular force fields due to the excessive computational cost. In this paper we present an efficient formalism for the treatment of higher order multipoles in Cartesian tensor formalism. The Ewald "direct sum" is evaluated through a McMurchie-Davidson formalism [L. McMurchie and E. Davidson, J. Comput. Phys. 26, 218 (1978)]. The "reciprocal sum" has been implemented in three different ways: using an Ewald scheme, a particle mesh Ewald (PME) method, and a multigrid-based approach. We find that even though the use of the McMurchie-Davidson formalism considerably reduces the cost of the calculation with respect to the standard matrix implementation of multipole interactions, the calculation in direct space remains expensive. When most of the calculation is moved to reciprocal space via the PME method, the cost of a calculation where all multipolar interactions (up to hexadecapole-hexadecapole) are included is only about 8.5 times more expensive than a regular AMBER 7 [D. A. Pearlman et al., Comput. Phys. Commun. 91, 1 (1995)] implementation with only charge-charge interactions. The multigrid implementation is slower but shows very promising results for parallelization. It provides a natural way to interface with continuous, Gaussian-based electrostatics in the future. It is
Fortenberry, Ryan C; Huang, Xinchuan; Schwenke, David W; Lee, Timothy J
2014-02-01
In this work, computational procedures are employed to compute the rotational and rovibrational spectra and line lists for H2O, CO2, and SO2. Building on the established use of quartic force fields, MP2 and CCSD(T) Dipole Moment Surfaces (DMSs) are computed for each system of study in order to produce line intensities as well as the transition energies. The computed results exhibit a clear correlation to reference data available in the HITRAN database. Additionally, even though CCSD(T) DMSs produce more accurate intensities as compared to experiment, the use of MP2 DMSs results in reliable line lists that are still comparable to experiment. The use of the less computationally costly MP2 method is beneficial in the study of larger systems where use of CCSD(T) would be more costly. PMID:23692860
Peng, Xiangda; Zhang, Yuebin; Chu, Huiying; Li, Yan; Zhang, Dinglin; Cao, Liaoran; Li, Guohui
2016-06-14
Classical molecular dynamic (MD) simulation of membrane proteins faces significant challenges in accurately reproducing and predicting experimental observables such as ion conductance and permeability due to its incapability of precisely describing the electronic interactions in heterogeneous systems. In this work, the free energy profiles of K(+) and Na(+) permeating through the gramicidin A channel are characterized by using the AMOEBA polarizable force field with a total sampling time of 1 μs. Our results indicated that by explicitly introducing the multipole terms and polarization into the electrostatic potentials, the permeation free energy barrier of K(+) through the gA channel is considerably reduced compared to the overestimated results obtained from the fixed-charge model. Moreover, the estimated maximum conductance, without any corrections, for both K(+) and Na(+) passing through the gA channel are much closer to the experimental results than any classical MD simulations, demonstrating the power of AMOEBA in investigating the membrane proteins. PMID:27171823
Huang, Xinchuan; Taylor, Peter R; Lee, Timothy J
2011-05-19
High levels of theory have been used to compute quartic force fields (QFFs) for the cyclic and linear forms of the C(3)H(3)(+) molecular cation, referred to as c-C(3)H(3)(+) and l-C(3)H(3)(+). Specifically, the singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations, CCSD(T), has been used in conjunction with extrapolation to the one-particle basis set limit, and corrections for scalar relativity and core correlation have been included. The QFFs have been used to compute highly accurate fundamental vibrational frequencies and other spectroscopic constants by use of both vibrational second-order perturbation theory and variational methods to solve the nuclear Schrödinger equation. Agreement between our best computed fundamental vibrational frequencies and recent infrared photodissociation experiments is reasonable for most bands, but there are a few exceptions. Possible sources for the discrepancies are discussed. We determine the energy difference between the cyclic and linear forms of C(3)H(3)(+), obtaining 27.9 kcal/mol at 0 K, which should be the most reliable available. It is expected that the fundamental vibrational frequencies and spectroscopic constants presented here for c-C(3)H(3)(+) and l-C(3)H(3)(+) are the most reliable available for the free gas-phase species, and it is hoped that these will be useful in the assignment of future high-resolution laboratory experiments or astronomical observations. PMID:21510653
NASA Technical Reports Server (NTRS)
Huang, Xinchuan; Taylor, Peter R.; Lee, Timothy J.
2011-01-01
High levels of theory have been used to compute quartic force fields (QFFs) for the cyclic and linear forms of the C H + molecular cation, referred to as c-C H + and I-C H +. Specifically the 33 3333 singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations, CCSD(T), has been used in conjunction with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core correlation have been included. The QFFs have been used to compute highly accurate fundamental vibrational frequencies and other spectroscopic constants using both vibrational 2nd-order perturbation theory and variational methods to solve the nuclear Schroedinger equation. Agreement between our best computed fundamental vibrational frequencies and recent infrared photodissociation experiments is reasonable for most bands, but there are a few exceptions. Possible sources for the discrepancies are discussed. We determine the energy difference between the cyclic and linear forms of C H +, 33 obtaining 27.9 kcal/mol at 0 K, which should be the most reliable available. It is expected that the fundamental vibrational frequencies and spectroscopic constants presented here for c-C H + 33 and I-C H + are the most reliable available for the free gas-phase species and it is hoped that 33 these will be useful in the assignment of future high-resolution laboratory experiments or astronomical observations.
Abramyan, Tigran M.; Snyder, James A.; Yancey, Jeremy A.; Thyparambil, Aby A.; Wei, Yang; Stuart, Steven J.; Latour, Robert A.
2015-01-01
Interfacial force field (IFF) parameters for use with the CHARMM force field have been developed for interactions between peptides and high-density polyethylene (HDPE). Parameterization of the IFF was performed to achieve agreement between experimental and calculated adsorption free energies of small TGTG–X–GTGT host–guest peptides (T = threonine, G = glycine, and X = variable amino-acid residue) on HDPE, with ±0.5 kcal/mol agreement. This IFF parameter set consists of tuned nonbonded parameters (i.e., partial charges and Lennard–Jones parameters) for use with an in-house-modified CHARMM molecular dynamic program that enables the use of an independent set of force field parameters to control molecular behavior at a solid–liquid interface. The R correlation coefficient between the simulated and experimental peptide adsorption free energies increased from 0.00 for the standard CHARMM force field parameters to 0.88 for the tuned IFF parameters. Subsequent studies are planned to apply the tuned IFF parameter set for the simulation of protein adsorption behavior on an HDPE surface for comparison with experimental values of adsorbed protein orientation and conformation. PMID:25818122
Baker, Christopher M.; Lopes, Pedro E. M.; Zhu, Xiao; Roux, Benoît; MacKerell, Alexander D.
2010-01-01
Lennard-Jones (LJ) parameters for a variety of model compounds have previously been optimized within the CHARMM Drude polarizable force field to reproduce accurately pure liquid phase thermodynamic properties as well as additional target data. While the polarizable force field resulting from this optimization procedure has been shown to satisfactorily reproduce a wide range of experimental reference data across numerous series of small molecules, a slight but systematic overestimate of the hydration free energies has also been noted. Here, the reproduction of experimental hydration free energies is greatly improved by the introduction of pair-specific LJ parameters between solute heavy atoms and water oxygen atoms that override the standard LJ parameters obtained from combining rules. The changes are small and a systematic protocol is developed for the optimization of pair-specific LJ parameters and applied to the development of pair-specific LJ parameters for alkanes, alcohols and ethers. The resulting parameters not only yield hydration free energies in good agreement with experimental values, but also provide a framework upon which other pair-specific LJ parameters can be added as new compounds are parametrized within the CHARMM Drude polarizable force field. Detailed analysis of the contributions to the hydration free energies reveals that the dispersion interaction is the main source of the systematic errors in the hydration free energies. This information suggests that the systematic error may result from problems with the LJ combining rules and is combined with analysis of the pair-specific LJ parameters obtained in this work to identify a preliminary improved combining rule. PMID:20401166
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Martin, Jan M. L.; Kwak, Dochan (Technical Monitor)
2002-01-01
The CCSD(T) method has been used to compute a highly accurate quartic force field and fundamental frequencies for all N-14 and N-15 isotopomers of the high energy density material T(sub d)N(sub 4). The computed fundamental frequencies show beyond doubt that the bands observed in a matrix isolation experiment by Radziszewski and coworkers are not due to different isotopomers of T(sub d)N(sub 4). The most sophisticated thermochemical calculations to date yield a N(sub 4) -> 2N(sub 2) heat of reaction of 182.22 +/- 0.5 kcal/mol at 0 K (180.64 +/- 0.5 at 298 K). It is hoped that the data reported herein will aid in the ultimate detection of T(sub d)N(sub 4).
Control Strategies for Accurate Force Generation and Relaxation.
Ohtaka, Chiaki; Fujiwara, Motoko
2016-10-01
Characteristics and motor strategies for force generation and force relaxation were examined using graded tasks during isometric force control. Ten female college students (M age = 20.2 yr., SD = 1.1) were instructed to accurately control the force of isometric elbow flexion using their right arm to match a target force level as quickly as possible. They performed: (1) a generation task, wherein they increased their force from 0% maximum voluntary force to 20% maximum voluntary force (0%-20%), 40% maximum voluntary force (0%-40%), or 60% maximum voluntary force (0%-60%) and (2) and a relaxation task, in which they decreased their force from 60% maximum voluntary force to 40% maximum voluntary force (60%-40%), 20% maximum voluntary force (60%-20%), or to 0% maximum voluntary force (60%-0%). Produced force parameters of point of accuracy (force level, error), quickness (reaction time, adjustment time, rate of force development), and strategy (force wave, rate of force development) were analyzed. Errors of force relaxation were all greater, and reaction times shorter, than those of force generation. Adjustment time depended on the magnitude of force and peak rates of force development and force relaxation differed. Controlled relaxation of force is more difficult with low magnitude of force control. PMID:27555365
Foley, B. Lachele; Tessier, Matthew B.; Woods, Robert J.
2014-01-01
Carbohydrates present a special set of challenges to the generation of force fields. First, the tertiary structures of monosaccharides are complex merely by virtue of their exceptionally high number of chiral centers. In addition, their electronic characteristics lead to molecular geometries and electrostatic landscapes that can be challenging to predict and model. The monosaccharide units can also interconnect in many ways, resulting in a large number of possible oligosaccharides and polysaccharides, both linear and branched. These larger structures contain a number of rotatable bonds, meaning they potentially sample an enormous conformational space. This article briefly reviews the history of carbohydrate force fields, examining and comparing their challenges, forms, philosophies, and development strategies. Then it presents a survey of recent uses of these force fields, noting trends, strengths, deficiencies, and possible directions for future expansion. PMID:25530813
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Martin, Jan M. L.; Dateo, Christopher E.; Taylor, Peter R.
1995-01-01
The XCN and XNC (X = F, Cl) isomers have been investigated using the CCSD(T) method in conjunction with correlation consistent basis sets. Equilibrium geometries, harmonic frequencies, anharmonic constants, fundamental frequencies, and heats of formation have been evaluated. Agreement with experiment for the fundamental frequencies is very good, even for nu(sub 2), for CICN, which is subject to a strong Fermi resonance with 2nu(sub 3). It is also shown that a second-order perturbation theory approach to solving the nuclear Schroedinger equation gives results in excellent agreement with essentially exact variational calculations. This is true even for nu(sub 2) of ClCN, provided that near-singular terms are eliminated from the perturbation theory formulas and the appropriate Fermi interaction energy matrix is then diagonalized. A band at 615/cm, tentatively assigned as the Cl-N stretch in ClNC in matrix isolation experiments, is shown not to be due to ClNC. Accurate atomization energies are determined and are used to evaluate accurate heats of formation (3.1 +/- 1.5, 33.2 +/- 1.5, 72.6 +/- 1.5, and 75.9 +/- 1.5 kcal/mol for FCN, ClCN, FNC, and ClNC, respectively). It is expected that the theoretical heats of formation for FCN, FNC, and ClNC are the most accurate available.
A general, accurate procedure for calculating molecular interaction force.
Yang, Pinghai; Qian, Xiaoping
2009-09-15
The determination of molecular interaction forces, e.g., van der Waals force, between macroscopic bodies is of fundamental importance for understanding sintering, adhesion and fracture processes. In this paper, we develop an accurate, general procedure for van der Waals force calculation. This approach extends a surface formulation that converts a six-dimensional (6D) volume integral into a 4D surface integral for the force calculation. It uses non-uniform rational B-spline (NURBS) surfaces to represent object surfaces. Surface integrals are then done on the parametric domain of the NURBS surfaces. It has combined advantages of NURBS surface representation and surface formulation, including (1) molecular interactions between arbitrary-shaped objects can be represented and evaluated by the NURBS model further common geometries such as spheres, cones, planes can be represented exactly and interaction forces are thus calculated accurately; (2) calculation efficiency is improved by converting the volume integral to the surface integral. This approach is implemented and validated via its comparison with analytical solutions for simple geometries. Calculation of van der Waals force between complex geometries with surface roughness is also demonstrated. A tutorial on the NURBS approach is given in Appendix A. PMID:19596335
Polarization effects in molecular mechanical force fields.
Cieplak, Piotr; Dupradeau, François-Yves; Duan, Yong; Wang, Junmei
2009-08-19
The focus here is on incorporating electronic polarization into classical molecular mechanical force fields used for macromolecular simulations. First, we briefly examine currently used molecular mechanical force fields and the current status of intermolecular forces as viewed by quantum mechanical approaches. Next, we demonstrate how some components of quantum mechanical energy are effectively incorporated into classical molecular mechanical force fields. Finally, we assess the modeling methods of one such energy component-polarization energy-and present an overview of polarizable force fields and their current applications. Incorporating polarization effects into current force fields paves the way to developing potentially more accurate, though more complex, parameterizations that can be used for more realistic molecular simulations. PMID:21828594
Polarization effects in molecular mechanical force fields
Cieplak, Piotr; Dupradeau, François-Yves; Duan, Yong; Wang, Junmei
2014-01-01
The focus here is on incorporating electronic polarization into classical molecular mechanical force fields used for macromolecular simulations. First, we briefly examine currently used molecular mechanical force fields and the current status of intermolecular forces as viewed by quantum mechanical approaches. Next, we demonstrate how some components of quantum mechanical energy are effectively incorporated into classical molecular mechanical force fields. Finally, we assess the modeling methods of one such energy component—polarization energy—and present an overview of polarizable force fields and their current applications. Incorporating polarization effects into current force fields paves the way to developing potentially more accurate, though more complex, parameterizations that can be used for more realistic molecular simulations. PMID:21828594
Energy Science and Technology Software Center (ESTSC)
2015-05-27
ParFit is a flexible and extendable framework and library of classes for fitting force-field parameters to data from high-level ab-initio calculations on the basis of deterministic and stochastic algorithms. Currently, the code is fitting MM3 and Merck force-field parameters but could easily extend to other force-field types.
Motor equivalence during multi-finger accurate force production
Mattos, Daniela; Schöner, Gregor; Zatsiorsky, Vladimir M.; Latash, Mark L.
2014-01-01
We explored stability of multi-finger cyclical accurate force production action by analysis of responses to small perturbations applied to one of the fingers and inter-cycle analysis of variance. Healthy subjects performed two versions of the cyclical task, with and without an explicit target. The “inverse piano” apparatus was used to lift/lower a finger by 1 cm over 0.5 s; the subjects were always instructed to perform the task as accurate as they could at all times. Deviations in the spaces of finger forces and modes (hypothetical commands to individual fingers) were quantified in directions that did not change total force (motor equivalent) and in directions that changed the total force (non-motor equivalent). Motor equivalent deviations started immediately with the perturbation and increased progressively with time. After a sequence of lifting-lowering perturbations leading to the initial conditions, motor equivalent deviations were dominating. These phenomena were less pronounced for analysis performed with respect to the total moment of force with respect to an axis parallel to the forearm/hand. Analysis of inter-cycle variance showed consistently higher variance in a subspace that did not change the total force as compared to the variance that affected total force. We interpret the results as reflections of task-specific stability of the redundant multi-finger system. Large motor equivalent deviations suggest that reactions of the neuromotor system to a perturbation involve large changes of neural commands that do not affect salient performance variables, even during actions with the purpose to correct those salient variables. Consistency of the analyses of motor equivalence and variance analysis provides additional support for the idea of task-specific stability ensured at a neural level. PMID:25344311
Accurate Satellite-Derived Estimates of Tropospheric Ozone Radiative Forcing
NASA Technical Reports Server (NTRS)
Joiner, Joanna; Schoeberl, Mark R.; Vasilkov, Alexander P.; Oreopoulos, Lazaros; Platnick, Steven; Livesey, Nathaniel J.; Levelt, Pieternel F.
2008-01-01
Estimates of the radiative forcing due to anthropogenically-produced tropospheric O3 are derived primarily from models. Here, we use tropospheric ozone and cloud data from several instruments in the A-train constellation of satellites as well as information from the GEOS-5 Data Assimilation System to accurately estimate the instantaneous radiative forcing from tropospheric O3 for January and July 2005. We improve upon previous estimates of tropospheric ozone mixing ratios from a residual approach using the NASA Earth Observing System (EOS) Aura Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) by incorporating cloud pressure information from OMI. Since we cannot distinguish between natural and anthropogenic sources with the satellite data, our estimates reflect the total forcing due to tropospheric O3. We focus specifically on the magnitude and spatial structure of the cloud effect on both the shortand long-wave radiative forcing. The estimates presented here can be used to validate present day O3 radiative forcing produced by models.
Accurate forced-choice recognition without awareness of memory retrieval.
Voss, Joel L; Baym, Carol L; Paller, Ken A
2008-06-01
Recognition confidence and the explicit awareness of memory retrieval commonly accompany accurate responding in recognition tests. Memory performance in recognition tests is widely assumed to measure explicit memory, but the generality of this assumption is questionable. Indeed, whether recognition in nonhumans is always supported by explicit memory is highly controversial. Here we identified circumstances wherein highly accurate recognition was unaccompanied by hallmark features of explicit memory. When memory for kaleidoscopes was tested using a two-alternative forced-choice recognition test with similar foils, recognition was enhanced by an attentional manipulation at encoding known to degrade explicit memory. Moreover, explicit recognition was most accurate when the awareness of retrieval was absent. These dissociations between accuracy and phenomenological features of explicit memory are consistent with the notion that correct responding resulted from experience-dependent enhancements of perceptual fluency with specific stimuli--the putative mechanism for perceptual priming effects in implicit memory tests. This mechanism may contribute to recognition performance in a variety of frequently-employed testing circumstances. Our results thus argue for a novel view of recognition, in that analyses of its neurocognitive foundations must take into account the potential for both (1) recognition mechanisms allied with implicit memory and (2) recognition mechanisms allied with explicit memory. PMID:18519546
Chung, Koo-Hyun; Pratt, Jon R; Reitsma, Mark G
2010-01-19
The colloidal probe technique for atomic force microscopy (AFM) has allowed the investigation of an extensive range of surface force phenomena, including the measurement of frictional (lateral) forces between numerous materials. The quantitative accuracy of such friction measurements is often debated, in part due to a lack of confidence in existing calibration strategies. Here we compare three in situ AFM lateral force calibration techniques using a single colloidal probe, seeking to establish a foundation for quantitative measurement by linking these techniques to accurate force references available at the National Institute of Standards and Technology. We introduce a procedure for calibrating the AFM lateral force response to known electrostatic forces applied directly to the conductive colloidal probe. In a second procedure, we apply known force directly to the colloidal probe using a precalibrated piezo-resistive reference cantilever. We found agreement between these direct methods on the order of 2% (within random uncertainty for both measurements). In a third procedure, we performed a displacement-based calibration using the piezo-resistive reference cantilever as a stiffness reference artifact. The method demonstrated agreement on the order of 7% with the direct force methods, with the difference attributed to an expected systematic uncertainty, caused by in-plane deflection in the cantilever during loading. The comparison establishes the existing limits of instrument accuracy and sets down a basis for selection criteria for materials and methods in colloidal probe friction (lateral) force measurements via atomic force microscopy. PMID:19827782
Huang, Xinchuan; Valeev, Edward F; Lee, Timothy J
2010-12-28
One-particle basis set extrapolation is compared with one of the new R12 methods for computing highly accurate quartic force fields (QFFs) and spectroscopic data, including molecular structures, rotational constants, and vibrational frequencies for the H(2)O, N(2)H(+), NO(2)(+), and C(2)H(2) molecules. In general, agreement between the spectroscopic data computed from the best R12 and basis set extrapolation methods is very good with the exception of a few parameters for N(2)H(+) where it is concluded that basis set extrapolation is still preferred. The differences for H(2)O and NO(2)(+) are small and it is concluded that the QFFs from both approaches are more or less equivalent in accuracy. For C(2)H(2), however, a known one-particle basis set deficiency for C-C multiple bonds significantly degrades the quality of results obtained from basis set extrapolation and in this case the R12 approach is clearly preferred over one-particle basis set extrapolation. The R12 approach used in the present study was modified in order to obtain high precision electronic energies, which are needed when computing a QFF. We also investigated including core-correlation explicitly in the R12 calculations, but conclude that current approaches are lacking. Hence core-correlation is computed as a correction using conventional methods. Considering the results for all four molecules, it is concluded that R12 methods will soon replace basis set extrapolation approaches for high accuracy electronic structure applications such as computing QFFs and spectroscopic data for comparison to high-resolution laboratory or astronomical observations, provided one uses a robust R12 method as we have done here. The specific R12 method used in the present study, CCSD(T)(R12), incorporated a reformulation of one intermediate matrix in order to attain machine precision in the electronic energies. Final QFFs for N(2)H(+) and NO(2)(+) were computed, including basis set extrapolation, core-correlation, scalar
NASA Astrophysics Data System (ADS)
Shavezipur, M.; Li, G. H.; Laboriante, I.; Gou, W. J.; Carraro, C.; Maboudian, R.
2011-11-01
This paper reports on accurate analysis of adhesion force between polysilicon-polysilicon surfaces in micro-/nanoelectromechanical systems (M/NEMS). The measurement is carried out using double-clamped beams. Electrostatic actuation and structural restoring force are exploited to respectively initiate and terminate the contact between the two surfaces under investigation. The adhesion force is obtained by balancing the electrostatic and mechanical forces acting on the beam just before the separation of the two surfaces. Different finite element models are developed to simulate the coupled-field multiphysics problem. The effects of fringing field in the electrostatic domain and geometric nonlinearity and residual stress in the structural domain are taken into consideration. Moreover, the beam stiffness is directly obtained for the case of combined loading (electrostatic and adhesion). Therefore, the overall electrostatic and structural forces used to extract the actual adhesion force from measured data are determined with high accuracy leading to accurate values for the adhesion force. The finite element simulations presented in this paper are not limited to adhesion force measurement and can be used to design or characterize electrostatically actuated devices such as MEM tunable capacitors and micromirrors, RF switches and M/NEM relays.
Accurate Extraction of Electrostatic Force by a Voltage-Pulse Force Spectroscopy.
Inami, Eiichi; Sugimoto, Yoshiaki
2015-06-19
The classification of interaction forces between two approaching bodies is important in a wide range of research fields. Here, we propose a method to unambiguously extract the electrostatic force (F(ele)), which is one of the most significant forces. This method is based on the measurement of the energy dissipation under applied voltage pulse between an atomic force microscopy (AFM) tip and sample. It allowed us to obtain F(ele) as a function of the tip-sample distance and voltage including the distance-independent part, to which conventional AFM is insensitive. The obtained F(ele) curves nicely fit the analytical model, enabling estimation of the geometry of the tip. The distance-dependent contact potential difference could also be correctly obtained by the measured F(ele), opening an alternative route to quantitative Kelvin probe force microscopy. PMID:26196989
Accurate Extraction of Electrostatic Force by a Voltage-Pulse Force Spectroscopy
NASA Astrophysics Data System (ADS)
Inami, Eiichi; Sugimoto, Yoshiaki
2015-06-01
The classification of interaction forces between two approaching bodies is important in a wide range of research fields. Here, we propose a method to unambiguously extract the electrostatic force (Fele ), which is one of the most significant forces. This method is based on the measurement of the energy dissipation under applied voltage pulse between an atomic force microscopy (AFM) tip and sample. It allowed us to obtain Fele as a function of the tip-sample distance and voltage including the distance-independent part, to which conventional AFM is insensitive. The obtained Fele curves nicely fit the analytical model, enabling estimation of the geometry of the tip. The distance-dependent contact potential difference could also be correctly obtained by the measured Fele, opening an alternative route to quantitative Kelvin probe force microscopy.
Accurate force spectroscopy in tapping mode atomic force microscopy in liquids
NASA Astrophysics Data System (ADS)
Xu, Xin; Melcher, John; Raman, Arvind
2010-01-01
Existing force spectroscopy methods in tapping mode atomic force microscopy (AFM) such as higher harmonic inversion [M. Stark, R. W. Stark, W. M. Heckl, and R. Guckenberger, Proc. Natl. Acad. Sci. U.S.A. 99, 8473 (2002)] or scanning probe acceleration microscopy [J. Legleiter, M. Park, B. Cusick, and T. Kowalewski, Proc. Natl. Acad. Sci. U.S.A. 103, 4813 (2006)] or integral relations [M. Lee and W. Jhe, Phys. Rev. Lett. 97, 036104 (2006); S. Hu and A. Raman, Nanotechnology 19, 375704 (2008); H. Hölscher, Appl. Phys. Lett. 89, 123109 (2006); A. J. Katan, Nanotechnology 20, 165703 (2009)] require and assume as an observable the tip dynamics in a single eigenmode of the oscillating microcantilever. We demonstrate that this assumption can distort significantly the extracted tip-sample interaction forces when applied to tapping mode AFM with soft cantilevers in liquid environments. This exception is due to the fact that under these conditions the second eigenmode is momentarily excited and the observed tip dynamics clearly contains contributions from the fundamental and second eigenmodes. To alleviate this problem, a simple experimental method is proposed to screen the second eigenmode contributions in the observed tip deflection signal to allow accurate tip-sample force reconstruction in liquids. The method is implemented experimentally to reconstruct interaction forces on polymer, bacteriorhodopsin membrane, and mica samples in buffer solutions.
NASA Astrophysics Data System (ADS)
Inami, Eiichi; Sugimoto, Yoshiaki
2016-08-01
We propose a new method of extracting electrostatic force. The technique is based on frequency modulation atomic force microscopy (FM-AFM) combined with a voltage pulse. In this method, the work that the electrostatic field does on the oscillating tip is measured through the cantilever energy dissipation. This allows us to directly extract capacitive forces including the longer range part, to which the conventional FM-AFM is insensitive. The distance-dependent contact potential difference, which is modulated by local charges distributed on the surfaces of the tip and/or sample, could also be correctly obtained. In the absence of local charges, our method can perfectly reproduce the electrostatic force as a function of the distance and the bias voltage. Furthermore, we demonstrate that the system serves as a sensitive sensor enabling us to check the existence of the local charges such as trapped charges and patch charges.
Devereux, Mike; Gresh, Nohad; Piquemal, Jean-Philip; Meuwly, Markus
2014-08-01
A supervised, semiautomated approach to force field parameter fitting is described and applied to the SIBFA polarizable force field. The I-NoLLS interactive, nonlinear least squares fitting program is used as an engine for parameter refinement while keeping parameter values within a physical range. Interactive fitting is shown to avoid many of the stability problems that frequently afflict highly correlated, nonlinear fitting problems occurring in force field parametrizations. The method is used to obtain parameters for the H2O, formamide, and imidazole molecular fragments and their complexes with the Mg(2+) cation. Reference data obtained from ab initio calculations using an auc-cc-pVTZ basis set exploit advances in modern computer hardware to provide a more accurate parametrization of SIBFA than has previously been available. PMID:24965869
Accurate computation of surface stresses and forces with immersed boundary methods
NASA Astrophysics Data System (ADS)
Goza, Andres; Liska, Sebastian; Morley, Benjamin; Colonius, Tim
2016-09-01
Many immersed boundary methods solve for surface stresses that impose the velocity boundary conditions on an immersed body. These surface stresses may contain spurious oscillations that make them ill-suited for representing the physical surface stresses on the body. Moreover, these inaccurate stresses often lead to unphysical oscillations in the history of integrated surface forces such as the coefficient of lift. While the errors in the surface stresses and forces do not necessarily affect the convergence of the velocity field, it is desirable, especially in fluid-structure interaction problems, to obtain smooth and convergent stress distributions on the surface. To this end, we show that the equation for the surface stresses is an integral equation of the first kind whose ill-posedness is the source of spurious oscillations in the stresses. We also demonstrate that for sufficiently smooth delta functions, the oscillations may be filtered out to obtain physically accurate surface stresses. The filtering is applied as a post-processing procedure, so that the convergence of the velocity field is unaffected. We demonstrate the efficacy of the method by computing stresses and forces that converge to the physical stresses and forces for several test problems.
TOPICAL REVIEW: Polarization effects in molecular mechanical force fields
NASA Astrophysics Data System (ADS)
Cieplak, Piotr; Dupradeau, François-Yves; Duan, Yong; Wang, Junmei
2009-08-01
The focus here is on incorporating electronic polarization into classical molecular mechanical force fields used for macromolecular simulations. First, we briefly examine currently used molecular mechanical force fields and the current status of intermolecular forces as viewed by quantum mechanical approaches. Next, we demonstrate how some components of quantum mechanical energy are effectively incorporated into classical molecular mechanical force fields. Finally, we assess the modeling methods of one such energy component—polarization energy—and present an overview of polarizable force fields and their current applications. Incorporating polarization effects into current force fields paves the way to developing potentially more accurate, though more complex, parameterizations that can be used for more realistic molecular simulations.
Accurate wavelength calibration method for flat-field grating spectrometers.
Du, Xuewei; Li, Chaoyang; Xu, Zhe; Wang, Qiuping
2011-09-01
A portable spectrometer prototype is built to study wavelength calibration for flat-field grating spectrometers. An accurate calibration method called parameter fitting is presented. Both optical and structural parameters of the spectrometer are included in the wavelength calibration model, which accurately describes the relationship between wavelength and pixel position. Along with higher calibration accuracy, the proposed calibration method can provide information about errors in the installation of the optical components, which will be helpful for spectrometer alignment. PMID:21929865
Light Field Imaging Based Accurate Image Specular Highlight Removal
Wang, Haoqian; Xu, Chenxue; Wang, Xingzheng; Zhang, Yongbing; Peng, Bo
2016-01-01
Specular reflection removal is indispensable to many computer vision tasks. However, most existing methods fail or degrade in complex real scenarios for their individual drawbacks. Benefiting from the light field imaging technology, this paper proposes a novel and accurate approach to remove specularity and improve image quality. We first capture images with specularity by the light field camera (Lytro ILLUM). After accurately estimating the image depth, a simple and concise threshold strategy is adopted to cluster the specular pixels into “unsaturated” and “saturated” category. Finally, a color variance analysis of multiple views and a local color refinement are individually conducted on the two categories to recover diffuse color information. Experimental evaluation by comparison with existed methods based on our light field dataset together with Stanford light field archive verifies the effectiveness of our proposed algorithm. PMID:27253083
Improving an all-atom force field.
Mohanty, Sandipan; Hansmann, U H E
2007-07-01
Experimentally well-characterized proteins that are small enough to be computationally tractable provide useful information for refining existing all-atom force fields. This is used by us for reparametrizing a recently developed all-atom force field. Relying on high statistics parallel tempering simulations of a designed 20 residue beta-sheet peptide, we propose incremental changes that improve the force field's range of applicability. PMID:17677516
A molecular mechanics force field for lignin
Petridis, Loukas; Smith, Jeremy C
2009-02-01
A CHARMM molecular mechanics force field for lignin is derived. Parameterization is based on reproducing quantum mechanical data of model compounds. Partial atomic charges are derived using the RESP electrostatic potential fitting method supplemented by the examination of methoxybenzene:water interactions. Dihedral parameters are optimized by fitting to critical rotational potentials and bonded parameters are obtained by optimizing vibrational frequencies and normal modes. Finally, the force field is validated by performing a molecular dynamics simulation of a crystal of a lignin fragment molecule and comparing simulation-derived structural features with experimental results. Together with the existing force field for polysaccharides, this lignin force field will enable full simulations of lignocellulose.
Accurate rotor loads prediction using the FLAP (Force and Loads Analysis Program) dynamics code
Wright, A.D.; Thresher, R.W.
1987-10-01
Accurately predicting wind turbine blade loads and response is very important in predicting the fatigue life of wind turbines. There is a clear need in the wind turbine community for validated and user-friendly structural dynamics codes for predicting blade loads and response. At the Solar Energy Research Institute (SERI), a Force and Loads Analysis Program (FLAP) has been refined and validated and is ready for general use. Currently, FLAP is operational on an IBM-PC compatible computer and can be used to analyze both rigid- and teetering-hub configurations. The results of this paper show that FLAP can be used to accurately predict the deterministic loads for rigid-hub rotors. This paper compares analytical predictions to field test measurements for a three-bladed, upwind turbine with a rigid-hub configuration. The deterministic loads predicted by FLAP are compared with 10-min azimuth averages of blade root flapwise bending moments for different wind speeds. 6 refs., 12 figs., 3 tabs.
Accurate noncontact calibration of colloidal probe sensitivities in atomic force microscopy
Chung, Koo-Hyun; Shaw, Gordon A.; Pratt, Jon R.
2009-06-15
The absolute force sensitivities of colloidal probes comprised of atomic force microscope, or AFM, cantilevers with microspheres attached to their distal ends are measured. The force sensitivities are calibrated through reference to accurate electrostatic forces, the realizations of which are described in detail. Furthermore, the absolute accuracy of a common AFM force calibration scheme, known as the thermal noise method, is evaluated. It is demonstrated that the thermal noise method can be applied with great success to colloidal probe calibration in air and in liquid to yield force measurements with relative standard uncertainties below 5%. Techniques to combine the electrostatics-based determination of the AFM force sensitivity with measurements of the colloidal probe's thermal noise spectrum to compute noncontact estimates of the displacement sensitivity and spring constant are also developed.
Force field development from first principles for materials design
NASA Astrophysics Data System (ADS)
Chan, Maria; Kinaci, Alper; Narayanan, Badri; Sen, Fatih; Gray, Stephen; Davis, Michael; Sankaranaryanan, Subramanian
2015-03-01
The ability to perform accurate calculations efficiently is crucial for computational materials design. In this talk, we will discuss a stream-lined approach to force field development using first principles density functional theory training data and machine learning algorithms. We will also discuss the validation of this approach on precious metal nanoparticles.
Force Field Development for Lipid Membrane Simulations.
Lyubartsev, Alexander P; Rabinovich, Alexander L
2016-10-01
With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg. PMID:26766518
Static and dynamical Meissner force fields
NASA Technical Reports Server (NTRS)
Weinberger, B. R.; Lynds, L.; Hull, J. R.; Mulcahy, T. M.
1991-01-01
The coupling between copper-based high temperature superconductors (HTS) and magnets is represented by a force field. Zero-field cooled experiments were performed with several forms of superconductors: 1) cold-pressed sintered cylindrical disks; 2) small particles fixed in epoxy polymers; and 3) small particles suspended in hydrocarbon waxes. Using magnets with axial field symmetries, direct spatial force measurements in the range of 0.1 to 10(exp 4) dynes were performed with an analytical balance and force constants were obtained from mechanical vibrational resonances. Force constants increase dramatically with decreasing spatial displacement. The force field displays a strong temperature dependence between 20 and 90 K and decreases exponentially with increasing distance of separation. Distinct slope changes suggest the presence of B-field and temperature-activated processes that define the forces. Hysteresis measurements indicated that the magnitude of force scales roughly with the volume fraction of HTS in composite structures. Thus, the net force resulting from the field interaction appears to arise from regions as small or smaller than the grain size and does not depend on contiguous electron transport over large areas. Results of these experiments are discussed.
Age-related changes in multifinger synergies in accurate moment of force production tasks
Olafsdottir, Halla; Zhang, Wei; Zatsiorsky, Vladimir M.; Latash, Mark L.
2010-01-01
The purpose of this investigation was to document and quantify age-related differences in the coordination of fingers during a task that required production of an accurate time profile of the total moment of force by the four fingers of a hand. We hypothesized that elderly subjects would show a decreased ability to stabilize a time profile of the total moment of force, leading to larger indexes of moment variability compared with young subjects. The subjects followed a trapezoidal template on a computer screen by producing a time profile of the total moment of force while pressing down on force sensors with the four fingers of the right (dominant) hand. To quantify synergies, we used the framework of the uncontrolled manifold hypothesis. The elderly subjects produced larger total force, larger variance of both total force and total moment of force, and larger involvement of fingers that produced moment of force against the required moment direction (antagonist moment). This was particularly prominent during supination efforts. Young subjects showed covariation of commands to fingers across trials that stabilized the moment of total force (moment-stabilizing synergy), while elderly subjects failed to do so. Both subject groups showed similar indexes of covariation of commands to the fingers that stabilized the time profile of the total force. The lack of moment-stabilizing synergies may be causally related to the documented impairment of hand function with age. PMID:17204576
Alternating Magnetic Field Forces for Satellite Formation Flying
NASA Technical Reports Server (NTRS)
Youngquist, Robert C.; Nurge, Mark A.; Starr, Stnaley O.
2012-01-01
Selected future space missions, such as large aperture telescopes and multi-component interferometers, will require the precise positioning of a number of isolated satellites, yet many of the suggested approaches for providing satellites positioning forces have serious limitations. In this paper we propose a new approach, capable of providing both position and orientation forces, that resolves or alleviates many of these problems. We show that by using alternating fields and currents that finely-controlled forces can be induced on the satellites, which can be individually selected through frequency allocation. We also show, through analysis and experiment, that near field operation is feasible and can provide sufficient force and the necessary degrees of freedom to accurately position and orient small satellites relative to one another. In particular, the case of a telescope with a large number of free mirrors is developed to provide an example of the concept. We. also discuss the far field extension of this concept.
NASA Astrophysics Data System (ADS)
Mead, A. J.; Heymans, C.; Lombriser, L.; Peacock, J. A.; Steele, O. I.; Winther, H. A.
2016-06-01
We present an accurate non-linear matter power spectrum prediction scheme for a variety of extensions to the standard cosmological paradigm, which uses the tuned halo model previously developed in Mead et al. We consider dark energy models that are both minimally and non-minimally coupled, massive neutrinos and modified gravitational forces with chameleon and Vainshtein screening mechanisms. In all cases, we compare halo-model power spectra to measurements from high-resolution simulations. We show that the tuned halo-model method can predict the non-linear matter power spectrum measured from simulations of parametrized w(a) dark energy models at the few per cent level for k < 10 h Mpc-1, and we present theoretically motivated extensions to cover non-minimally coupled scalar fields, massive neutrinos and Vainshtein screened modified gravity models that result in few per cent accurate power spectra for k < 10 h Mpc-1. For chameleon screened models, we achieve only 10 per cent accuracy for the same range of scales. Finally, we use our halo model to investigate degeneracies between different extensions to the standard cosmological model, finding that the impact of baryonic feedback on the non-linear matter power spectrum can be considered independently of modified gravity or massive neutrino extensions. In contrast, considering the impact of modified gravity and massive neutrinos independently results in biased estimates of power at the level of 5 per cent at scales k > 0.5 h Mpc-1. An updated version of our publicly available HMCODE can be found at https://github.com/alexander-mead/hmcode.
Comparing molecular dynamics force fields in the essential subspace.
Martín-García, Fernando; Papaleo, Elena; Gomez-Puertas, Paulino; Boomsma, Wouter; Lindorff-Larsen, Kresten
2015-01-01
The continued development and utility of molecular dynamics simulations requires improvements in both the physical models used (force fields) and in our ability to sample the Boltzmann distribution of these models. Recent developments in both areas have made available multi-microsecond simulations of two proteins, ubiquitin and Protein G, using a number of different force fields. Although these force fields mostly share a common mathematical form, they differ in their parameters and in the philosophy by which these were derived, and previous analyses showed varying levels of agreement with experimental NMR data. To complement the comparison to experiments, we have performed a structural analysis of and comparison between these simulations, thereby providing insight into the relationship between force-field parameterization, the resulting ensemble of conformations and the agreement with experiments. In particular, our results show that, at a coarse level, many of the motional properties are preserved across several, though not all, force fields. At a finer level of detail, however, there are distinct differences in both the structure and dynamics of the two proteins, which can, together with comparison with experimental data, help to select force fields for simulations of proteins. A noteworthy observation is that force fields that have been reparameterized and improved to provide a more accurate energetic description of the balance between helical and coil structures are difficult to distinguish from their "unbalanced" counterparts in these simulations. This observation implies that simulations of stable, folded proteins, even those reaching 10 microseconds in length, may provide relatively little information that can be used to modify torsion parameters to achieve an accurate balance between different secondary structural elements. PMID:25811178
Accurate calculation of field and carrier distributions in doped semiconductors
NASA Astrophysics Data System (ADS)
Yang, Wenji; Tang, Jianping; Yu, Hongchun; Wang, Yanguo
2012-06-01
We use the numerical squeezing algorithm(NSA) combined with the shooting method to accurately calculate the built-in fields and carrier distributions in doped silicon films (SFs) in the micron and sub-micron thickness range and results are presented in graphical form for variety of doping profiles under different boundary conditions. As a complementary approach, we also present the methods and the results of the inverse problem (IVP) - finding out the doping profile in the SFs for given field distribution. The solution of the IVP provides us the approach to arbitrarily design field distribution in SFs - which is very important for low dimensional (LD) systems and device designing. Further more, the solution of the IVP is both direct and much easy for all the one-, two-, and three-dimensional semiconductor systems. With current efforts focused on the LD physics, knowing of the field and carrier distribution details in the LD systems will facilitate further researches on other aspects and hence the current work provides a platform for those researches.
Equifinality and its violations in a redundant system: multifinger accurate force production
Wilhelm, Luke; Zatsiorsky, Vladimir M.
2013-01-01
We explored a hypothesis that transient perturbations applied to a redundant system result in equifinality in the space of task-related performance variables but not in the space of elemental variables. The subjects pressed with four fingers and produced an accurate constant total force level. The “inverse piano” device was used to lift and lower one of the fingers smoothly. The subjects were instructed “not to intervene voluntarily” with possible force changes. Analysis was performed in spaces of finger forces and finger modes (hypothetical neural commands to fingers) as elemental variables. Lifting a finger led to an increase in its force and a decrease in the forces of the other three fingers; the total force increased. Lowering the finger back led to a drop in the force of the perturbed finger. At the final state, the sum of the variances of finger forces/modes computed across repetitive trials was significantly higher than the variance of the total force/mode. Most variance of the individual finger force/mode changes between the preperturbation and postperturbation states was compatible with constant total force. We conclude that a transient perturbation applied to a redundant system leads to relatively small variance in the task-related performance variable (equifinality), whereas in the space of elemental variables much more variance occurs that does not lead to total force changes. We interpret the results within a general theoretical scheme that incorporates the ideas of hierarchically organized control, control with referent configurations, synergic control, and the uncontrolled manifold hypothesis. PMID:23904497
Three-Dimensional Force Field Spectroscopy
NASA Astrophysics Data System (ADS)
Schwarz, Alexander; Hölscher, Hendrik; Langkat, S. M.; Wiesendanger, R.
2003-12-01
A method is presented that utilizes the frequency modulation technique in ultra-high vacuum to measure the tip-sample force field in all three dimensions with atomic resolution. It is based on a systematic procedure to record frequency shift versus distance curves. After their conversion into the tip-surface potential landscape the complete force field in all three dimensions can be calculated. Experimental results obtained in the non-contact regime on NiO(001) with an iron-coated silicon tip are presented to demonstrate that interatomic vertical and lateral forces can be determined and assigned to specific sites within the surface unit cell.
NASA Astrophysics Data System (ADS)
Hedrick, A. R.; Marks, D. G.; Winstral, A. H.; Marshall, H. P.
2014-12-01
The ability to forecast snow water equivalent, or SWE, in mountain catchments would benefit many different communities ranging from avalanche hazard mitigation to water resource management. Historical model runs of Isnobal, the physically based energy balance snow model, have been produced over the 2150 km2 Boise River Basin for water years 2012 - 2014 at 100-meter resolution. Spatially distributed forcing parameters such as precipitation, wind, and relative humidity are generated from automated weather stations located throughout the watershed, and are supplied to Isnobal at hourly timesteps. Similarly, the Weather Research & Forecasting (WRF) Model provides hourly predictions of the same forcing parameters from an atmospheric physics perspective. This work aims to quantitatively compare WRF model output to the spatial meteorologic fields developed to force Isnobal, with the hopes of eventually using WRF predictions to create accurate hourly forecasts of SWE over a large mountainous basin.
Brief: Field measurements of casing tension forces
Quigley, M.S.; Lewis, D.B.; Boswell, R.S.
1995-02-01
Tension forces acting on individual casing joints were accurately measured during installation of 10,158 ft of 9 5/8-in. {times} 47-lbm/ft casing and 11,960 ft of 11 7/8-in. {times} 71.8-lbm/ft casing. A unique casing load table (CLT) weighed the casing string after the addition of each casing joint. Strain gauges attached inside the pin ends of instrumented casing joints (ICJ`s) directly measured tension force on those joints. A high-speed computer data-acquisition system (DAS) automatically recorded data from all the sensors. Several casing joints were intentionally subjected to extreme deceleration to determine upper limits for dynamic tension forces. Data from these tests clearly show effects of wellbore friction and casing handling conditions. In every case, tension forces in the casing during maximum deceleration were considerably less than expected. In some cases, the highest tension forces occurred when the casing lifted out of the slips. Peak tension forces caused by setting the casing slips were typically no more than 5% greater than tension forces in the casing at rest. This dynamic amplification was far less than the 60% value used in the previous casing design method. Reducing the safety factor for installation loads has permitted use of lighter, less-expensive casing than dictated by older design criteria.
Lipid14: The Amber Lipid Force Field.
Dickson, Callum J; Madej, Benjamin D; Skjevik, Age A; Betz, Robin M; Teigen, Knut; Gould, Ian R; Walker, Ross C
2014-02-11
The AMBER lipid force field has been updated to create Lipid14, allowing tensionless simulation of a number of lipid types with the AMBER MD package. The modular nature of this force field allows numerous combinations of head and tail groups to create different lipid types, enabling the easy insertion of new lipid species. The Lennard-Jones and torsion parameters of both the head and tail groups have been revised and updated partial charges calculated. The force field has been validated by simulating bilayers of six different lipid types for a total of 0.5 μs each without applying a surface tension; with favorable comparison to experiment for properties such as area per lipid, volume per lipid, bilayer thickness, NMR order parameters, scattering data, and lipid lateral diffusion. As the derivation of this force field is consistent with the AMBER development philosophy, Lipid14 is compatible with the AMBER protein, nucleic acid, carbohydrate, and small molecule force fields. PMID:24803855
Controlling Casimir force via coherent driving field
NASA Astrophysics Data System (ADS)
Ahmad, Rashid; Abbas, Muqaddar; Ahmad, Iftikhar; Qamar, Sajid
2016-04-01
A four level atom-field configuration is used to investigate the coherent control of Casimir force between two identical plates made up of chiral atomic media and separated by vacuum of width d. The electromagnetic chirality-induced negative refraction is obtained via atomic coherence. The behavior of Casimir force is investigated using Casimir-Lifshitz formula. It is noticed that Casimir force can be switched from repulsive to attractive and vice versa via coherent control of the driving field. This switching feature provides new possibilities of using the repulsive Casimir force in the development of new emerging technologies, such as, micro-electro-mechanical and nano-electro-mechanical systems, i.e., MEMS and NEMS, respectively.
Improved dynamic compensation for accurate cutting force measurements in milling applications
NASA Astrophysics Data System (ADS)
Scippa, A.; Sallese, L.; Grossi, N.; Campatelli, G.
2015-03-01
Accurate cutting-force measurements appear to be the key information in most of the machining related studies as they are fundamental in understanding the cutting processes, optimizing the cutting operations and evaluating the presence of instabilities that could affect the effectiveness of cutting processes. A variety of specifically designed transducers are commercially available nowadays and many different approaches in measuring cutting forces are presented in literature. The available transducers, though, express some limitations since they are conditioned by the vibration of the surrounding system and by the transducer's natural frequency. These parameters can drastically affect the measurement accuracy in some cases; hence an effective and accurate tool is required to compensate those dynamically induced errors in cutting force measurements. This work is aimed at developing and testing a compensation technique based on Kalman filter estimator. Two different approaches named "band-fitting" and "parallel elaboration" methods, have been developed to extend applications of this compensation technique, especially for milling purpose. The compensation filter has been designed upon the experimentally identified system's dynamic and its accuracy and effectiveness has been evaluated by numerical and experimental tests. Finally its specific application in cutting force measurements compensation is described.
Current Status of Protein Force Fields for Molecular Dynamics
Lopes, Pedro E.M.; Guvench, Olgun
2015-01-01
Summary The current status of classical force fields for proteins is reviewed. These include additive force fields as well as the latest developments in the Drude and AMOEBA polarizable force fields. Parametrization strategies developed specifically for the Drude force field are described and compared with the additive CHARMM36 force field. Results from molecular simulations of proteins and small peptides are summarized to illustrate the performance of the Drude and AMOEBA force fields. PMID:25330958
Tailor-made force fields for crystal-structure prediction.
Neumann, Marcus A
2008-08-14
A general procedure is presented to derive a complete set of force-field parameters for flexible molecules in the crystalline state on a case-by-case basis. The force-field parameters are fitted to the electrostatic potential as well as to accurate energies and forces generated by means of a hybrid method that combines solid-state density functional theory (DFT) calculations with an empirical van der Waals correction. All DFT calculations are carried out with the VASP program. The mathematical structure of the force field, the generation of reference data, the choice of the figure of merit, the optimization algorithm, and the parameter-refinement strategy are discussed in detail. The approach is applied to cyclohexane-1,4-dione, a small flexible ring. The tailor-made force field obtained for cyclohexane-1,4-dione is used to search for low-energy crystal packings in all 230 space groups with one molecule per asymmetric unit, and the most stable crystal structures are reoptimized in a second step with the hybrid method. The experimental crystal structure is found as the most stable predicted crystal structure both with the tailor-made force field and the hybrid method. The same methodology has also been applied successfully to the four compounds of the fourth CCDC blind test on crystal-structure prediction. For the five aforementioned compounds, the root-mean-square deviations between lattice energies calculated with the tailor-made force fields and the hybrid method range from 0.024 to 0.053 kcal/mol per atom around an average value of 0.034 kcal/mol per atom. PMID:18642947
Li, Rui; Ye, Hongfei; Zhang, Weisheng; Ma, Guojun; Su, Yewang
2015-01-01
Spring constant calibration of the atomic force microscope (AFM) cantilever is of fundamental importance for quantifying the force between the AFM cantilever tip and the sample. The calibration within the framework of thin plate theory undoubtedly has a higher accuracy and broader scope than that within the well-established beam theory. However, thin plate theory-based accurate analytic determination of the constant has been perceived as an extremely difficult issue. In this paper, we implement the thin plate theory-based analytic modeling for the static behavior of rectangular AFM cantilevers, which reveals that the three-dimensional effect and Poisson effect play important roles in accurate determination of the spring constants. A quantitative scaling law is found that the normalized spring constant depends only on the Poisson’s ratio, normalized dimension and normalized load coordinate. Both the literature and our refined finite element model validate the present results. The developed model is expected to serve as the benchmark for accurate calibration of rectangular AFM cantilevers. PMID:26510769
Li, Rui; Ye, Hongfei; Zhang, Weisheng; Ma, Guojun; Su, Yewang
2015-01-01
Spring constant calibration of the atomic force microscope (AFM) cantilever is of fundamental importance for quantifying the force between the AFM cantilever tip and the sample. The calibration within the framework of thin plate theory undoubtedly has a higher accuracy and broader scope than that within the well-established beam theory. However, thin plate theory-based accurate analytic determination of the constant has been perceived as an extremely difficult issue. In this paper, we implement the thin plate theory-based analytic modeling for the static behavior of rectangular AFM cantilevers, which reveals that the three-dimensional effect and Poisson effect play important roles in accurate determination of the spring constants. A quantitative scaling law is found that the normalized spring constant depends only on the Poisson's ratio, normalized dimension and normalized load coordinate. Both the literature and our refined finite element model validate the present results. The developed model is expected to serve as the benchmark for accurate calibration of rectangular AFM cantilevers. PMID:26510769
NASA Astrophysics Data System (ADS)
Li, Rui; Ye, Hongfei; Zhang, Weisheng; Ma, Guojun; Su, Yewang
2015-10-01
Spring constant calibration of the atomic force microscope (AFM) cantilever is of fundamental importance for quantifying the force between the AFM cantilever tip and the sample. The calibration within the framework of thin plate theory undoubtedly has a higher accuracy and broader scope than that within the well-established beam theory. However, thin plate theory-based accurate analytic determination of the constant has been perceived as an extremely difficult issue. In this paper, we implement the thin plate theory-based analytic modeling for the static behavior of rectangular AFM cantilevers, which reveals that the three-dimensional effect and Poisson effect play important roles in accurate determination of the spring constants. A quantitative scaling law is found that the normalized spring constant depends only on the Poisson’s ratio, normalized dimension and normalized load coordinate. Both the literature and our refined finite element model validate the present results. The developed model is expected to serve as the benchmark for accurate calibration of rectangular AFM cantilevers.
NASA Technical Reports Server (NTRS)
Arnaiz, H. H.
1975-01-01
As part of a NASA program to evaluate current methods of predicting the performance of large, supersonic airplanes, the drag of the XB-70 airplane was measured accurately in flight at Mach numbers from 0.75 to 2.5. This paper describes the techniques used to determine engine net thrust and the drag forces charged to the propulsion system that were required for the in-flight drag measurements. The accuracy of the measurements and the application of the measurement techniques to aircraft with different propulsion systems are discussed. Examples of results obtained for the XB-70 airplane are presented.
Accurate formula for dissipative interaction in frequency modulation atomic force microscopy
Suzuki, Kazuhiro; Matsushige, Kazumi; Yamada, Hirofumi; Kobayashi, Kei; Labuda, Aleksander
2014-12-08
Much interest has recently focused on the viscosity of nano-confined liquids. Frequency modulation atomic force microscopy (FM-AFM) is a powerful technique that can detect variations in the conservative and dissipative forces between a nanometer-scale tip and a sample surface. We now present an accurate formula to convert the dissipation power of the cantilever measured during the experiment to damping of the tip-sample system. We demonstrated the conversion of the dissipation power versus tip-sample separation curve measured using a colloidal probe cantilever on a mica surface in water to the damping curve, which showed a good agreement with the theoretical curve. Moreover, we obtained the damping curve from the dissipation power curve measured on the hydration layers on the mica surface using a nanometer-scale tip, demonstrating that the formula allows us to quantitatively measure the viscosity of a nano-confined liquid using FM-AFM.
Accurate formula for dissipative interaction in frequency modulation atomic force microscopy
NASA Astrophysics Data System (ADS)
Suzuki, Kazuhiro; Kobayashi, Kei; Labuda, Aleksander; Matsushige, Kazumi; Yamada, Hirofumi
2014-12-01
Much interest has recently focused on the viscosity of nano-confined liquids. Frequency modulation atomic force microscopy (FM-AFM) is a powerful technique that can detect variations in the conservative and dissipative forces between a nanometer-scale tip and a sample surface. We now present an accurate formula to convert the dissipation power of the cantilever measured during the experiment to damping of the tip-sample system. We demonstrated the conversion of the dissipation power versus tip-sample separation curve measured using a colloidal probe cantilever on a mica surface in water to the damping curve, which showed a good agreement with the theoretical curve. Moreover, we obtained the damping curve from the dissipation power curve measured on the hydration layers on the mica surface using a nanometer-scale tip, demonstrating that the formula allows us to quantitatively measure the viscosity of a nano-confined liquid using FM-AFM.
Jaramillo-Botero, Andres; Naserifar, Saber; Goddard, William A
2014-04-01
First-principles-based force fields prepared from large quantum mechanical data sets are now the norm in predictive molecular dynamics simulations for complex chemical processes, as opposed to force fields fitted solely from phenomenological data. In principle, the former allow improved accuracy and transferability over a wider range of molecular compositions, interactions, and environmental conditions unexplored by experiments. That is, assuming they have been optimally prepared from a diverse training set. The trade-off has been force field engines that are functionally complex, with a large number of nonbonded and bonded analytical forms that give rise to rather large parameter search spaces. To address this problem, we have developed GARFfield (genetic algorithm-based reactive force field optimizer method), a hybrid multiobjective Pareto-optimal parameter development scheme based on genetic algorithms, hill-climbing routines and conjugate-gradient minimization. To demonstrate the capabilities of GARFfield we use it to develop two very different force fields: (1) the ReaxFF reactive force field for modeling the adiabatic reactive dynamics of silicon carbide growth from an methyltrichlorosilane precursor and (2) the SiC electron force field with effective core pseudopotentials for modeling nonadiabatic dynamic phenomena with highly excited electronic states. The flexible and open architecture of GARFfield enables efficient and fast parallel optimization of parameters from quantum mechanical data sets for demanding applications like ReaxFF, electronic fast forward (or electron force field), and others including atomistic reactive charge-optimized many-body interatomic potentials, Morse, and coarse-grain force fields. PMID:26580361
Planar Near-Field Phase Retrieval Using GPUs for Accurate THz Far-Field Prediction
NASA Astrophysics Data System (ADS)
Junkin, Gary
2013-04-01
With a view to using Phase Retrieval to accurately predict Terahertz antenna far-field from near-field intensity measurements, this paper reports on three fundamental advances that achieve very low algorithmic error penalties. The first is a new Gaussian beam analysis that provides accurate initial complex aperture estimates including defocus and astigmatic phase errors, based only on first and second moment calculations. The second is a powerful noise tolerant near-field Phase Retrieval algorithm that combines Anderson's Plane-to-Plane (PTP) with Fienup's Hybrid-Input-Output (HIO) and Successive Over-Relaxation (SOR) to achieve increased accuracy at reduced scan separations. The third advance employs teraflop Graphical Processing Units (GPUs) to achieve practically real time near-field phase retrieval and to obtain the optimum aperture constraint without any a priori information.
Current Status of the AMOEBA Polarizable Force Field
Ponder, Jay W.; Wu, Chuanjie; Ren, Pengyu; Pande, Vijay S.; Chodera, John D.; Schnieders, Michael J.; Haque, Imran; Mobley, David L.; Lambrecht, Daniel S.; DiStasio, Robert A.; Head-Gordon, Martin; Clark, Gary N. I.; Johnson, Margaret E.
2010-01-01
Molecular force fields have been approaching a generational transition over the past several years, moving away from well-established and well-tuned, but intrinsically limited, fixed point charge models towards more intricate and expensive polarizable models that should allow more accurate description of molecular properties. The recently introduced AMOEBA force field is a leading publicly available example of this next generation of theoretical model, but to date has only received relatively limited validation, which we address here. We show that the AMOEBA force field is in fact a significant improvement over fixed charge models for small molecule structural and thermodynamic observables in particular, although further fine-tuning is necessary to describe solvation free energies of drug-like small molecules, dynamical properties away from ambient conditions, and possible improvements in aromatic interactions. State of the art electronic structure calculations reveal generally very good agreement with AMOEBA for demanding problems such as relative conformational energies of the alanine tetrapeptide and isomers of water sulfate complexes. AMOEBA is shown to be especially successful on protein-ligand binding and computational X-ray crystallography where polarization and accurate electrostatics are critical. PMID:20136072
An improved generalized AMBER force field (GAFF) for urea.
Ozpinar, Gül Altinbaş; Peukert, Wolfgang; Clark, Timothy
2010-09-01
We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data. PMID:20162312
Methodology to set up accurate OPC model using optical CD metrology and atomic force microscopy
NASA Astrophysics Data System (ADS)
Shim, Yeon-Ah; Kang, Jaehyun; Lee, Sang-Uk; Kim, Jeahee; Kim, Keeho
2007-03-01
For the 90nm node and beyond, smaller Critical Dimension(CD) control budget is required and the ways to control good CD uniformity are needed. Moreover Optical Proximity Correction(OPC) for the sub-90nm node demands more accurate wafer CD data in order to improve accuracy of OPC model. Scanning Electron Microscope (SEM) is the typical method for measuring CD until ArF process. However SEM can give serious attack such as shrinkage of Photo Resist(PR) by burning of weak chemical structure of ArF PR due to high energy electron beam. In fact about 5nm CD narrowing occur when we measure CD by using CD-SEM in ArF photo process. Optical CD Metrology(OCD) and Atomic Force Microscopy(AFM) has been considered to the method for measuring CD without attack of organic materials. Also the OCD and AFM measurement system have the merits of speed, easiness and accurate data. For model-based OPC, the model is generated using CD data of test patterns transferred onto the wafer. In this study we discuss to generate accurate OPC model using OCD and AFM measurement system.
Accurate integral equation theory for the central force model of liquid water and ionic solutions
NASA Astrophysics Data System (ADS)
Ichiye, Toshiko; Haymet, A. D. J.
1988-10-01
The atom-atom pair correlation functions and thermodynamics of the central force model of water, introduced by Lemberg, Stillinger, and Rahman, have been calculated accurately by an integral equation method which incorporates two new developments. First, a rapid new scheme has been used to solve the Ornstein-Zernike equation. This scheme combines the renormalization methods of Allnatt, and Rossky and Friedman with an extension of the trigonometric basis-set solution of Labik and co-workers. Second, by adding approximate ``bridge'' functions to the hypernetted-chain (HNC) integral equation, we have obtained predictions for liquid water in which the hydrogen bond length and number are in good agreement with ``exact'' computer simulations of the same model force laws. In addition, for dilute ionic solutions, the ion-oxygen and ion-hydrogen coordination numbers display both the physically correct stoichiometry and good agreement with earlier simulations. These results represent a measurable improvement over both a previous HNC solution of the central force model and the ex-RISM integral equation solutions for the TIPS and other rigid molecule models of water.
Sultan - forced flow, high field test facility
Horvath, I.; Vecsey, G.; Weymuth, P.; Zellweger, J.
1981-09-01
Three European laboratories: CNEN (Frascati, I) ECN (Petten, NL) and SIN (Villigen, CH) decided to coordinate their development efforts and to install a common high field forced flow test facility at Villigen Switzerland. The test facility SULTAN (Supraleiter Testanlage) is presently under construction. As a first step, an 8T/1m bore solenoid with cryogenic periphery will be ready in 1981. The cryogenic system, data acquisition system and power supplies which are contributed by SIN are described. Experimental feasibilities, including cooling, and instrumentation are reviewed. Progress of components and facility construction is described. Planned extension of the background field up to 12T by insert coils is outlined. 5 refs.
CHARMM additive and polarizable force fields for biophysics and computer-aided drug design
Vanommeslaeghe, K.
2014-01-01
Background Molecular Mechanics (MM) is the method of choice for computational studies of biomolecular systems owing to its modest computational cost, which makes it possible to routinely perform molecular dynamics (MD) simulations on chemical systems of biophysical and biomedical relevance. Scope of Review As one of the main factors limiting the accuracy of MD results is the empirical force field used, the present paper offers a review of recent developments in the CHARMM additive force field, one of the most popular bimolecular force fields. Additionally, we present a detailed discussion of the CHARMM Drude polarizable force field, anticipating a growth in the importance and utilization of polarizable force fields in the near future. Throughout the discussion emphasis is placed on the force fields’ parametrization philosophy and methodology. Major Conclusions Recent improvements in the CHARMM additive force field are mostly related to newly found weaknesses in the previous generation of additive force fields. Beyond the additive approximation is the newly available CHARMM Drude polarizable force field, which allows for MD simulations of up to 1 microsecond on proteins, DNA, lipids and carbohydrates. General Significance Addressing the limitations ensures the reliability of the new CHARMM36 additive force field for the types of calculations that are presently coming into routine computational reach while the availability of the Drude polarizable force fields offers a model that is an inherently more accurate model of the underlying physical forces driving macromolecular structures and dynamics. PMID:25149274
The importance of accurate experimental data to marginal field development
Overa, S.J.; Lingelem, M.N.
1997-12-31
Since exploration started in the Norwegian North Sea in 1965 a total of 196 fields have been discovered. Less than one-third of these fields have been developed. The marginal fields can not be developed economically with current technology even though some of those fields have significant reserves. The total cost to develop one of those large installations is estimated to be 2--5 billion US dollars. Therefore new technology is needed to lower the designed and installed costs of each unit. The need for new physical property data is shown. The value of valid operating data from present units is also pointed out.
The variance needed to accurately describe jump height from vertical ground reaction force data.
Richter, Chris; McGuinness, Kevin; O'Connor, Noel E; Moran, Kieran
2014-12-01
In functional principal component analysis (fPCA) a threshold is chosen to define the number of retained principal components, which corresponds to the amount of preserved information. A variety of thresholds have been used in previous studies and the chosen threshold is often not evaluated. The aim of this study is to identify the optimal threshold that preserves the information needed to describe a jump height accurately utilizing vertical ground reaction force (vGRF) curves. To find an optimal threshold, a neural network was used to predict jump height from vGRF curve measures generated using different fPCA thresholds. The findings indicate that a threshold from 99% to 99.9% (6-11 principal components) is optimal for describing jump height, as these thresholds generated significantly lower jump height prediction errors than other thresholds. PMID:25010220
Magnetic field models of nine CP stars from "accurate" measurements
NASA Astrophysics Data System (ADS)
Glagolevskij, Yu. V.
2013-01-01
The dipole models of magnetic fields in nine CP stars are constructed based on the measurements of metal lines taken from the literature, and performed by the LSD method with an accuracy of 10-80 G. The model parameters are compared with the parameters obtained for the same stars from the hydrogen line measurements. For six out of nine stars the same type of structure was obtained. Some parameters, such as the field strength at the poles B p and the average surface magnetic field B s differ considerably in some stars due to differences in the amplitudes of phase dependences B e (Φ) and B s (Φ), obtained by different authors. It is noted that a significant increase in the measurement accuracy has little effect on the modelling of the large-scale structures of the field. By contrast, it is more important to construct the shape of the phase dependence based on a fairly large number of field measurements, evenly distributed by the rotation period phases. It is concluded that the Zeeman component measurement methods have a strong effect on the shape of the phase dependence, and that the measurements of the magnetic field based on the lines of hydrogen are more preferable for modelling the large-scale structures of the field.
Modeling of Non-Gravitational Forces for Precise and Accurate Orbit Determination
NASA Astrophysics Data System (ADS)
Hackel, Stefan; Gisinger, Christoph; Steigenberger, Peter; Balss, Ulrich; Montenbruck, Oliver; Eineder, Michael
2014-05-01
Remote sensing satellites support a broad range of scientific and commercial applications. The two radar imaging satellites TerraSAR-X and TanDEM-X provide spaceborne Synthetic Aperture Radar (SAR) and interferometric SAR data with a very high accuracy. The precise reconstruction of the satellite's trajectory is based on the Global Positioning System (GPS) measurements from a geodetic-grade dual-frequency Integrated Geodetic and Occultation Receiver (IGOR) onboard the spacecraft. The increasing demand for precise radar products relies on validation methods, which require precise and accurate orbit products. An analysis of the orbit quality by means of internal and external validation methods on long and short timescales shows systematics, which reflect deficits in the employed force models. Following the proper analysis of this deficits, possible solution strategies are highlighted in the presentation. The employed Reduced Dynamic Orbit Determination (RDOD) approach utilizes models for gravitational and non-gravitational forces. A detailed satellite macro model is introduced to describe the geometry and the optical surface properties of the satellite. Two major non-gravitational forces are the direct and the indirect Solar Radiation Pressure (SRP). The satellite TerraSAR-X flies on a dusk-dawn orbit with an altitude of approximately 510 km above ground. Due to this constellation, the Sun almost constantly illuminates the satellite, which causes strong across-track accelerations on the plane rectangular to the solar rays. The indirect effect of the solar radiation is called Earth Radiation Pressure (ERP). This force depends on the sunlight, which is reflected by the illuminated Earth surface (visible spectra) and the emission of the Earth body in the infrared spectra. Both components of ERP require Earth models to describe the optical properties of the Earth surface. Therefore, the influence of different Earth models on the orbit quality is assessed. The scope of
Nonlinear Force-Free Field Extrapolation of NOAA AR 0696
NASA Astrophysics Data System (ADS)
Thalmann, J. K.; Wiegelmann, T.
2007-12-01
We investigate the 3D coronal magnetic field structure of NOAA AR 0696 in the period of November 09-11, 2004, before and after an X2.5 flare (occurring around 02:13 UT on November 10, 2004). The coronal magnetic field dominates the structure of the solar corona and consequently plays a key role for the understanding of the initiation of flares. The most accurate presently available method to derive the coronal magnetic field is currently the nonlinear force-free field extrapolation from measurements of the photospheric magnetic field vector. These vector-magnetograms were processed from stokes I, Q, U, and V measurements of the Big Bear Solar Observatory and extrapolated into the corona with the nonlinear force-free optimization code developed by Wiegelmann (2004). We analyze the corresponding time series of coronal equilibria regarding topology changes of the 3D coronal magnetic field during the flare. Furthermore, quantities such as the temporal evolution of the magnetic energy and helicity are computed.
Force field dependence of riboswitch dynamics.
Hanke, Christian A; Gohlke, Holger
2015-01-01
Riboswitches are noncoding regulatory elements that control gene expression in response to the presence of metabolites, which bind to the aptamer domain. Metabolite binding appears to occur through a combination of conformational selection and induced fit mechanism. This demands to characterize the structural dynamics of the apo state of aptamer domains. In principle, molecular dynamics (MD) simulations can give insights at the atomistic level into the dynamics of the aptamer domain. However, it is unclear to what extent contemporary force fields can bias such insights. Here, we show that the Amber force field ff99 yields the best agreement with detailed experimental observations on differences in the structural dynamics of wild type and mutant aptamer domains of the guanine-sensing riboswitch (Gsw), including a pronounced influence of Mg2+. In contrast, applying ff99 with parmbsc0 and parmχOL modifications (denoted ff10) results in strongly damped motions and overly stable tertiary loop-loop interactions. These results are based on 58 MD simulations with an aggregate simulation time>11 μs, careful modeling of Mg2+ ions, and thorough statistical testing. Our results suggest that the moderate stabilization of the χ-anti region in ff10 can have an unwanted damping effect on functionally relevant structural dynamics of marginally stable RNA systems. This suggestion is supported by crystal structure analyses of Gsw aptamer domains that reveal χ torsions with high-anti values in the most mobile regions. We expect that future RNA force field development will benefit from considering marginally stable RNA systems and optimization toward good representations of dynamics in addition to structural characteristics. PMID:25726465
Rigorous force field optimization principles based on statistical distance minimization
Vlcek, Lukas; Chialvo, Ariel A.
2015-10-14
We use the concept of statistical distance to define a measure of distinguishability between a pair of statistical mechanical systems, i.e., a model and its target, and show that its minimization leads to general convergence of the model’s static measurable properties to those of the target. We exploit this feature to define a rigorous basis for the development of accurate and robust effective molecular force fields that are inherently compatible with coarse-grained experimental data. The new model optimization principles and their efficient implementation are illustrated through selected examples, whose outcome demonstrates the higher robustness and predictive accuracy of the approach compared to other currently used methods, such as force matching and relative entropy minimization. We also discuss relations between the newly developed principles and established thermodynamic concepts, which include the Gibbs-Bogoliubov inequality and the thermodynamic length.
Conformal field theory of critical Casimir forces
NASA Astrophysics Data System (ADS)
Emig, Thorsten; Bimonte, Giuseppe; Kardar, Mehran
2015-03-01
Thermal fluctuations of a critical system induce long-ranged Casimir forces between objects that couple to the underlying field. For two dimensional conformal field theories (CFT) we derive exact results for the Casimir interaction for a deformed strip and for two compact objects of arbitrary shape in terms of the free energy of a standard region (circular ring or flat strip) whose dimension is determined by the mutual capacitance of two conductors with the objects' shape; and a purely geometric energy that is proportional to conformal charge of the CFT, but otherwise super-universal in that it depends only on the shapes and is independent of boundary conditions and other details. The effect of inhomogenous boundary conditions is also discussed.
Bauer, Sebastian; Mathias, Gerald; Tavan, Paul
2014-03-14
We present a reaction field (RF) method which accurately solves the Poisson equation for proteins embedded in dielectric solvent continua at a computational effort comparable to that of an electrostatics calculation with polarizable molecular mechanics (MM) force fields. The method combines an approach originally suggested by Egwolf and Tavan [J. Chem. Phys. 118, 2039 (2003)] with concepts generalizing the Born solution [Z. Phys. 1, 45 (1920)] for a solvated ion. First, we derive an exact representation according to which the sources of the RF potential and energy are inducible atomic anti-polarization densities and atomic shielding charge distributions. Modeling these atomic densities by Gaussians leads to an approximate representation. Here, the strengths of the Gaussian shielding charge distributions are directly given in terms of the static partial charges as defined, e.g., by standard MM force fields for the various atom types, whereas the strengths of the Gaussian anti-polarization densities are calculated by a self-consistency iteration. The atomic volumes are also described by Gaussians. To account for covalently overlapping atoms, their effective volumes are calculated by another self-consistency procedure, which guarantees that the dielectric function ε(r) is close to one everywhere inside the protein. The Gaussian widths σ{sub i} of the atoms i are parameters of the RF approximation. The remarkable accuracy of the method is demonstrated by comparison with Kirkwood's analytical solution for a spherical protein [J. Chem. Phys. 2, 351 (1934)] and with computationally expensive grid-based numerical solutions for simple model systems in dielectric continua including a di-peptide (Ac-Ala-NHMe) as modeled by a standard MM force field. The latter example shows how weakly the RF conformational free energy landscape depends on the parameters σ{sub i}. A summarizing discussion highlights the achievements of the new theory and of its approximate solution
NASA Astrophysics Data System (ADS)
Bauer, Sebastian; Mathias, Gerald; Tavan, Paul
2014-03-01
We present a reaction field (RF) method which accurately solves the Poisson equation for proteins embedded in dielectric solvent continua at a computational effort comparable to that of an electrostatics calculation with polarizable molecular mechanics (MM) force fields. The method combines an approach originally suggested by Egwolf and Tavan [J. Chem. Phys. 118, 2039 (2003)] with concepts generalizing the Born solution [Z. Phys. 1, 45 (1920)] for a solvated ion. First, we derive an exact representation according to which the sources of the RF potential and energy are inducible atomic anti-polarization densities and atomic shielding charge distributions. Modeling these atomic densities by Gaussians leads to an approximate representation. Here, the strengths of the Gaussian shielding charge distributions are directly given in terms of the static partial charges as defined, e.g., by standard MM force fields for the various atom types, whereas the strengths of the Gaussian anti-polarization densities are calculated by a self-consistency iteration. The atomic volumes are also described by Gaussians. To account for covalently overlapping atoms, their effective volumes are calculated by another self-consistency procedure, which guarantees that the dielectric function ɛ(r) is close to one everywhere inside the protein. The Gaussian widths σi of the atoms i are parameters of the RF approximation. The remarkable accuracy of the method is demonstrated by comparison with Kirkwood's analytical solution for a spherical protein [J. Chem. Phys. 2, 351 (1934)] and with computationally expensive grid-based numerical solutions for simple model systems in dielectric continua including a di-peptide (Ac-Ala-NHMe) as modeled by a standard MM force field. The latter example shows how weakly the RF conformational free energy landscape depends on the parameters σi. A summarizing discussion highlights the achievements of the new theory and of its approximate solution particularly by
Bauer, Sebastian; Mathias, Gerald; Tavan, Paul
2014-03-14
We present a reaction field (RF) method which accurately solves the Poisson equation for proteins embedded in dielectric solvent continua at a computational effort comparable to that of an electrostatics calculation with polarizable molecular mechanics (MM) force fields. The method combines an approach originally suggested by Egwolf and Tavan [J. Chem. Phys. 118, 2039 (2003)] with concepts generalizing the Born solution [Z. Phys. 1, 45 (1920)] for a solvated ion. First, we derive an exact representation according to which the sources of the RF potential and energy are inducible atomic anti-polarization densities and atomic shielding charge distributions. Modeling these atomic densities by Gaussians leads to an approximate representation. Here, the strengths of the Gaussian shielding charge distributions are directly given in terms of the static partial charges as defined, e.g., by standard MM force fields for the various atom types, whereas the strengths of the Gaussian anti-polarization densities are calculated by a self-consistency iteration. The atomic volumes are also described by Gaussians. To account for covalently overlapping atoms, their effective volumes are calculated by another self-consistency procedure, which guarantees that the dielectric function ε(r) is close to one everywhere inside the protein. The Gaussian widths σ(i) of the atoms i are parameters of the RF approximation. The remarkable accuracy of the method is demonstrated by comparison with Kirkwood's analytical solution for a spherical protein [J. Chem. Phys. 2, 351 (1934)] and with computationally expensive grid-based numerical solutions for simple model systems in dielectric continua including a di-peptide (Ac-Ala-NHMe) as modeled by a standard MM force field. The latter example shows how weakly the RF conformational free energy landscape depends on the parameters σ(i). A summarizing discussion highlights the achievements of the new theory and of its approximate solution particularly by
Transferable force field for alcohols and polyalcohols.
Ferrando, Nicolas; Lachet, Véronique; Teuler, Jean-Marie; Boutin, Anne
2009-04-30
A new force field has been developed for alcohol and polyalcohol molecules. Based on the anisotropic united-atom force field AUA4 developed for hydrocarbons, it only introduces one new anisotropic united atom corresponding to the hydroxyl group OH. In the case of polyalcohols and complex molecules, the calculation of the intramolecular electrostatic energy is revisited. These interactions are calculated between charges belonging to the different local dipoles of the molecule, one dipole being defined as a group of consecutive charges globally neutral. This new method allows avoiding the use of empirical scaling parameters commonly introduced to calculate 1-4 electrostatic interactions. The transferability of the proposed potential is demonstrated through the simulation of a wide variety of alcohol families: primary alcohols (methanol, ethanol, propan-1-ol, hexan-1-ol, octan-1-ol), secondary alcohols (propan-2-ol), tertiary alcohols (2-methylpropan-2-ol), phenol, and diols (1,2-ethanediol, 1,3-propanediol, 1,5-pentanediol). Monte Carlo simulations carried out in the Gibbs ensemble lead to a good agreement between calculated and experimental data for the thermodynamic properties along the liquid/vapor saturation curve, for the critical point coordinates and for the liquid structure at room temperature. Additional simulations were performed on the methanol + n-butane system showing the capability of the proposed potential to reproduce the azeotropic behavior of such mixtures with a good agreement. PMID:19344171
Electromagnetic unification of matter and force fields
NASA Astrophysics Data System (ADS)
John, Sarah
2004-05-01
Special relativity and quantum mechanics are descriptive of electromagnetic propagation in waveguides, with mass analogous to the cutoff frequency of a waveguide mode [S.John, Bull.Am.Phys.Soc. vol.39,no.2,1254 (1994)]. It is further postulated herein that all spin 1/2 matter (necessarily massive) and spin 1 force fields have their origin in the electromagnetic fields E and B. This concept is not new. Majorana, among others have obtained electromagnetic representations of Dirac-like equations valid for the zero-mass case. Here, the spinor representation of the Maxwell equations, as given by Sallhofer, is extended to oscillatory fields with propagation constant m to obtain, in the absence of charge and current densities, the coupled equation (M. hatp + β E)ψ = 0 , where M = diag[ M σ, M^* σ ] , β = offdiag[I,I] , ψ ^ = i ^dag ( σ. B0 ( p), σ. E_0(p)), and M=m+ip, with the energy-mass relation given by E^2 = M M . Further, it is shown that the interaction term of QED is a direct consequence of including the sources and currents of Maxwell equations. Qualitative field patterns for spin 1/2 and spin 1 states, such as the electron, neutrino, magnetic monopole, quarks, photon, and massive gauge bosons are suggested.
Cosmic Ray Acceleration in Force Free Fields
NASA Astrophysics Data System (ADS)
Colgate, Stirling; Li, Hui; Kronberg, Philipp
2002-11-01
Galactic, extragalactic, and cluster magnetic fields are in apparent pressure equilibrium with the in-fall pressure of matter from the external medium, IGM, onto the Galaxies and clusters, and from the voids onto the galaxy sheets, (walls), implying fields of 5 , 0.5, & 20 μG respectively. Equipartition or minimum energy, implies β_CR=n_CRm_pc^2/(B^2/8π)˜= 1. The total energy in field and CRs is then ˜= 10^55 ergs Galactic and ˜= 4 ot 10^60 ergs per galaxy in the IGM and less within clusters, e.g., radio lobes, synchrotron "glow" in the IGM (Kronberg), and the UHECRs spectrum, Γ =-2.6. CRs escape from the Galaxy to the IGM, τ˜=10^7y, and similarly from the walls to the voids, ˜=10^8y, less than the GZK cut-off time provided B_galaxy>B_IGM>B_voids. The free energy of black hole formation, The Los Alamos model, is just sufficient. The lack of shocks at the boundaries of over pressured radio lobes and the need for high acceleration efficiency suggests eE_allel˜= eη_reconJ_allel, acceleration by reconnection of these force-free fields.
Weber's gravitational force as static weak field approximation
NASA Astrophysics Data System (ADS)
Tiandho, Yuant
2016-02-01
Weber's gravitational force (WGF) is one of gravitational model that can accommodate a non-static system because it depends not only on the distance but also on the velocity and the acceleration. Unlike Newton's law of gravitation, WGF can predict the anomalous of Mercury and gravitational bending of light near massive object very well. Then, some researchers use WGF as an alternative model of gravitation and propose a new mechanics theory namely the relational mechanics theory. However, currently we have known that the theory of general relativity which proposed by Einstein can explain gravity with very accurate. Through the static weak field approximation for the non-relativistic object, we also have known that the theory of general relativity will reduce to Newton's law of gravity. In this work, we expand the static weak field approximation that compatible with relativistic object and we obtain a force equation which correspond to WGF. Therefore, WGF is more precise than Newton's gravitational law. The static-weak gravitational field that we used is a solution of the Einstein's equation in the vacuum that satisfy the linear field approximation. The expression of WGF with ξ = 1 and satisfy the requirement of energy conservation are obtained after resolving the geodesic equation. By this result, we can conclude that WGF can be derived from the general relativity.
Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; MacKerell, A. D.
2010-01-01
The widely used CHARMM additive all-atom force field includes parameters for proteins, nucleic acids, lipids and carbohydrates. In the present paper an extension of the CHARMM force field to drug-like molecules is presented. The resulting CHARMM General Force Field (CGenFF) covers a wide range of chemical groups present in biomolecules and drug-like molecules, including a large number of heterocyclic scaffolds. The parametrization philosophy behind the force field focuses on quality at the expense of transferability, with the implementation concentrating on an extensible force field. Statistics related to the quality of the parametrization with a focus on experimental validation are presented. Additionally, the parametrization procedure, described fully in the present paper in the context of the model systems, pyrrolidine, and 3-phenoxymethylpyrrolidine will allow users to readily extend the force field to chemical groups that are not explicitly covered in the force field as well as add functional groups to and link together molecules already available in the force field. CGenFF thus makes it possible to perform “all-CHARMM” simulations on drug-target interactions thereby extending the utility of CHARMM force fields to medicinally relevant systems. PMID:19575467
Slattery, Ashley D; Blanch, Adam J; Quinton, Jamie S; Gibson, Christopher T
2013-08-01
calibrate the cantilever spring constant using the thermal noise method, allowing complete force calibration to be accurately performed without tip-sample contact. PMID:23685172
Nonlinear gravitational self-force: Field outside a small body
NASA Astrophysics Data System (ADS)
Pound, Adam
2012-10-01
A small extended body moving through an external spacetime gαβ creates a metric perturbation hαβ, which forces the body away from geodesic motion in gαβ. The foundations of this effect, called the gravitational self-force, are now well established, but concrete results have mostly been limited to linear order. Accurately modeling the dynamics of compact binaries requires proceeding to nonlinear orders. To that end, I show how to obtain the metric perturbation outside the body at all orders in a class of generalized wave gauges. In a small buffer region surrounding the body, the form of the perturbation can be found analytically as an expansion for small distances r from a representative worldline. Given only a specification of the body’s multipole moments, the field obtained in the buffer region suffices to find the metric everywhere outside the body via a numerical puncture scheme. Following this procedure at first and second order, I calculate the field in the buffer region around an arbitrarily structured compact body at sufficiently high order in r to numerically implement a second-order puncture scheme, including effects of the body’s spin. I also define nth-order (local) generalizations of the Detweiler-Whiting singular and regular fields and show that in a certain sense, the body can be viewed as a skeleton of multipole moments.
The Energetics of Motivated Cognition: A Force-Field Analysis
ERIC Educational Resources Information Center
Kruglanski, Arie W.; Belanger, Jocelyn J.; Chen, Xiaoyan; Kopetz, Catalina; Pierro, Antonio; Mannetti, Lucia
2012-01-01
A force-field theory of motivated cognition is presented and applied to a broad variety of phenomena in social judgment and self-regulation. Purposeful cognitive activity is assumed to be propelled by a "driving force" and opposed by a "restraining force". "Potential" driving force represents the maximal amount of energy an individual is prepared…
Motor imagery facilitates force field learning.
Anwar, Muhammad Nabeel; Tomi, Naoki; Ito, Koji
2011-06-13
Humans have the ability to produce an internal reproduction of a specific motor action without any overt motor output. Recent findings show that the processes underlying motor imagery are similar to those active during motor execution and both share common neural substrates. This suggests that the imagery of motor movements might play an important role in acquiring new motor skills. In this study we used haptic robot in conjunction with motor imagery technique to improve learning in a robot-based adaptation task. Two groups of subjects performed reaching movements with or without motor imagery in a velocity-dependent and position-dependent mixed force field. The groups performed movements with motor imagery produced higher after effects and decreased muscle co-contraction with respect to no-motor imagery group. These results showed a positive influence of motor imagery on acquiring new motor skill and suggest that motor learning can be facilitated by mental practice and could be used to increase the rate of adaptation. PMID:21555118
The Introduction of Fields in Relation to Force
ERIC Educational Resources Information Center
Brunt, Marjorie; Brunt, Geoff
2012-01-01
The introduction of force at age 14-16 years is considered, starting with elementary student experiments using magnetic force fields. The meaningless use of terms such as "action" and "reaction", or "agent" and "receiver" is discussed. (Contains 6 figures.)
Development of force field parameters for molecular simulation of polylactide
McAliley, James H.; Bruce, David A.
2011-01-01
Polylactide is a biodegradable polymer that is widely used for biomedical applications, and it is a replacement for some petroleum based polymers in applications that range from packaging to carpeting. Efforts to characterize and further enhance polylactide based systems using molecular simulations have to this point been hindered by the lack of accurate atomistic models for the polymer. Thus, we present force field parameters specifically suited for molecular modeling of PLA. The model, which we refer to as PLAFF3, is based on a combination of the OPLS and CHARMM force fields, with modifications to bonded and nonbonded parameters. Dihedral angle parameters were adjusted to reproduce DFT data using newly developed CMAP dihedral cross terms, and the model was further adjusted to reproduce experimentally resolved crystal structure conformations, melt density, volume expansivity, and the glass transition temperature of PLA. We recommend the use of PLAFF3 in modeling PLA in its crystalline or amorphous states and have provided the necessary input files required for the publicly available molecular dynamics code GROMACS. PMID:22180734
NASA Astrophysics Data System (ADS)
Luo, Ye; Sorella, Sandro
2014-03-01
We introduce a general and efficient method for the calculation of vibrational frequencies of electronic systems, ranging from molecules to solids. By performing damped molecular dynamics with ab initio forces, we show that quantum vibrational frequencies can be evaluated by diagonalizing the time averaged position-position or force-force correlation matrices, although the ionic motion is treated on the classical level within the Born-Oppenheimer approximation. The novelty of our approach is to evaluate atomic forces with QMC by means of a highly accurate and correlated variational wave function which is optimized simultaneously during the dynamics. QMC is an accurate and promising many-body technique for electronic structure calculation thanks to massively parallel computers. However, since infinite statistics is not feasible, property evaluation may be affected by large noise that is difficult to harness. Our approach controls the QMC stochastic bias systematically and gives very accurate results with moderate computational effort, namely even with noisy forces. We prove the accuracy and efficiency of our method on the water monomer[A. Zen et al., JCTC 9 (2013) 4332] and dimer. We are currently working on the challenging problem of simulating liquid water at ambient conditions.
Dislocation core fields and forces in FCC metals
Henager, Charles H.; Hoagland, Richard G.
2004-04-01
Atomistic models were used to obtain dislocation core fields for edge, screw, and mixed dislocations in Al and Cu using EAM. Core fields are analyzed using a line force dipole representation, with dilatant and dipole terms. The core field contribution to the force between dislocations is shown to be significant for interactions within 50b.
Accurate spring constant calibration for very stiff atomic force microscopy cantilevers
Grutzik, Scott J.; Zehnder, Alan T.; Gates, Richard S.; Gerbig, Yvonne B.; Smith, Douglas T.; Cook, Robert F.
2013-11-15
There are many atomic force microscopy (AFM) applications that rely on quantifying the force between the AFM cantilever tip and the sample. The AFM does not explicitly measure force, however, so in such cases knowledge of the cantilever stiffness is required. In most cases, the forces of interest are very small, thus compliant cantilevers are used. A number of methods have been developed that are well suited to measuring low stiffness values. However, in some cases a cantilever with much greater stiffness is required. Thus, a direct, traceable method for calibrating very stiff (approximately 200 N/m) cantilevers is presented here. The method uses an instrumented and calibrated nanoindenter to determine the stiffness of a reference cantilever. This reference cantilever is then used to measure the stiffness of a number of AFM test cantilevers. This method is shown to have much smaller uncertainty than previously proposed methods. An example application to fracture testing of nanoscale silicon beam specimens is included.
NASA Astrophysics Data System (ADS)
Zhang, Xiang; Vu-Quoc, Loc
2007-07-01
We present in this paper the displacement-driven version of a tangential force-displacement (TFD) model that accounts for both elastic and plastic deformations together with interfacial friction occurring in collisions of spherical particles. This elasto-plastic frictional TFD model, with its force-driven version presented in [L. Vu-Quoc, L. Lesburg, X. Zhang. An accurate tangential force-displacement model for granular-flow simulations: contacting spheres with plastic deformation, force-driven formulation, Journal of Computational Physics 196(1) (2004) 298-326], is consistent with the elasto-plastic frictional normal force-displacement (NFD) model presented in [L. Vu-Quoc, X. Zhang. An elasto-plastic contact force-displacement model in the normal direction: displacement-driven version, Proceedings of the Royal Society of London, Series A 455 (1991) 4013-4044]. Both the NFD model and the present TFD model are based on the concept of additive decomposition of the radius of contact area into an elastic part and a plastic part. The effect of permanent indentation after impact is represented by a correction to the radius of curvature. The effect of material softening due to plastic flow is represented by a correction to the elastic moduli. The proposed TFD model is accurate, and is validated against nonlinear finite element analyses involving plastic flows in both the loading and unloading conditions. The proposed consistent displacement-driven, elasto-plastic NFD and TFD models are designed for implementation in computer codes using the discrete-element method (DEM) for granular-flow simulations. The model is shown to be accurate and is validated against nonlinear elasto-plastic finite-element analysis.
Towards a force field based on density fitting
Piquemal, Jean-Philip; Cisneros, G. Andrés; Reinhardt, Peter; Gresh, Nohad; Darden, Thomas A.
2007-01-01
Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n=16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed. PMID:16542062
Yin, Jian; Fenley, Andrew T.; Henriksen, Niel M.; Gilson, Michael K.
2015-01-01
Improving the capability of atomistic computer models to predict the thermodynamics of noncovalent binding is critical for successful structure-based drug design, and the accuracy of such calculations remains limited by non-optimal force field parameters. Ideally, one would incorporate protein-ligand affinity data into force field parametrization, but this would be inefficient and costly. We now demonstrate that sensitivity analysis can be used to efficiently tune Lennard-Jones parameters of aqueous host-guest systems for increasingly accurate calculations of binding enthalpy. These results highlight the promise of a comprehensive use of calorimetric host-guest binding data, along with existing validation data sets, to improve force field parameters for the simulation of noncovalent binding, with the ultimate goal of making protein-ligand modeling more accurate and hence speeding drug discovery. PMID:26181208
Lorentz Body Force Induced by Traveling Magnetic Fields
NASA Technical Reports Server (NTRS)
Volz, M. P.; Mazuruk, K.
2003-01-01
The Lorentz force induced by a traveling magnetic field (TMF) in a cylindrical container has been calculated. The force can be used to control flow in dectrically conducting melts and the direction of the magnetic field and resulting flow can be reversed. A TMF can be used to partially cancel flow driven by buoyancy. The penetration of the field into the cylinder decreases as the frequency increases, and there exists an optimal value of frequency for which the resulting force is a maximum. Expressions for the Lorentz force in the limiting cases of low frequency and infinite cylinder are also given and compared to the numerical calculations.
Theoretical study on the HIV-1 integrase-5CITEP complex based on polarized force fields
NASA Astrophysics Data System (ADS)
Wei, Caiyi; Mei, Ye; Zhang, Dawei
2010-07-01
Molecular dynamics studies of 5CITEP binding with HIV-1 integrase (IN) are presented using both polarized and nonpolarized force fields. When nonpolarized force field is used, the ligand drifts away from the original binding site. However, this depressing behavior can be curbed by introducing electronic polarization effect into the force field that stabilizes the protein structure and keeps the ligand in the binding pocket. Moreover, simulation under polarized force field gives a binding energy of -4.85 kcal/mol which is in excellent agreement with the experimental Δ G of -4.38 kcal/mol. The results demonstrate the importance of intra-protein electronic polarization in stabilizing the binding complex of IN-5CITEP and accurately predicting the binding energy.
Assessing the Current State of Amber Force Field Modifications for DNA.
Galindo-Murillo, Rodrigo; Robertson, James C; Zgarbová, Marie; Šponer, Jiří; Otyepka, Michal; Jurečka, Petr; Cheatham, Thomas E
2016-08-01
The utility of molecular dynamics (MD) simulations to model biomolecular structure, dynamics, and interactions has witnessed enormous advances in recent years due to the availability of optimized MD software and access to significant computational power, including GPU multicore computing engines and other specialized hardware. This has led researchers to routinely extend conformational sampling times to the microsecond level and beyond. The extended sampling time has allowed the community not only to converge conformational ensembles through complete sampling but also to discover deficiencies and overcome problems with the force fields. Accuracy of the force fields is a key component, along with sampling, toward being able to generate accurate and stable structures of biopolymers. The Amber force field for nucleic acids has been used extensively since the 1990s, and multiple artifacts have been discovered, corrected, and reassessed by different research groups. We present a direct comparison of two of the most recent and state-of-the-art Amber force field modifications, bsc1 and OL15, that focus on accurate modeling of double-stranded DNA. After extensive MD simulations with five test cases and two different water models, we conclude that both modifications are a remarkable improvement over the previous bsc0 force field. Both force field modifications show better agreement when compared to experimental structures. To ensure convergence, the Drew-Dickerson dodecamer (DDD) system was simulated using 100 independent MD simulations, each extended to at least 10 μs, and the independent MD simulations were concatenated into a single 1 ms long trajectory for each combination of force field and water model. This is significantly beyond the time scale needed to converge the conformational ensemble of the internal portions of a DNA helix absent internal base pair opening. Considering all of the simulations discussed in the current work, the MD simulations performed to
Assessing the Current State of Amber Force Field Modifications for DNA
2016-01-01
The utility of molecular dynamics (MD) simulations to model biomolecular structure, dynamics, and interactions has witnessed enormous advances in recent years due to the availability of optimized MD software and access to significant computational power, including GPU multicore computing engines and other specialized hardware. This has led researchers to routinely extend conformational sampling times to the microsecond level and beyond. The extended sampling time has allowed the community not only to converge conformational ensembles through complete sampling but also to discover deficiencies and overcome problems with the force fields. Accuracy of the force fields is a key component, along with sampling, toward being able to generate accurate and stable structures of biopolymers. The Amber force field for nucleic acids has been used extensively since the 1990s, and multiple artifacts have been discovered, corrected, and reassessed by different research groups. We present a direct comparison of two of the most recent and state-of-the-art Amber force field modifications, bsc1 and OL15, that focus on accurate modeling of double-stranded DNA. After extensive MD simulations with five test cases and two different water models, we conclude that both modifications are a remarkable improvement over the previous bsc0 force field. Both force field modifications show better agreement when compared to experimental structures. To ensure convergence, the Drew–Dickerson dodecamer (DDD) system was simulated using 100 independent MD simulations, each extended to at least 10 μs, and the independent MD simulations were concatenated into a single 1 ms long trajectory for each combination of force field and water model. This is significantly beyond the time scale needed to converge the conformational ensemble of the internal portions of a DNA helix absent internal base pair opening. Considering all of the simulations discussed in the current work, the MD simulations performed to
Catch trials in force field learning influence adaptation and consolidation of human motor memory
Stockinger, Christian; Focke, Anne; Stein, Thorsten
2014-01-01
Force field studies are a common tool to investigate motor adaptation and consolidation. Thereby, subjects usually adapt their reaching movements to force field perturbations induced by a robotic device. In this context, so-called catch trials, in which the disturbing forces are randomly turned off, are commonly used to detect after-effects of motor adaptation. However, catch trials also produce sudden large motor errors that might influence the motor adaptation and the consolidation process. Yet, the detailed influence of catch trials is far from clear. Thus, the aim of this study was to investigate the influence of catch trials on motor adaptation and consolidation in force field experiments. Therefore, 105 subjects adapted their reaching movements to robot-generated force fields. The test groups adapted their reaching movements to a force field A followed by learning a second interfering force field B before retest of A (ABA). The control groups were not exposed to force field B (AA). To examine the influence of diverse catch trial ratios, subjects received catch trials during force field adaptation with a probability of either 0, 10, 20, 30, or 40%, depending on the group. First, the results on motor adaptation revealed significant differences between the diverse catch trial ratio groups. With increasing amount of catch trials, the subjects' motor performance decreased and subjects' ability to accurately predict the force field—and therefore internal model formation—was impaired. Second, our results revealed that adapting with catch trials can influence the following consolidation process as indicated by a partial reduction to interference. Here, the optimal catch trial ratio was 30%. However, detection of consolidation seems to be biased by the applied measure of performance. PMID:24795598
NASA Astrophysics Data System (ADS)
Lonardoni, D.; Pederiva, F.; Gandolfi, S.
2014-01-01
Background: An accurate assessment of the hyperon-nucleon interaction is of great interest in view of recent observations of very massive neutron stars. The challenge is to build a realistic interaction that can be used over a wide range of masses and in infinite matter starting from the available experimental data on the binding energy of light hypernuclei. To this end, accurate calculations of the hyperon binding energy in a hypernucleus are necessary. Purpose: We present a quantum Monte Carlo study of Λ and ΛΛ hypernuclei up to A =91. We investigate the contribution of two- and three-body Λ-nucleon forces to the Λ binding energy. Method: Ground state energies are computed solving the Schrödinger equation for nonrelativistic baryons by means of the auxiliary field diffusion Monte Carlo algorithm extended to the hypernuclear sector. Results: We show that a simple adjustment of the parameters of the ΛNN three-body force yields a very good agreement with available experimental data over a wide range of hypernuclear masses. In some cases no experiments have been performed yet, and we give new predictions. Conclusions: The newly fitted ΛNN force properly describes the physics of medium-heavy Λ hypernuclei, correctly reproducing the saturation property of the hyperon separation energy.
Ponderomotive Force in the Presence of Electric Fields
NASA Technical Reports Server (NTRS)
Khazanov, G. V.; Krivorutsky, E. N.
2013-01-01
This paper presents averaged equations of particle motion in an electromagnetic wave of arbitrary frequency with its wave vector directed along the ambient magnetic field. The particle is also subjected to an E cross B drift and a background electric field slowly changing in space and acting along the magnetic field line. The fields, wave amplitude, and the wave vector depend on the coordinate along the magnetic field line. The derivations of the ponderomotive forces are done by assuming that the drift velocity in the ambient magnetic field is comparable to the particle velocity. Such a scenario leads to new ponderomotive forces, dependent on the wave magnetic field intensity, and, as a result, to the additional energy exchange between the wave and the plasma particles. It is found that the parallel electric field can lead to the change of the particle-wave energy exchange rate comparable to that produced by the previously discussed ponderomotive forces.
Novel concepts in near-field optics: from magnetic near-field to optical forces
NASA Astrophysics Data System (ADS)
Yang, Honghua
Driven by the progress in nanotechnology, imaging and spectroscopy tools with nanometer spatial resolution are needed for in situ material characterizations. Near-field optics provides a unique way to selectively excite and detect elementary electronic and vibrational interactions at the nanometer scale, through interactions of light with matter in the near-field region. This dissertation discusses the development and applications of near-field optical imaging techniques, including plasmonic material characterization, optical spectral nano-imaging and magnetic field detection using scattering-type scanning near-field optical microscopy (s-SNOM), and exploring new modalities of optical spectroscopy based on optical gradient force detection. Firstly, the optical dielectric functions of one of the most common plasmonic materials---silver is measured with ellipsometry, and analyzed with the Drude model over a broad spectral range from visible to mid-infrared. This work was motivated by the conflicting results of previous measurements, and the need for accurate values for a wide range of applications of silver in plasmonics, optical antennas, and metamaterials. This measurement provides a reference for dielectric functions of silver used in metamaterials, plasmonics, and nanophotonics. Secondly, I implemented an infrared s-SNOM instrument for spectroscopic nano-imaging at both room temperature and low temperature. As one of the first cryogenic s-SNOM instruments, the novel design concept and key specifications are discussed. Initial low-temperature and high-temperature performances of the instrument are examined by imaging of optical conductivity of vanadium oxides (VO2 and V2O 3) across their phase transitions. The spectroscopic imaging capability is demonstrated on chemical vibrational resonances of Poly(methyl methacrylate) (PMMA) and other samples. The third part of this dissertation explores imaging of optical magnetic fields. As a proof-of-principle, the magnetic
reaxFF Reactive Force Field for Disulfide Mechanochemistry, Fitted to Multireference ab Initio Data.
Müller, Julian; Hartke, Bernd
2016-08-01
Mechanochemistry, in particular in the form of single-molecule atomic force microscopy experiments, is difficult to model theoretically, for two reasons: Covalent bond breaking is not captured accurately by single-determinant, single-reference quantum chemistry methods, and experimental times of milliseconds or longer are hard to simulate with any approach. Reactive force fields have the potential to alleviate both problems, as demonstrated in this work: Using nondeterministic global parameter optimization by evolutionary algorithms, we have fitted a reaxFF force field to high-level multireference ab initio data for disulfides. The resulting force field can be used to reliably model large, multifunctional mechanochemistry units with disulfide bonds as designed breaking points. Explorative calculations show that a significant part of the time scale gap between AFM experiments and dynamical simulations can be bridged with this approach. PMID:27415976
Comparison of Cellulose Ib Simulations with Three Carbohydrate Force Fields
Matthews, J. F.; Beckham, G. T.; Bergenstrahle, M.; Brady, J. W.; Himmel, M. E.; Crowley, M. F.
2012-02-14
Molecular dynamics simulations of cellulose have recently become more prevalent due to increased interest in renewable energy applications, and many atomistic and coarse-grained force fields exist that can be applied to cellulose. However, to date no systematic comparison between carbohydrate force fields has been conducted for this important system. To that end, we present a molecular dynamics simulation study of hydrated, 36-chain cellulose I{beta} microfibrils at room temperature with three carbohydrate force fields (CHARMM35, GLYCAM06, and Gromos 45a4) up to the near-microsecond time scale. Our results indicate that each of these simulated microfibrils diverge from the cellulose I{beta} crystal structure to varying degrees under the conditions tested. The CHARMM35 and GLYCAM06 force fields eventually result in structures similar to those observed at 500 K with the same force fields, which are consistent with the experimentally observed high-temperature behavior of cellulose I. The third force field, Gromos 45a4, produces behavior significantly different from experiment, from the other two force fields, and from previously reported simulations with this force field using shorter simulation times and constrained periodic boundary conditions. For the GLYCAM06 force field, initial hydrogen-bond conformations and choice of electrostatic scaling factors significantly affect the rate of structural divergence. Our results suggest dramatically different time scales for convergence of properties of interest, which is important in the design of computational studies and comparisons to experimental data. This study highlights that further experimental and theoretical work is required to understand the structure of small diameter cellulose microfibrils typical of plant cellulose.
PARMBSC1: A REFINED FORCE-FIELD FOR DNA SIMULATIONS
Ivani, Ivan; Dans, Pablo D.; Noy, Agnes; Pérez, Alberto; Faustino, Ignacio; Hospital, Adam; Walther, Jürgen; Andrio, Pau; Goñi, Ramon; Balaceanu, Alexandra; Portella, Guillem; Battistini, Federica; Gelpí, Josep Lluis; González, Carlos; Vendruscolo, Michele; Laughton, Charles A.; Harris, Sarah A.; Case, David A.; Orozco, Modesto
2015-01-01
We present parmbsc1, a new force-field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (~140 μs) covering most of the DNA structural space. Parmbsc1 provides high quality results in diverse systems, solving problems of previous force-fields. Parmbsc1 aims to be a reference force-field for the study of DNA in the next decade. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/. PMID:26569599
Approximate photochemical dynamics of azobenzene with reactive force fields
Li, Yan; Hartke, Bernd
2013-12-14
We have fitted reactive force fields of the ReaxFF type to the ground and first excited electronic states of azobenzene, using global parameter optimization by genetic algorithms. Upon coupling with a simple energy-gap transition probability model, this setup allows for completely force-field-based simulations of photochemical cis→trans- and trans→cis-isomerizations of azobenzene, with qualitatively acceptable quantum yields. This paves the way towards large-scale dynamics simulations of molecular machines, including bond breaking and formation (via the reactive force field) as well as photochemical engines (presented in this work)
Accurate near-field calculation in the rigorous coupled-wave analysis method
NASA Astrophysics Data System (ADS)
Weismann, Martin; Gallagher, Dominic F. G.; Panoiu, Nicolae C.
2015-12-01
The rigorous coupled-wave analysis (RCWA) is one of the most successful and widely used methods for modeling periodic optical structures. It yields fast convergence of the electromagnetic far-field and has been adapted to model various optical devices and wave configurations. In this article, we investigate the accuracy with which the electromagnetic near-field can be calculated by using RCWA and explain the observed slow convergence and numerical artifacts from which it suffers, namely unphysical oscillations at material boundaries due to the Gibbs phenomenon. In order to alleviate these shortcomings, we also introduce a mathematical formulation for accurate near-field calculation in RCWA, for one- and two-dimensional straight and slanted diffraction gratings. This accurate near-field computational approach is tested and evaluated for several representative test-structures and configurations in order to illustrate the advantages provided by the proposed modified formulation of the RCWA.
The interoperability force in the ERP field
NASA Astrophysics Data System (ADS)
Boza, Andrés; Cuenca, Llanos; Poler, Raúl; Michaelides, Zenon
2015-04-01
Enterprise resource planning (ERP) systems participate in interoperability projects and this participation sometimes leads to new proposals for the ERP field. The aim of this paper is to identify the role that interoperability plays in the evolution of ERP systems. To go about this, ERP systems have been first identified within interoperability frameworks. Second, the initiatives in the ERP field driven by interoperability requirements have been identified from two perspectives: technological and business. The ERP field is evolving from classical ERP as information system integrators to a new generation of fully interoperable ERP. Interoperability is changing the way of running business, and ERP systems are changing to adapt to the current stream of interoperability.
Prediction of Mechanical Properties of Polymers With Various Force Fields
NASA Technical Reports Server (NTRS)
Odegard, Gregory M.; Clancy, Thomas C.; Gates, Thomas S.
2005-01-01
The effect of force field type on the predicted elastic properties of a polyimide is examined using a multiscale modeling technique. Molecular Dynamics simulations are used to predict the atomic structure and elastic properties of the polymer by subjecting a representative volume element of the material to bulk and shear finite deformations. The elastic properties of the polyimide are determined using three force fields: AMBER, OPLS-AA, and MM3. The predicted values of Young s modulus and shear modulus of the polyimide are compared with experimental values. The results indicate that the mechanical properties of the polyimide predicted with the OPLS-AA force field most closely matched those from experiment. The results also indicate that while the complexity of the force field does not have a significant effect on the accuracy of predicted properties, small differences in the force constants and the functional form of individual terms in the force fields determine the accuracy of the force field in predicting the elastic properties of the polyimide.
Vibrational spectrum and force field of dimethyldimethoxysilane
Tenisheva, T.F.; Lazarev, A.N.
1986-01-01
Experimental data is presented on the spectra of (CH/sub 3/)/sub 2/Si(OCH/sub 3/)/sub 2/ (I), (CH/sub 3/)/sub 2/Si(OCD/sub 3/)/sub 2/ (II), and (CD/sub 3/)/sub 2/Si(OCH/sub 3/)/sub 2/ (III). The results of the determination of the force constants on the basis of the optimization of the solution of the inverse mechanical problem of the theory of molecular vibrations with consideration of all the internal degrees of freedom with the exception of the coordinates corresponding to internal rotations are discussed. Raman spectra of I, II, and III in the liquid phase is shown and the IR spectra of amorphous films of I, II, and III are illustrated.
A transferable force field for CdS-CdSe-PbS-PbSe solid systems
Fan, Zhaochuan; Vlugt, Thijs J. H.; Koster, Rik S.; Fang, Changming; Huis, Marijn A. van; Wang, Shuaiwei; Yalcin, Anil O.; Tichelaar, Frans D.; Zandbergen, Henny W.
2014-12-28
A transferable force field for the PbSe-CdSe solid system using the partially charged rigid ion model has been successfully developed and was used to study the cation exchange in PbSe-CdSe heteronanocrystals [A. O. Yalcin et al., “Atomic resolution monitoring of cation exchange in CdSe-PbSe heteronanocrystals during epitaxial solid-solid-vapor growth,” Nano Lett. 14, 3661–3667 (2014)]. In this work, we extend this force field by including another two important binary semiconductors, PbS and CdS, and provide detailed information on the validation of this force field. The parameterization combines Bader charge analysis, empirical fitting, and ab initio energy surface fitting. When compared with experimental data and density functional theory calculations, it is shown that a wide range of physical properties of bulk PbS, PbSe, CdS, CdSe, and their mixed phases can be accurately reproduced using this force field. The choice of functional forms and parameterization strategy is demonstrated to be rational and effective. This transferable force field can be used in various studies on II-VI and IV-VI semiconductor materials consisting of CdS, CdSe, PbS, and PbSe. Here, we demonstrate the applicability of the force field model by molecular dynamics simulations whereby transformations are initiated by cation exchange.
Molecular dynamics simulations of methane hydrate using polarizable force fields
Jiang, H.N.; Jordan, K.D.; Taylor, C.E.
2007-06-14
Molecular dynamics simulations of methane hydrate have been carried out using the polarizable AMOEBA and COS/G2 force fields. Properties calculated include the temperature dependence of the lattice constant, the OC and OO radial distribution functions, and the vibrational spectra. Both the AMOEBA and COS/G2 force fields are found to successfully account for the available experimental data, with overall somewhat better agreement with experiment being found for the AMOEBA model. Comparison is made with previous results obtained using TIP4P and SPC/E effective two-body force fields and the polarizable TIP4P-FQ force field, which allows for in-plane polarization only. Significant differences are found between the properties calculated using the TIP4P-FQ model and those obtained using the other models, indicating an inadequacy of restricting explicit polarization to in-plane onl
Wallace, W C; Ghafur, O; Khurmi, C; Sainadh U, Satya; Calvert, J E; Laban, D E; Pullen, M G; Bartschat, K; Grum-Grzhimailo, A N; Wells, D; Quiney, H M; Tong, X M; Litvinyuk, I V; Sang, R T; Kielpinski, D
2016-07-29
Ionization of atoms and molecules in strong laser fields is a fundamental process in many fields of research, especially in the emerging field of attosecond science. So far, demonstrably accurate data have only been acquired for atomic hydrogen (H), a species that is accessible to few investigators. Here, we present measurements of the ionization yield for argon, krypton, and xenon with percent-level accuracy, calibrated using H, in a laser regime widely used in attosecond science. We derive a transferable calibration standard for laser peak intensity, accurate to 1.3%, that is based on a simple reference curve. In addition, our measurements provide a much needed benchmark for testing models of ionization in noble-gas atoms, such as the widely employed single-active electron approximation. PMID:27517769
NASA Astrophysics Data System (ADS)
Wallace, W. C.; Ghafur, O.; Khurmi, C.; Sainadh U, Satya; Calvert, J. E.; Laban, D. E.; Pullen, M. G.; Bartschat, K.; Grum-Grzhimailo, A. N.; Wells, D.; Quiney, H. M.; Tong, X. M.; Litvinyuk, I. V.; Sang, R. T.; Kielpinski, D.
2016-07-01
Ionization of atoms and molecules in strong laser fields is a fundamental process in many fields of research, especially in the emerging field of attosecond science. So far, demonstrably accurate data have only been acquired for atomic hydrogen (H), a species that is accessible to few investigators. Here, we present measurements of the ionization yield for argon, krypton, and xenon with percent-level accuracy, calibrated using H, in a laser regime widely used in attosecond science. We derive a transferable calibration standard for laser peak intensity, accurate to 1.3%, that is based on a simple reference curve. In addition, our measurements provide a much needed benchmark for testing models of ionization in noble-gas atoms, such as the widely employed single-active electron approximation.
Mitigated-force carriage for high magnetic field environments
Ludtka, Gerard M.; Ludtka, Gail M.; Wilgen, John B.; Murphy, Bart L.
2015-05-19
A carriage for high magnetic field environments includes a plurality of work-piece separators disposed in an operable relationship with a work-piece processing magnet having a magnetic field strength of at least 1 Tesla for supporting and separating a plurality of work-pieces by a preselected, essentially equal spacing, so that, as a first work-piece is inserted into the magnetic field, a second work-piece is simultaneously withdrawn from the magnetic field, so that an attractive magnetic force imparted on the first work-piece offsets a resistive magnetic force imparted on the second work-piece.
NASA Astrophysics Data System (ADS)
Greiner, Maximilian; Elts, Ekaterina; Schneider, Julian; Reuter, Karsten; Briesen, Heiko
2014-11-01
The CHARMM, general Amber and OPLS force fields are evaluated for their suitability in simulating the molecular dynamics of the dissolution of the hydrophobic, small-molecule active pharmaceutical ingredients aspirin, ibuprofen, and paracetamol in aqueous media. The force fields are evaluated by comparison with quantum chemical simulations or experimental references on the basis of the following capabilities: accurately representing intra- and intermolecular interactions, appropriately reproducing crystal lattice parameters, adequately describing thermodynamic properties, and the qualitative description of the dissolution behavior. To make this approach easily accessible for evaluating the dissolution properties of novel drug candidates in the early stage of drug development, the force field parameter files are generated using online resources such as the SWISS PARAM servers, and the software packages ACPYPE and Maestro. All force fields are found to reproduce the intermolecular interactions with a reasonable degree of accuracy, with the general Amber and CHARMM force fields showing the best agreement with quantum mechanical calculations. A stable crystal bulk structure is obtained for all model substances, except for ibuprofen, where the reproductions of the lattice parameters and observed crystal stability are considerably poor for all force fields. The heat of solution used to evaluate the solid-to-solution phase transitions is found to be in qualitative agreement with the experimental data for all combinations tested, with the results being quantitatively optimum for the general Amber and CHARMM force fields. For aspirin and paracetamol, stable crystal-water interfaces were obtained. The (100), (110), (011) and (001) interfaces of aspirin or paracetamol and water were simulated for each force field for 30 ns. Although generally expected as a rare event, in some of the simulations, dissolution is observed at 310 K and ambient pressure conditions.
Transferable next-generation force fields from simple liquids to complex materials.
Schmidt, J R; Yu, Kuang; McDaniel, Jesse G
2015-03-17
Molecular simulations have had a transformative impact on chemists' understanding of the structure and dynamics of molecular systems. Simulations can both explain and predict chemical phenomena, and they provide a unique bridge between the microscopic and macroscopic regimes. The input for such simulations is the intermolecular interactions, which then determine the forces on the constituent atoms and therefore the time evolution and equilibrium properties of the system. However, in practice, accuracy and reliability are often limited by the fidelity of the description of those very same interactions, most typically embodied approximately in mathematical form in what are known as force fields. Force fields most often utilize conceptually simple functional forms that have been parametrized to reproduce existing experimental gas phase or bulk data. Yet, reliance on empirical parametrization can sometimes introduce limitations with respect to novel chemical systems or uncontrolled errors when moving to temperatures, pressures, or environments that differ from those for which they were developed. Alternatively, it is possible to develop force fields entirely from first principles, using accurate electronic structure calculations to determine the intermolecular interactions. This introduces a new set of challenges, including the transferability of the resulting force field to related chemical systems. In response, we recently developed an alternative approach to develop force fields entirely from first-principles electronic structure calculations based on intermolecular perturbation theory. Making use of an energy decomposition analysis ensures, by construction, that the resulting force fields contain the correct balance of the various components of intermolecular interaction (exchange repulsion, electrostatics, induction, and dispersion), each treated by a functional form that reflects the underlying physics. We therefore refer to the resulting force fields as
NASA Astrophysics Data System (ADS)
Jenke, Martin Günter; Santschi, Christian; Hoffmann, Patrik
2008-02-01
Accurate simultaneous measurements on the topography and electrostatic force field of 500nm pitch interdigitated electrodes embedded in a thin SiO2 layer in a plane perpendicular to the orientation of the electrodes are shown for the first time. A static force distance curve (FDC) based method has been developed, which allows a lateral and vertical resolution of 25 and 2nm, respectively. The measured force field distribution remains stable as result of the well controlled fabrication procedure of Pt cantilever tips that allows thousands of FDC measurements. A numerical model is established as well which demonstrates good agreement with the experimental results.
Apparatus having reduced mechanical forces for supporting high magnetic fields
Prueitt, Melvin L.; Mueller, Fred M.; Smith, James L.
1991-01-01
The present invention identifies several configurations of conducting elements capable of supporting extremely high magnetic fields suitable for plasma confinement, wherein forces experienced by the conducting elements are significantly reduced over those which are present as a result of the generation of such high fields by conventional techniques. It is anticipated that the use of superconducting materials will both permit the attainment of such high fields and further permit such fields to be generated with vastly improved efficiency.
Balancing the interactions of ions, water, and DNA in the Drude polarizable force field.
Savelyev, Alexey; MacKerell, Alexander D
2014-06-19
Recently we presented a first-generation all-atom Drude polarizable force field for DNA based on the classical Drude oscillator model, focusing on optimization of key dihedral angles followed by extensive validation of the force field parameters. Presently, we describe the procedure for balancing the electrostatic interactions between ions, water, and DNA as required for development of the Drude force field for DNA. The proper balance of these interactions is shown to impact DNA stability and subtler conformational properties, including the conformational equilibrium between the BI and BII states, and the A and B forms of DNA. The parametrization efforts were simultaneously guided by gas-phase quantum mechanics (QM) data on small model compounds and condensed-phase experimental data on the hydration and osmotic properties of biologically relevant ions and their solutions, as well as theoretical predictions for ionic distribution around DNA oligomer. In addition, fine-tuning of the internal base parameters was performed to obtain the final DNA model. Notably, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by the condensed ions as compared to the CHARMM36 additive DNA force field, indicating an improved physical description of the forces dictating the ionic solvation of DNA due to the explicit treatment of electronic polarizability. In combination with the polarizable DNA force field, the availability of Drude polarizable parameters for proteins, lipids, and carbohydrates will allow for simulation studies of heterogeneous biological systems. PMID:24874104
Balancing the Interactions of Ions, Water, and DNA in the Drude Polarizable Force Field
2015-01-01
Recently we presented a first-generation all-atom Drude polarizable force field for DNA based on the classical Drude oscillator model, focusing on optimization of key dihedral angles followed by extensive validation of the force field parameters. Presently, we describe the procedure for balancing the electrostatic interactions between ions, water, and DNA as required for development of the Drude force field for DNA. The proper balance of these interactions is shown to impact DNA stability and subtler conformational properties, including the conformational equilibrium between the BI and BII states, and the A and B forms of DNA. The parametrization efforts were simultaneously guided by gas-phase quantum mechanics (QM) data on small model compounds and condensed-phase experimental data on the hydration and osmotic properties of biologically relevant ions and their solutions, as well as theoretical predictions for ionic distribution around DNA oligomer. In addition, fine-tuning of the internal base parameters was performed to obtain the final DNA model. Notably, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by the condensed ions as compared to the CHARMM36 additive DNA force field, indicating an improved physical description of the forces dictating the ionic solvation of DNA due to the explicit treatment of electronic polarizability. In combination with the polarizable DNA force field, the availability of Drude polarizable parameters for proteins, lipids, and carbohydrates will allow for simulation studies of heterogeneous biological systems. PMID:24874104
NASA Technical Reports Server (NTRS)
El-Kaddah, N.; Szekely, J.
1982-01-01
A mathematical representation was developed for the electromagnetic force field, the flow field, the temperature field (and for transport controlled kinetics), in a levitation melted metal droplet. The technique of mutual inductances was employed for the calculation of the electromagnetic force field, while the turbulent Navier - Stokes equations and the turbulent convective transport equations were used to represent the fluid flow field, the temperature field and the concentration field. The governing differential equations, written in spherical coordinates, were solved numerically. The computed results were in good agreement with measurements, regarding the lifting force, and the average temperature of the specimen and carburization rates, which were transport controlled.
Force field parameter estimation of functional perfluoropolyether lubricants
Smith, Robert; Seung Chung, Pil; Steckel, Janice A.; Jhon, Myung S.; Biegler, Lorenz T.
2011-01-01
The head disk interface in hard disk drive can be considered one of the hierarchical multiscale systems, which require the hybridization of multiscale modeling methods with coarse-graining procedure. However, the fundamental force field parameters are required to enable the coarse-graining procedure from atomistic/molecular scale to mesoscale models .In this paper, we investigate beyond molecular level and perform ab-initio calculations to obtain the force field parameters. Intramolecular force field parameters for the Zdol and Ztetraol were evaluated with truncated PFPE molecules to allow for feasible quantum calculations while still maintaining the characteristic chemical structure of the end groups. Using the harmonic approximation to the bond and angle potentials, the parameters were derived from the Hessian matrix, and the dihedral force constants are fit to the torsional energy profiles generated by a series of constrained molecular geometry optimization.
Force Field Parameter Estimation of Functional Perfluoropolyether Lubricants
Smith, R.; Chung, P.S.; Steckel, J; Jhon, M.S.; Biegler, L.T.
2011-01-01
The head disk interface in hard disk drive can be considered one of the hierarchical multiscale systems, which require the hybridization of multiscale modeling methods with coarse-graining procedure. However, the fundamental force field parameters are required to enable the coarse-graining procedure from atomistic/molecular scale to mesoscale models .In this paper, we investigate beyond molecular level and perform ab-initio calculations to obtain the force field parameters. Intramolecular force field parameters for the Zdol and Ztetraol were evaluated with truncated PFPE molecules to allow for feasible quantum calculations while still maintaining the characteristic chemical structure of the end groups. Using the harmonic approximation to the bond and angle potentials, the parameters were derived from the Hessian matrix, and the dihedral force constants are fit to the torsional energy profiles generated by a series of constrained molecular geometry optimization.
Force Field Parameter Estimation of Functional Perfluoropolyether Lubricants
Smith, R.; Chung, P.S.; Steckel, J; Jhon, M.S.; Biegler, L.T.
2011-01-01
The head disk interface in a hard disk drive can be considered to be one of the hierarchical multiscale systems, which require the hybridization of multiscale modeling methods with coarse-graining procedure. However, the fundamental force field parameters are required to enable the coarse-graining procedure from atomistic/molecular scale to mesoscale models. In this paper, we investigate beyond molecular level and perform ab initio calculations to obtain the force field parameters. Intramolecular force field parameters for Zdol and Ztetraol were evaluated with truncated PFPE molecules to allow for feasible quantum calculations while still maintaining the characteristic chemical structure of the end groups. Using the harmonic approximation to the bond and angle potentials, the parameters were derived from the Hessian matrix, and the dihedral force constants are fit to the torsional energy profiles generated by a series of constrained molecular geometry optimization.
Force Field Parameter Estimation of Functional Perfluoropolyether Lubricants
Smith, R; Chung, P S; Steckel, J A; Jhon, M S; Biegler, L T
2011-01-01
The head disk interface in hard disk drive can be considered one of the hierarchical multiscale systems, which require the hybridization of multiscale modeling methods with coarse-graining procedure. However, the fundamental force field parameters are required to enable the coarse-graining procedure from atomistic/molecular scale to mesoscale models .In this paper, we investigate beyond molecular level and perform ab-initio calculations to obtain the force field parameters. Intramolecular force field parameters for the Zdol and Ztetraolwere evaluated with truncated PFPE molecules to allow for feasible quantum calculations while still maintaining the characteristic chemical structure of the end groups. Using the harmonic approximation to the bond and angle potentials, the parameters were derived from the Hessian matrix, and the dihedral force constants are fit to the torsional energy profiles generated by a series of constrained molecular geometry optimization.
NASA Astrophysics Data System (ADS)
Hackel, Stefan; Montenbruck, Oliver; Steigenberger, -Peter; Eineder, Michael; Gisinger, Christoph
Remote sensing satellites support a broad range of scientific and commercial applications. The two radar imaging satellites TerraSAR-X and TanDEM-X provide spaceborne Synthetic Aperture Radar (SAR) and interferometric SAR data with a very high accuracy. The increasing demand for precise radar products relies on sophisticated validation methods, which require precise and accurate orbit products. Basically, the precise reconstruction of the satellite’s trajectory is based on the Global Positioning System (GPS) measurements from a geodetic-grade dual-frequency receiver onboard the spacecraft. The Reduced Dynamic Orbit Determination (RDOD) approach utilizes models for the gravitational and non-gravitational forces. Following a proper analysis of the orbit quality, systematics in the orbit products have been identified, which reflect deficits in the non-gravitational force models. A detailed satellite macro model is introduced to describe the geometry and the optical surface properties of the satellite. Two major non-gravitational forces are the direct and the indirect Solar Radiation Pressure (SRP). Due to the dusk-dawn orbit configuration of TerraSAR-X, the satellite is almost constantly illuminated by the Sun. Therefore, the direct SRP has an effect on the lateral stability of the determined orbit. The indirect effect of the solar radiation principally contributes to the Earth Radiation Pressure (ERP). The resulting force depends on the sunlight, which is reflected by the illuminated Earth surface in the visible, and the emission of the Earth body in the infrared spectra. Both components of ERP require Earth models to describe the optical properties of the Earth surface. Therefore, the influence of different Earth models on the orbit quality is assessed within the presentation. The presentation highlights the influence of non-gravitational force and satellite macro models on the orbit quality of TerraSAR-X.
NASA Astrophysics Data System (ADS)
Williams, Paul; Simonds, Brian; Sowards, Jeffrey; Hadler, Joshua
2016-03-01
In laser manufacturing operations, accurate measurement of laser power is important for product quality, operational repeatability, and process validation. Accurate real-time measurement of high-power lasers, however, is difficult. Typical thermal power meters must absorb all the laser power in order to measure it. This constrains power meters to be large, slow and exclusive (that is, the laser cannot be used for its intended purpose during the measurement). To address these limitations, we have developed a different paradigm in laser power measurement where the power is not measured according to its thermal equivalent but rather by measuring the laser beam's momentum (radiation pressure). Very simply, light reflecting from a mirror imparts a small force perpendicular to the mirror which is proportional to the optical power. By mounting a high-reflectivity mirror on a high-sensitivity force transducer (scale), we are able to measure laser power in the range of tens of watts up to ~ 100 kW. The critical parameters for such a device are mirror reflectivity, angle of incidence, and scale sensitivity and accuracy. We will describe our experimental characterization of a radiation-pressure-based optical power meter. We have tested it for modulated and CW laser powers up to 92 kW in the laboratory and up to 20 kW in an experimental laser welding booth. We will describe present accuracy, temporal response, sources of measurement uncertainty, and hurdles which must be overcome to have an accurate power meter capable of routine operation as a turning mirror within a laser delivery head.
Free energy simulations with the AMOEBA polarizable force field and metadynamics on GPU platform.
Peng, Xiangda; Zhang, Yuebin; Chu, Huiying; Li, Guohui
2016-03-01
The free energy calculation library PLUMED has been incorporated into the OpenMM simulation toolkit, with the purpose to perform enhanced sampling MD simulations using the AMOEBA polarizable force field on GPU platform. Two examples, (I) the free energy profile of water pair separation (II) alanine dipeptide dihedral angle free energy surface in explicit solvent, are provided here to demonstrate the accuracy and efficiency of our implementation. The converged free energy profiles could be obtained within an affordable MD simulation time when the AMOEBA polarizable force field is employed. Moreover, the free energy surfaces estimated using the AMOEBA polarizable force field are in agreement with those calculated from experimental data and ab initio methods. Hence, the implementation in this work is reliable and would be utilized to study more complicated biological phenomena in both an accurate and efficient way. © 2015 Wiley Periodicals, Inc. PMID:26493154
Development of a True Transition State Force Field from Quantum Mechanical Calculations.
Madarász, Ádám; Berta, Dénes; Paton, Robert S
2016-04-12
Transition state force fields (TSFF) treated the TS structure as an artificial minimum on the potential energy surface in the past decades. The necessary parameters were developed either manually or by the Quantum-to-molecular mechanics method (Q2MM). In contrast with these approaches, here we propose to model the TS structures as genuine saddle points at the molecular mechanics level. Different methods were tested on small model systems of general chemical reactions such as protonation, nucleophilic attack, and substitution, and the new procedure led to more accurate models than the Q2MM-type parametrization. To demonstrate the practicality of our approach, transferrable parameters have been developed for Mo-catalyzed olefin metathesis using quantum mechanical properties as reference data. Based on the proposed strategy, any force field can be extended with true transition state force field (TTSFF) parameters, and they can be readily applied in several molecular mechanics programs as well. PMID:26925858
Reactive Force Field for Liquid Hydrazoic Acid with Applications to Detonation Chemistry
NASA Astrophysics Data System (ADS)
Furman, David; Dubnikova, Faina; van Duin, Adri; Zeiri, Yehuda; Kosloff, Ronnie
The development of a reactive force field (ReaxFF formalism) for Hydrazoic acid (HN3), a highly sensitive liquid energetic material, is reported. The force field accurately reproduces results of density functional theory (DFT) calculations. The quality and performance of the force field are examined by detailed comparison with DFT calculations related to uni, bi and trimolecular thermal decomposition routes. Reactive molecular dynamics (RMD) simulations are performed to reveal the initial chemical events governing the detonation chemistry of liquid HN3. The outcome of these simulations compares very well with recent results of tight-binding DFT molecular dynamics and thermodynamic calculations. Based on our RMD simulations, predictions were made for the activation energies and volumes in a broad range of temperatures and initial material compressions. Work Supported by The Center of Excellence for Explosives Detection, Mitigation and Response, Department of Homeland Security.
Force-Field Compensation in a Manual Tracking Task
Squeri, Valentina; Masia, Lorenzo; Casadio, Maura; Morasso, Pietro; Vergaro, Elena
2010-01-01
This study addresses force/movement control in a dynamic “hybrid” task: the master sub-task is continuous manual tracking of a target moving along an eight-shaped Lissajous figure, with the tracking error as the primary performance index; the slave sub-task is compensation of a disturbing curl viscous field, compatibly with the primary performance index. The two sub-tasks are correlated because the lateral force the subject must exert on the eight-shape must be proportional to the longitudinal movement speed in order to perform a good tracking. The results confirm that visuo-manual tracking is characterized by an intermittent control mechanism, in agreement with previous work; the novel finding is that the overall control patterns are not altered by the presence of a large deviating force field, if compared with the undisturbed condition. It is also found that the control of interaction-forces is achieved by a combination of arm stiffness properties and direct force control, as suggested by the systematic lateral deviation of the trajectories from the nominal path and the comparison between perturbed trials and catch trials. The coordination of the two sub-tasks is quickly learnt after the activation of the deviating force field and is achieved by a combination of force and the stiffness components (about 80% vs. 20%), which is a function of the implicit accuracy of the tracking task. PMID:20567516
Visualization of Force Fields in Protein Structure Prediction
Crawford, Clark; Kreylos, Oliver; Hamann, Bernd; Crivelli, Silvia
2005-04-26
The force fields used in molecular computational biology are not mathematically defined in such a way that their mathematical representation would facilitate the straightforward application of volume visualization techniques. To visualize energy, it is necessary to define a spatial mapping for these fields. Equipped with such a mapping, we can generate volume renderings of the internal energy states in a molecule. We describe our force field, the spatial mapping that we used for energy, and the visualizations that we produced from this mapping. We provide images and animations that offer insight into the computational behavior of the energy optimization algorithms that we employ.
Highly accurate analytical energy of a two-dimensional exciton in a constant magnetic field
NASA Astrophysics Data System (ADS)
Hoang, Ngoc-Tram D.; Nguyen, Duy-Anh P.; Hoang, Van-Hung; Le, Van-Hoang
2016-08-01
Explicit expressions are given for analytically describing the dependence of the energy of a two-dimensional exciton on magnetic field intensity. These expressions are highly accurate with the precision of up to three decimal places for the whole range of the magnetic field intensity. The results are shown for the ground state and some excited states; moreover, we have all formulae to obtain similar expressions of any excited state. Analysis of numerical results shows that the precision of three decimal places is maintained for the excited states with the principal quantum number of up to n=100.
Qiu, Wentao; Ndao, Abdoulaye; Lu, Huihui; Bernal, Maria-Pilar; Baida, Fadi Issam
2016-09-01
We present a theoretical study of guided resonances (GR) on a thin film lithium niobate rectangular lattice photonic crystal by band diagram calculations and 3D Finite Difference Time Domain (FDTD) transmission investigations which cover a broad range of parameters. A photonic crystal with an active zone as small as 13μm×13μm×0.7μm can be easily designed to obtain a resonance Q value in the order of 1000. These resonances are then employed in electric field (E-field) sensing applications exploiting the electro optic (EO) effect of lithium niobate. A local field factor that is calculated locally for each FDTD cell is proposed to accurately estimate the sensitivity of GR based E-field sensor. The local field factor allows well agreement between simulations and reported experimental data therefore providing a valuable method in optimizing the GR structure to obtain high sensitivities. When these resonances are associated with sub-picometer optical spectrum analyzer and high field enhancement antenna design, an E-field probe with a sensitivity of 50 μV/m could be achieved. The results of our simulations could be also exploited in other EO based applications such as EEG (Electroencephalography) or ECG (Electrocardiography) probe and E-field frequency detector with an 'invisible' probe to the field being detected etc. PMID:27607627
Force field dependent solution properties of glycine oligomers
Drake, Justin A.
2015-01-01
Molecular simulations can be used to study disordered polypeptide systems and to generate hypotheses on the underlying structural and thermodynamic mechanisms that govern their function. As the number of disordered protein systems investigated with simulations increase, it is important to understand how particular force fields affect the structural properties of disordered polypeptides in solution. To this end, we performed a comparative structural analysis of Gly3 and Gly10 in aqueous solution from all-atom, microsecond MD simulations using the CHARMM 27 (C27), CHARMM 36 (C36), and Amber ff12SB force fields. For each force field, Gly3 and Gly10 were simulated for at least 300 ns and 1 μs, respectively. Simulating oligoglycines of two different lengths allows us to evaluate how force field effects depend on polypeptide length. Using a variety of structural metrics (e.g. end-to-end distance, radius of gyration, dihedral angle distributions), we characterize the distribution of oligoglycine conformers for each force field and show that each sample conformation space differently, yielding considerably different structural tendencies of the same oligoglycine model in solution. Notably, we find that C36 samples more extended oligoglycine structures than both C27 and ff12SB. PMID:25952623
Force field-dependent solution properties of glycine oligomers.
Drake, Justin A; Pettitt, B Montgomery
2015-06-30
Molecular simulations can be used to study disordered polypeptide systems and to generate hypotheses on the underlying structural and thermodynamic mechanisms that govern their function. As the number of disordered protein systems investigated with simulations increase, it is important to understand how particular force fields affect the structural properties of disordered polypeptides in solution. To this end, we performed a comparative structural analysis of Gly(3) and Gly(10) in aqueous solution from all atom, microsecond molecular dynamics (MD) simulations using the CHARMM 27 (C27), CHARMM 36 (C36), and Amber ff12SB force fields. For each force field, Gly(3) and Gly(10) were simulated for at least 300 ns and 1 μs, respectively. Simulating oligoglycines of two different lengths allows us to evaluate how force field effects depend on polypeptide length. Using a variety of structural metrics (e.g., end-to-end distance, radius of gyration, dihedral angle distributions), we characterize the distribution of oligoglycine conformers for each force field and show that each sample conformation space differently, yielding considerably different structural tendencies of the same oligoglycine model in solution. Notably, we find that C36 samples more extended oligoglycine structures than both C27 and ff12SB. PMID:25952623
Benchmarking of Force Fields for Molecule-Membrane Interactions.
Paloncýová, Markéta; Fabre, Gabin; DeVane, Russell H; Trouillas, Patrick; Berka, Karel; Otyepka, Michal
2014-09-01
Studies of drug-membrane interactions witness an ever-growing interest, as penetration, accumulation, and positioning of drugs play a crucial role in drug delivery and metabolism in human body. Molecular dynamics simulations complement nicely experimental measurements and provide us with new insight into drug-membrane interactions; however, the quality of the theoretical data dramatically depends on the quality of the force field used. We calculated the free energy profiles of 11 molecules through a model dimyristoylphosphatidylcholine (DMPC) membrane bilayer using five force fields, namely Berger, Slipids, CHARMM36, GAFFlipids, and GROMOS 43A1-S3. For the sake of comparison, we also employed the semicontinuous tool COSMOmic. High correlation was observed between theoretical and experimental partition coefficients (log K). Partition coefficients calculated by all-atomic force fields (Slipids, CHARMM36, and GAFFlipids) and COSMOmic differed by less than 0.75 log units from the experiment and Slipids emerged as the best performing force field. This work provides the following recommendations (i) for a global, systematic and high throughput thermodynamic evaluations (e.g., log K) of drugs COSMOmic is a tool of choice due to low computational costs; (ii) for studies of the hydrophilic molecules CHARMM36 should be considered; and (iii) for studies of more complex systems, taking into account all pros and cons, Slipids is the force field of choice. PMID:26588554
Force-free magnetic fields - The magneto-frictional method
NASA Technical Reports Server (NTRS)
Yang, W. H.; Sturrock, P. A.; Antiochos, S. K.
1986-01-01
The problem under discussion is that of calculating magnetic field configurations in which the Lorentz force j x B is everywhere zero, subject to specified boundary conditions. We choose to represent the magnetic field in terms of Clebsch variables in the form B = grad alpha x grad beta. These variables are constant on any field line so that each field line is labeled by the corresponding values of alpha and beta. When the field is described in this way, the most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. We show that such field configurations may be calculated by a magneto-frictional method. We imagine that the field lines move through a stationary medium, and that each element of magnetic field is subject to a frictional force parallel to and opposing the velocity of the field line. This concept leads to an iteration procedure for modifying the variables alpha and beta, that tends asymptotically towards the force-free state. We apply the method first to a simple problem in two rectangular dimensions, and then to a problem of cylindrical symmetry that was previously discussed by Barnes and Sturrock (1972). In one important respect, our new results differ from the earlier results of Barnes and Sturrock, and we conclude that the earlier article was in error.
Physical Limitations of Empirical Field Models: Force Balance and Plasma Pressure
Sorin Zaharia; C.Z. Cheng
2002-06-18
In this paper, we study whether the magnetic field of the T96 empirical model can be in force balance with an isotropic plasma pressure distribution. Using the field of T96, we obtain values for the pressure P by solving a Poisson-type equation {del}{sup 2}P = {del} {center_dot} (J x B) in the equatorial plane, and 1-D profiles on the Sun-Earth axis by integrating {del}P = J x B. We work in a flux coordinate system in which the magnetic field is expressed in terms of Euler potentials. Our results lead to the conclusion that the T96 model field cannot be in equilibrium with an isotropic pressure. We also analyze in detail the computation of Birkeland currents using the Vasyliunas relation and the T96 field, which yields unphysical results, again indicating the lack of force balance in the empirical model. The underlying reason for the force imbalance is likely the fact that the derivatives of the least-square fitted model B are not accurate predictions of the actual magnetospheric field derivatives. Finally, we discuss a possible solution to the problem of lack of force balance in empirical field models.
Polarizable Force Field with a σ-Hole for Liquid and Aqueous Bromomethane.
Adluri, Archita N S; Murphy, Jennifer N; Tozer, Tiffany; Rowley, Christopher N
2015-10-22
Bromomethane (CH3Br) is an acutely toxic environmental pollutant that contributes to ozone depletion. Molecular simulation could be a valuable tool for studying its partitioning and transport in the environment if an accurate molecular model was available. The generalized Amber force field (GAFF), OPLS (optimized potentials for liquid simulations) force field, and CHARMM general force field (CGenFF) were tested for their ability to model the physical properties of liquid bromomethane. The OPLS force field was in fairly good agreement with experiment, while CGenFF and GAFF were significantly in error. The Br Lennard-Jones parameters of the GAFF and CGenFF models were reparameterized, but their radial distribution functions still have significant deviations from those calculated by ab initio molecular dynamics (AIMD). A Drude polarizable force field for bromomethane was parametrized with an off-center positively charged site to represent the C-Br σ-hole. This model is in good agreement with the bulk physical properties and the AIMD RDFs. The modest solubility of bromomethane was reproduced by this model, with dispersion interactions being the dominant water-solute interaction. The water-solute electrostatic interactions are a smaller factor in solubility. This model predicts bromomethane to have a 13 kJ mol(-1) surface excess potential at the water-vapor interface. PMID:26419599
Parameterization and optimization of the menthol force field for molecular dynamics simulations.
Jasik, Mateusz; Szefczyk, Borys
2016-10-01
Menthol's various biological properties render it a useful component for medical and cosmetological applications, while its three centers of asymmetry mean that it can be used in a range of organic reactions. Menthol-substituted ionic liquids (ILs) have been found to exhibit promising antimicrobial and antielectrostatic properties, as well as being useful in organic catalysis and biochemical studies. However, so far, a force field designed and validated specifically for the menthol molecule has not been constructed. In the present work, the validation and optimization of force field parameters with regard to the ability to reproduce the macroscopic properties of menthol is presented. The set of optimized potentials for liquid simulations all atom (OPLS-AA) compatible parameters was tested and carefully tuned. The refinement of parameters included fitting of partial atomic charges, optimization of Lennard-Jones parameters, and recalculation of the dihedral angle parameters needed to reproduce quantum energy profiles. To validate the force field, a variety of physicochemical properties were calculated for liquid menthol. Both thermodynamic and kinetic properties were taken into account, including density, surface tension, enthalpy of vaporization, and shear viscosity. The obtained force field was proven to accurately reproduce the properties of the investigated compound while being fully compatible with the OPLS-AA force field. PMID:27604277
Hölzl, Christoph; Kibies, Patrick; Imoto, Sho; Frach, Roland; Suladze, Saba; Winter, Roland; Marx, Dominik; Horinek, Dominik; Kast, Stefan M
2016-04-14
Accurate force fields are one of the major pillars on which successful molecular dynamics simulations of complex biomolecular processes rest. They have been optimized for ambient conditions, whereas high-pressure simulations become increasingly important in pressure perturbation studies, using pressure as an independent thermodynamic variable. Here, we explore the design of non-polarizable force fields tailored to work well in the realm of kilobar pressures - while avoiding complete reparameterization. Our key is to first compute the pressure-induced electronic and structural response of a solute by combining an integral equation approach to include pressure effects on solvent structure with a quantum-chemical treatment of the solute within the embedded cluster reference interaction site model (EC-RISM) framework. Next, the solute's response to compression is taken into account by introducing pressure-dependence into selected parameters of a well-established force field. In our proof-of-principle study, the full machinery is applied to N,N,N-trimethylamine-N-oxide (TMAO) in water being a potent osmolyte that counteracts pressure denaturation. EC-RISM theory is shown to describe well the charge redistribution upon compression of TMAO(aq) to 10 kbar, which is then embodied in force field molecular dynamics by pressure-dependent partial charges. The performance of the high pressure force field is assessed by comparing to experimental and ab initio molecular dynamics data. Beyond its broad usefulness for designing non-polarizable force fields for extreme thermodynamic conditions, a good description of the pressure-response of solutions is highly recommended when constructing and validating polarizable force fields. PMID:27083705
NASA Astrophysics Data System (ADS)
Hölzl, Christoph; Kibies, Patrick; Imoto, Sho; Frach, Roland; Suladze, Saba; Winter, Roland; Marx, Dominik; Horinek, Dominik; Kast, Stefan M.
2016-04-01
Accurate force fields are one of the major pillars on which successful molecular dynamics simulations of complex biomolecular processes rest. They have been optimized for ambient conditions, whereas high-pressure simulations become increasingly important in pressure perturbation studies, using pressure as an independent thermodynamic variable. Here, we explore the design of non-polarizable force fields tailored to work well in the realm of kilobar pressures - while avoiding complete reparameterization. Our key is to first compute the pressure-induced electronic and structural response of a solute by combining an integral equation approach to include pressure effects on solvent structure with a quantum-chemical treatment of the solute within the embedded cluster reference interaction site model (EC-RISM) framework. Next, the solute's response to compression is taken into account by introducing pressure-dependence into selected parameters of a well-established force field. In our proof-of-principle study, the full machinery is applied to N,N,N-trimethylamine-N-oxide (TMAO) in water being a potent osmolyte that counteracts pressure denaturation. EC-RISM theory is shown to describe well the charge redistribution upon compression of TMAO(aq) to 10 kbar, which is then embodied in force field molecular dynamics by pressure-dependent partial charges. The performance of the high pressure force field is assessed by comparing to experimental and ab initio molecular dynamics data. Beyond its broad usefulness for designing non-polarizable force fields for extreme thermodynamic conditions, a good description of the pressure-response of solutions is highly recommended when constructing and validating polarizable force fields.
Imaging Nanoscale Electromagnetic Near-Field Distributions Using Optical Forces.
Huang, Fei; Tamma, Venkata Ananth; Mardy, Zahra; Burdett, Jonathan; Wickramasinghe, H Kumar
2015-01-01
We demonstrate the application of Atomic Force Microscopy (AFM) for mapping optical near-fields with nanometer resolution, limited only by the AFM probe geometry. By detecting the optical force between a gold coated AFM probe and its image dipole on a glass substrate, we profile the electric field distributions of tightly focused laser beams with different polarizations. The experimentally recorded focal force maps agree well with theoretical predictions based on a dipole-dipole interaction model. We experimentally estimate the aspect ratio of the apex of gold coated AFM probe using only optical forces. We also show that the optical force between a sharp gold coated AFM probe and a spherical gold nanoparticle of radius 15 nm, is indicative of the electric field distribution between the two interacting particles. Photo Induced Force Microscopy (PIFM) allows for background free, thermal noise limited mechanical imaging of optical phenomenon over wide range of wavelengths from Visible to RF with detection sensitivity limited only by AFM performance. PMID:26073331
Force-free field model of ball lightning
NASA Astrophysics Data System (ADS)
Tsui, K. H.
2001-03-01
Due to the nature that the force-free magnetic field, whose current carried by the conducting plasma is everywhere parallel to the magnetic field it generates, is the minimum energy configuration under the constraint of magnetic helicity conservation, ball lightning is considered as a self-organized phenomenon with a plasma fireball immersed in a spherical force-free magnetic field. Since this field does not exert force on the plasma, the plasma pressure, by itself, is in equilibrium with the surrounding environment, and the force-free magnetic field can take on any value without affecting the plasma. Due to this second feature, singular solutions of the magnetic field that are otherwise excluded are allowed, which enable a large amount of energy to be stored to sustain the ball lightning. The singularity is truncated only by the physical limit of current density that a plasma can carry. Scaling the customary soccer-size fireball to larger dimensions could account for day and night sightings of luminous objects in the sky.
Toroidal linear force-free magnetic fields with axial symmetry
NASA Astrophysics Data System (ADS)
Vandas, M.; Romashets, E.
2016-01-01
Aims: Interplanetary magnetic flux ropes are often described as linear force-free fields. To account for their curvature, toroidal configurations must be used. The aim is to find an analytic description of a linear force-free magnetic field of the toroidal geometry in which the cross section of flux ropes can be controlled. Methods: The solution is found as a superposition of fields given by linear force-free cylinders tangential to a generating toroid. The cylindrical field is expressed in a series of terms that are not all cylindrically symmetric. Results: We found the general form of a toroidal linear force-free magnetic field. The field is azimuthally symmetric with respect to the torus axis. It depends on a set of coefficients that enables controlling the flux rope shape (cross section) to some extent. By varying the coefficients, flux ropes with circular and elliptic cross sections were constructed. Numerical comparison suggests that the simple analytic formula for calculating the helicity in toroidal flux ropes of the circular cross section can be used for flux ropes with elliptic cross sections if the minor radius in the formula is set to the geometric mean of the semi-axes of the elliptic cross section.
Interaction Forces Between Multiple Bodies in a Magnetic Field
NASA Technical Reports Server (NTRS)
Joffe, Benjamin
1996-01-01
Some of the results from experiments to determine the interaction forces between multiple bodies in a magnetic field are presented in this paper. It is shown how the force values and the force directions depend on the configuration of the bodies, their relative positions to each other, and the vector of the primary magnetic field. A number of efficient new automatic loading and assembly machines, as well as manipulators and robots, have been created based on the relationship between bodies and magnetic fields. A few of these patented magnetic devices are presented. The concepts involved open a new way to design universal grippers for robot and other kinds of mechanisms for the manipulation of objects. Some of these concepts can be used for space applications.
An implicit divalent counterion force field for RNA molecular dynamics
NASA Astrophysics Data System (ADS)
Henke, Paul S.; Mak, Chi H.
2016-03-01
How to properly account for polyvalent counterions in a molecular dynamics simulation of polyelectrolytes such as nucleic acids remains an open question. Not only do counterions such as Mg2+ screen electrostatic interactions, they also produce attractive intrachain interactions that stabilize secondary and tertiary structures. Here, we show how a simple force field derived from a recently reported implicit counterion model can be integrated into a molecular dynamics simulation for RNAs to realistically reproduce key structural details of both single-stranded and base-paired RNA constructs. This divalent counterion model is computationally efficient. It works with existing atomistic force fields, or coarse-grained models may be tuned to work with it. We provide optimized parameters for a coarse-grained RNA model that takes advantage of this new counterion force field. Using the new model, we illustrate how the structural flexibility of RNA two-way junctions is modified under different salt conditions.
Additive CHARMM force field for naturally occurring modified ribonucleotides.
Xu, You; Vanommeslaeghe, Kenno; Aleksandrov, Alexey; MacKerell, Alexander D; Nilsson, Lennart
2016-04-15
More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all-atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. PMID:26841080
Additive CHARMM force field for naturally occurring modified ribonucleotides
Xu, You; Vanommeslaeghe, Kenno; Aleksandrov, Alexey; MacKerell, Alexander D.
2016-01-01
More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc. PMID:26841080
Transferable force fields for adsorption of small gases in zeolites.
Martin-Calvo, A; Gutiérrez-Sevillano, J J; Parra, J B; Ania, C O; Calero, S
2015-10-01
We provide transferable force fields for oxygen, nitrogen, and carbon monoxide that are able to reproduce experimental adsorption in both pure silica and alumino-substituted zeolites at cryogenic and high temperatures. The force field parameters can be combined with those previously reported for carbon dioxide, methane, and argon, opening the possibility for studying mixtures of interest containing the six components. Using these force field parameters we obtained some adsorption isotherms at cryogenic temperatures that at first sight were in discrepancies with experimental values for certain molecules and structures. We attribute these discrepancies to the sensitiveness of the equipment and to kinetic impedimenta that can lead to erratic results. Additional problems can be found during simulations when extra-framework cations are present in the system as their lack of mobility at low temperatures could lead to kinetic effects that hinder experimental adsorption. PMID:26313242
Water-Silica Force Field for Simulating Nanodevices
Cruz-Chu, Eduardo R.; Aksimentiev, Aleksei; Schulten, Klaus
2008-01-01
Amorphous silica is an inorganic material that is central for many nanotechnology appplications, such as nanoelectronics, microfluidics, and nanopore technology. In order to use molecular dynamics (MD) simulations to study the behavior of biomolecules with silica, we developed a force field for amorphous silica surfaces based on their macroscopic wetting properties that is compatible with the CHARMM force field and TIP3P water model. The contact angle of a water droplet with silica served as a criterion to tune the intermolecular interactions. The resulting force field was used to study the permeation of water through silica nanopores, illustrating the influence of the surface topography and the intermolecular parameters on permeation kinetics. We find that minute modeling of the amorphous surface is critical for MD studies, since the particular arrangement of surface atoms controls sensitively electrostatic interactions between silica and water. PMID:17064100
2012-01-01
Atmospheric humidity strongly influences the interactions between dry granular particles in process containers. To reduce the energy loss in industrial production processes caused by particle agglomeration, a basic understanding of the dependence of particle interactions on humidity is necessary. Hence, in this study, molecular dynamic simulations were carried out to calculate the adhesion between silica surfaces in the presence of adsorbed water. For a realistic description, the choice of force field is crucial. Because of their frequent use and transferability to biochemical systems, the Clay and CWCA force fields were investigated with respect to their ability to describe the water–silica interface in comparison to the more advanced Reax force field, ab initio calculations, and experiments. PMID:23378869
Scattered field generation and optical forces in transformation optics
NASA Astrophysics Data System (ADS)
Novitsky, A. V.
2016-04-01
In this paper we develop an approach for making various scattered electromagnetic fields on the transformation-optics ground. To do so, we use the a special coordinate transformation from the a vacuum virtual space to physical space, which changes the boundary of the scattering device upon transformation. We explore this approach for small scatterers compared with radiation wavelength, which allows us to predict the arbitrarily directed optical forces. Obtaining scattered fields and optical forces can be useful in nano-optics and optical micromanipulation.
WIPPER: an accurate and efficient field phenotyping platform for large-scale applications
Utsushi, Hiroe; Abe, Akira; Tamiru, Muluneh; Ogasawara, Yumiko; Obara, Tsutomu; Sato, Emiko; Ochiai, Yusuke; Terauchi, Ryohei; Takagi, Hiroki
2015-01-01
More accurate, rapid, and easy phenotyping tools are required to match the recent advances in high-throughput genotyping for accelerating breeding and genetic analysis. The conventional data recording in field notebooks and then inputting data to computers for further analysis is inefficient, time-consuming, laborious, and prone to human error. Here, we report WIPPER (for Wireless Plant Phenotyper), a new phenotyping platform that combines field phenotyping and data recording with the aid of Bluetooth communication, thus saving time and labor not only for field data recoding but also for inputting data to computers. Additionally, it eliminates the risk of human error associated with phenotyping and inputting data. We applied WIPPER to 100 individuals of a rice recombinant inbred line (RIL) for measuring leaf width and relative chlorophyll content (SPAD value), and were able to record an accurate data in a significantly reduced time compared with the conventional method of data collection. We are currently using WIPPER for routine management of rice germplasm including recording and documenting information on phenotypic data, seeds, and DNA for their accelerated utilization in crop breeding. PMID:26175626
Mitigated-force carriage for high magnetic field environments
Ludtka, Gerard M; Ludtka, Gail M; Wilgen, John B; Murphy, Bart L
2014-05-20
A carriage for high magnetic field environments includes a first work-piece holding means for holding a first work-piece, the first work-piece holding means being disposed in an operable relationship with a work-piece processing magnet having a magnetic field strength of at least 1 Tesla. The first work-piece holding means is further disposed in operable connection with a second work-piece holding means for holding a second work-piece so that, as the first work-piece is inserted into the magnetic field, the second work-piece is simultaneously withdrawn from the magnetic field, so that an attractive magnetic force imparted on the first work-piece offsets a resistive magnetic force imparted on the second work-piece.
Energy buildup in sheared force-free magnetic fields
NASA Technical Reports Server (NTRS)
Wolfson, Richard; Low, Boon C.
1992-01-01
Photospheric displacement of the footpoints of solar magnetic field lines results in shearing and twisting of the field, and consequently in the buildup of electric currents and magnetic free energy in the corona. The sudden release of this free energy may be the origin of eruptive events like coronal mass ejections, prominence eruptions, and flares. An important question is whether such an energy release may be accompanied by the opening of magnetic field lines that were previously closed, for such open field lines can provide a route for matter frozen into the field to escape the sun altogether. This paper presents the results of numerical calculations showing that opening of the magnetic field is permitted energetically, in that it is possible to build up more free energy in a sheared, closed, force-free magnetic field than is in a related magnetic configuration having both closed and open field lines. Whether or not the closed force-free field attains enough energy to become partially open depends on the form of the shear profile; the results presented compare the energy buildup for different shear profiles. Implications for solar activity are discussed briefly.
Comparison of different force fields for the study of disaccharides
Technology Transfer Automated Retrieval System (TEKTRAN)
Eighteen empirical force fields and the semi-empirical quantum method PM3CARB-1 were compared for studying ß-cellobiose, a-maltose, and a-galabiose [a-D-Galp-(1'4)-a-D-Galp]. For each disaccharide, the energies of 54 conformers with differing hydroxymethyl, hydroxyl and glycosidic linkage orientatio...
Frequency-dependent force fields for QMMM calculations.
Harczuk, Ignat; Vahtras, Olav; Ågren, Hans
2015-03-28
We outline the construction of frequency-dependent polarizable force fields. The force fields are derived from analytic response theory for different frequencies using a generalization of the LoProp algorithm giving a decomposition of a molecular dynamical polarizability to localized atomic dynamical polarizabilities. These force fields can enter in a variety of applications - we focus on two such applications in this work: firstly, they can be incorporated in a physical, straightforward, way for current existing methods that use polarizable embeddings, and we can show, for the first time, the effect of the frequency dispersion within the classical environment of a quantum mechanics-molecular mechanics (QMMM) method. Our methodology is here evaluated for some test cases comprising water clusters and organic residues. Secondly, together with a modified Silberstein-Applequist procedure for interacting inducible point-dipoles, these frequency-dependent polarizable force fields can be used for a classical determination of frequency-dependent cluster polarizabilities. We evaluate this methodology by comparing with the corresponding results obtained from quantum mechanics or QMMM where the absolute mean [small alpha, Greek, macron] is determined with respect to the size of the QM and MM parts of the total system. PMID:25714984
Symmetrization of the AMBER and CHARMM force fields.
Małolepsza, Edyta; Strodel, Birgit; Khalili, Mey; Trygubenko, Semen; Fejer, Szilard N; Wales, David J
2010-05-01
The AMBER and CHARMM force fields are analyzed from the viewpoint of the permutational symmetry of the potential for feasible exchanges of identical atoms and chemical groups in amino and nucleic acids. In each case, we propose schemes for symmetrizing the potentials, which greatly facilitate the bookkeeping associated with constructing kinetic transition networks via geometry optimization. PMID:20082393
Generating distributed forcing fields for spatial hydrologic modeling
Technology Transfer Automated Retrieval System (TEKTRAN)
Spatial hydrologic modeling requires the development of distributed forcing fields of weather and precipitation. This is particularly difficult in mountainous regions of the western US, where measurement sites are limited and the landscape is dominated by complex terrain and variations in vegetatio...
ON THE FORCE-FREE NATURE OF PHOTOSPHERIC SUNSPOT MAGNETIC FIELDS AS OBSERVED FROM HINODE (SOT/SP)
Tiwari, Sanjiv Kumar
2012-01-01
A magnetic field is force-free if there is no interaction between it and the plasma in the surrounding atmosphere, i.e., electric currents are aligned with the magnetic field, giving rise to zero Lorentz force. The computation of various magnetic parameters, such as magnetic energy (using the virial theorem), gradient of twist of sunspot magnetic fields (computed from the force-free parameter {alpha}), and any kind of extrapolation, heavily hinges on the force-free approximation of the photospheric sunspot magnetic fields. Thus, it is of vital importance to inspect the force-free behavior of sunspot magnetic fields. The force-free nature of sunspot magnetic fields has been examined earlier by some researchers, ending with incoherent results. Accurate photospheric vector field measurements with high spatial resolution are required to inspect the force-free nature of sunspots. For this purpose, we use several vector magnetograms of high spatial resolution obtained from the Solar Optical Telescope/Spectro-Polarimeter on board Hinode. Both the necessary and sufficient conditions for force-free nature are examined by checking the global and local nature of equilibrium magnetic forces over sunspots. We find that sunspot magnetic fields are not very far from the force-free configuration, although they are not completely force-free on the photosphere. The umbral and inner penumbral fields are more force-free than the middle and outer penumbral fields. During their evolution, sunspot magnetic fields are found to maintain their proximity to force-free field behavior. Although a dependence of net Lorentz force components is seen on the evolutionary stages of the sunspots, we do not find a systematic relationship between the nature of sunspot magnetic fields and the associated flare activity. Further, we examine whether the fields at the photosphere follow linear or nonlinear force-free conditions. After examining this in various complex and simple sunspots, we conclude that
The effect of gravitational tidal forces on renormalized quantum fields
NASA Astrophysics Data System (ADS)
Hollowood, Timothy J.; Shore, Graham M.
2012-02-01
The effect of gravitational tidal forces on renormalized quantum fields propagating in curved spacetime is investigated and a generalisation of the optical theorem to curved spacetime is proved. In the case of QED, the interaction of tidal forces with the vacuum polarization cloud of virtual e + e - pairs dressing the renormalized photon has been shown to produce several novel phenomena. In particular, the photon field amplitude can locally increase as well as decrease, corresponding to a negative imaginary part of the refractive index, in apparent violation of unitarity and the optical theorem. Below threshold decays into e + e - pairs may also occur. In this paper, these issues are studied from the point of view of a non-equilibrium initial-value problem, with the field evolution from an initial null surface being calculated for physically distinct initial conditions and for both scalar field theories and QED. It is shown how a generalised version of the optical theorem, valid in curved spacetime, allows a local increase in amplitude while maintaining consistency with unitarity. The picture emerges of the field being dressed and undressed as it propagates through curved spacetime, with the local gravitational tidal forces determining the degree of dressing and hence the amplitude of the renormalized quantum field. These effects are illustrated with many examples, including a description of the undressing of a photon in the vicinity of a black hole singularity.
Generation of accurate integral surfaces in time-dependent vector fields.
Garth, Christoph; Krishnan, Han; Tricoche, Xavier; Bobach, Tom; Joy, Kenneth I
2008-01-01
We present a novel approach for the direct computation of integral surfaces in time-dependent vector fields. As opposed to previous work, which we analyze in detail, our approach is based on a separation of integral surface computation into two stages: surface approximation and generation of a graphical representation. This allows us to overcome several limitations of existing techniques. We first describe an algorithm for surface integration that approximates a series of time lines using iterative refinement and computes a skeleton of the integral surface. In a second step, we generate a well-conditioned triangulation. Our approach allows a highly accurate treatment of very large time-varying vector fields in an efficient, streaming fashion. We examine the properties of the presented methods on several example datasets and perform a numerical study of its correctness and accuracy. Finally, we investigate some visualization aspects of integral surfaces. PMID:18988990
Accurate, efficient, and (iso)geometrically flexible collocation methods for phase-field models
NASA Astrophysics Data System (ADS)
Gomez, Hector; Reali, Alessandro; Sangalli, Giancarlo
2014-04-01
We propose new collocation methods for phase-field models. Our algorithms are based on isogeometric analysis, a new technology that makes use of functions from computational geometry, such as, for example, Non-Uniform Rational B-Splines (NURBS). NURBS exhibit excellent approximability and controllable global smoothness, and can represent exactly most geometries encapsulated in Computer Aided Design (CAD) models. These attributes permitted us to derive accurate, efficient, and geometrically flexible collocation methods for phase-field models. The performance of our method is demonstrated by several numerical examples of phase separation modeled by the Cahn-Hilliard equation. We feel that our method successfully combines the geometrical flexibility of finite elements with the accuracy and simplicity of pseudo-spectral collocation methods, and is a viable alternative to classical collocation methods.
NASA Astrophysics Data System (ADS)
Zacharias, Panagiotis P.; Chatzineofytou, Elpida G.; Spantideas, Sotirios T.; Capsalis, Christos N.
2016-07-01
In the present work, the determination of the magnetic behavior of localized magnetic sources from near-field measurements is examined. The distance power law of the magnetic field fall-off is used in various cases to accurately predict the magnetic signature of an equipment under test (EUT) consisting of multiple alternating current (AC) magnetic sources. Therefore, parameters concerning the location of the observation points (magnetometers) are studied towards this scope. The results clearly show that these parameters are independent of the EUT's size and layout. Additionally, the techniques developed in the present study enable the placing of the magnetometers close to the EUT, thus achieving high signal-to-noise ratio (SNR). Finally, the proposed method is verified by real measurements, using a mobile phone as an EUT.
NASA Astrophysics Data System (ADS)
Bucci, Ovidio M.; Gennarelli, Claudio; Savarese, Catello
1991-01-01
An optimal sampling interpolation algorithm which allows the accurate recovery of plane-rectangular near-field samples from the knowledge of the plane-polar ones is developed. This enables the standard near field-far field (NF-FF) transformation, which takes full advantage of the FFT algorithm, to be applied to plane-polar scanning. The maximum allowable sample spacing is also rigorously derived, and it is shown that it can be significantly greater than lambda/2 as the measurement place moves away from the source. This allows a remarkable reduction of both measurement time and memory storage requirements. The sampling approach is compared with that based on the bivariate Lagrange interpolation (BLI) method. The sampling reconstruction agrees with the exact results significantly better than the BLI, in spite of the significantly lower number of required measurements.
GLYCAM06: a generalizable biomolecular force field. Carbohydrates.
Kirschner, Karl N; Yongye, Austin B; Tschampel, Sarah M; González-Outeiriño, Jorge; Daniels, Charlisa R; Foley, B Lachele; Woods, Robert J
2008-03-01
A new derivation of the GLYCAM06 force field, which removes its previous specificity for carbohydrates, and its dependency on the AMBER force field and parameters, is presented. All pertinent force field terms have been explicitly specified and so no default or generic parameters are employed. The new GLYCAM is no longer limited to any particular class of biomolecules, but is extendible to all molecular classes in the spirit of a small-molecule force field. The torsion terms in the present work were all derived from quantum mechanical data from a collection of minimal molecular fragments and related small molecules. For carbohydrates, there is now a single parameter set applicable to both alpha- and beta-anomers and to all monosaccharide ring sizes and conformations. We demonstrate that deriving dihedral parameters by fitting to QM data for internal rotational energy curves for representative small molecules generally leads to correct rotamer populations in molecular dynamics simulations, and that this approach removes the need for phase corrections in the dihedral terms. However, we note that there are cases where this approach is inadequate. Reported here are the basic components of the new force field as well as an illustration of its extension to carbohydrates. In addition to reproducing the gas-phase properties of an array of small test molecules, condensed-phase simulations employing GLYCAM06 are shown to reproduce rotamer populations for key small molecules and representative biopolymer building blocks in explicit water, as well as crystalline lattice properties, such as unit cell dimensions, and vibrational frequencies. PMID:17849372
GLYCAM06: A Generalizable Biomolecular Force Field. Carbohydrates
KIRSCHNER, KARL N.; YONGYE, AUSTIN B.; TSCHAMPEL, SARAH M.; GONZÁLEZ-OUTEIRIÑO, JORGE; DANIELS, CHARLISA R.; FOLEY, B. LACHELE; WOODS, ROBERT J.
2015-01-01
A new derivation of the GLYCAM06 force field, which removes its previous specificity for carbohydrates, and its dependency on the AMBER force field and parameters, is presented. All pertinent force field terms have been explicitly specified and so no default or generic parameters are employed. The new GLYCAM is no longer limited to any particular class of biomolecules, but is extendible to all molecular classes in the spirit of a small-molecule force field. The torsion terms in the present work were all derived from quantum mechanical data from a collection of minimal molecular fragments and related small molecules. For carbohydrates, there is now a single parameter set applicable to both α- and β-anomers and to all monosaccharide ring sizes and conformations. We demonstrate that deriving dihedral parameters by fitting to QM data for internal rotational energy curves for representative small molecules generally leads to correct rotamer populations in molecular dynamics simulations, and that this approach removes the need for phase corrections in the dihedral terms. However, we note that there are cases where this approach is inadequate. Reported here are the basic components of the new force field as well as an illustration of its extension to carbohydrates. In addition to reproducing the gas-phase properties of an array of small test molecules, condensed-phase simulations employing GLYCAM06 are shown to reproduce rotamer populations for key small molecules and representative biopolymer building blocks in explicit water, as well as crystalline lattice properties, such as unit cell dimensions, and vibrational frequencies. PMID:17849372
Accurate 2d finite element calculations for hydrogen in magnetic fields of arbitrary strength
NASA Astrophysics Data System (ADS)
Schimeczek, C.; Wunner, G.
2014-02-01
Recent observations of hundreds of hydrogen-rich magnetic white dwarf stars with magnetic fields up to 105 T (103 MG) have called for more comprehensive and accurate databases for wavelengths and oscillator strengths of the H atom in strong magnetic fields for all states evolving from the field-free levels with principal quantum numbers n≤10. We present a code to calculate the energy eigenvalues and wave functions of such states which is capable of covering the entire regime of field strengths B=0 T to B˜109 T. We achieve this high flexibility by using a two-dimensional finite element expansion of the wave functions in terms of B-splines in the directions parallel and perpendicular to the magnetic field, instead of using asymptotically valid basis expansions in terms of spherical harmonics or Landau orbitals. We have paid special attention to the automation of the program such that the data points for the magnetic field strengths at which the energy of a given state are calculated can be selected automatically. Furthermore, an elaborate method for varying the basis parameters is applied to ensure that the results reach a pre-selected precision, which also can be adjusted freely. Energies and wave functions are stored in a convenient format for further analysis, e.g. for the calculation of transition energies and oscillator strengths. The code has been tested to work for 300 states with an accuracy of better than 10-6 Rydberg across several symmetry subspaces over the entire regime of magnetic field strengths.
Direct measurement of optical force induced by near-field plasmonic cavity using dynamic mode AFM
Guan, Dongshi; Hang, Zhi Hong; Marset, Zsolt; Liu, Hui; Kravchenko, Ivan I.; Chan, Ho Bun; Chan, C. T.; Tong, Penger
2015-11-20
Plasmonic nanostructures have attracted much attention in recent years because of their potential applications in optical manipulation through near-field enhancement. Continuing experimental efforts have been made to develop accurate techniques to directly measure the near-field optical force induced by the plasmonic nanostructures in the visible frequency range. In this work, we report a new application of dynamic mode atomic force microscopy (DM-AFM) in the measurement of the enhanced optical force acting on a nano-structured plasmonic resonant cavity. The plasmonic cavity is made of an upper gold-coated glass sphere and a lower quartz substrate patterned with an array of subwavelength gold disks. In the near-field when the sphere is positioned close to the disk array, plasmonic resonance is excited in the cavity and the induced force by a 1550 nm infrared laser is found to be increased by an order of magnitude compared with the photon pressure generated by the same laser light. Lastly, the experiment demonstrates that DM-AFM is a powerful tool for the study of light induced forces and their enhancement in plasmonic nanostructures.
Direct measurement of optical force induced by near-field plasmonic cavity using dynamic mode AFM
Guan, Dongshi; Hang, Zhi Hong; Marset, Zsolt; Liu, Hui; Kravchenko, Ivan I.; Chan, Ho Bun; Chan, C. T.; Tong, Penger
2015-11-20
Plasmonic nanostructures have attracted much attention in recent years because of their potential applications in optical manipulation through near-field enhancement. Continuing experimental efforts have been made to develop accurate techniques to directly measure the near-field optical force induced by the plasmonic nanostructures in the visible frequency range. In this work, we report a new application of dynamic mode atomic force microscopy (DM-AFM) in the measurement of the enhanced optical force acting on a nano-structured plasmonic resonant cavity. The plasmonic cavity is made of an upper gold-coated glass sphere and a lower quartz substrate patterned with an array of subwavelength goldmore » disks. In the near-field when the sphere is positioned close to the disk array, plasmonic resonance is excited in the cavity and the induced force by a 1550 nm infrared laser is found to be increased by an order of magnitude compared with the photon pressure generated by the same laser light. Lastly, the experiment demonstrates that DM-AFM is a powerful tool for the study of light induced forces and their enhancement in plasmonic nanostructures.« less
Direct Measurement of Optical Force Induced by Near-Field Plasmonic Cavity Using Dynamic Mode AFM.
Guan, Dongshi; Hang, Zhi Hong; Marcet, Zsolt; Liu, Hui; Kravchenko, I I; Chan, C T; Chan, H B; Tong, Penger
2015-01-01
Plasmonic nanostructures have attracted much attention in recent years because of their potential applications in optical manipulation through near-field enhancement. Continuing experimental efforts have been made to develop accurate techniques to directly measure the near-field optical force induced by the plasmonic nanostructures in the visible frequency range. In this work, we report a new application of dynamic mode atomic force microscopy (DM-AFM) in the measurement of the enhanced optical force acting on a nano-structured plasmonic resonant cavity. The plasmonic cavity is made of an upper gold-coated glass sphere and a lower quartz substrate patterned with an array of subwavelength gold disks. In the near-field when the sphere is positioned close to the disk array, plasmonic resonance is excited in the cavity and the induced force by a 1550 nm infrared laser is found to be increased by an order of magnitude compared with the photon pressure generated by the same laser light. The experiment demonstrates that DM-AFM is a powerful tool for the study of light induced forces and their enhancement in plasmonic nanostructures. PMID:26586455
Direct Measurement of Optical Force Induced by Near-Field Plasmonic Cavity Using Dynamic Mode AFM
Guan, Dongshi; Hang, Zhi Hong; Marcet, Zsolt; Liu, Hui; Kravchenko, I. I.; Chan, C. T.; Chan, H. B.; Tong, Penger
2015-01-01
Plasmonic nanostructures have attracted much attention in recent years because of their potential applications in optical manipulation through near-field enhancement. Continuing experimental efforts have been made to develop accurate techniques to directly measure the near-field optical force induced by the plasmonic nanostructures in the visible frequency range. In this work, we report a new application of dynamic mode atomic force microscopy (DM-AFM) in the measurement of the enhanced optical force acting on a nano-structured plasmonic resonant cavity. The plasmonic cavity is made of an upper gold-coated glass sphere and a lower quartz substrate patterned with an array of subwavelength gold disks. In the near-field when the sphere is positioned close to the disk array, plasmonic resonance is excited in the cavity and the induced force by a 1550 nm infrared laser is found to be increased by an order of magnitude compared with the photon pressure generated by the same laser light. The experiment demonstrates that DM-AFM is a powerful tool for the study of light induced forces and their enhancement in plasmonic nanostructures. PMID:26586455
Direct Measurement of Optical Force Induced by Near-Field Plasmonic Cavity Using Dynamic Mode AFM
NASA Astrophysics Data System (ADS)
Guan, Dongshi; Hang, Zhi Hong; Marcet, Zsolt; Liu, Hui; Kravchenko, I. I.; Chan, C. T.; Chan, H. B.; Tong, Penger
2015-11-01
Plasmonic nanostructures have attracted much attention in recent years because of their potential applications in optical manipulation through near-field enhancement. Continuing experimental efforts have been made to develop accurate techniques to directly measure the near-field optical force induced by the plasmonic nanostructures in the visible frequency range. In this work, we report a new application of dynamic mode atomic force microscopy (DM-AFM) in the measurement of the enhanced optical force acting on a nano-structured plasmonic resonant cavity. The plasmonic cavity is made of an upper gold-coated glass sphere and a lower quartz substrate patterned with an array of subwavelength gold disks. In the near-field when the sphere is positioned close to the disk array, plasmonic resonance is excited in the cavity and the induced force by a 1550 nm infrared laser is found to be increased by an order of magnitude compared with the photon pressure generated by the same laser light. The experiment demonstrates that DM-AFM is a powerful tool for the study of light induced forces and their enhancement in plasmonic nanostructures.
Automatic molecular structure perception for the universal force field.
Artemova, Svetlana; Jaillet, Léonard; Redon, Stephane
2016-05-15
The Universal Force Field (UFF) is a classical force field applicable to almost all atom types of the periodic table. Such a flexibility makes this force field a potential good candidate for simulations involving a large spectrum of systems and, indeed, UFF has been applied to various families of molecules. Unfortunately, initializing UFF, that is, performing molecular structure perception to determine which parameters should be used to compute the UFF energy and forces, appears to be a difficult problem. Although many perception methods exist, they mostly focus on organic molecules, and are thus not well-adapted to the diversity of systems potentially considered with UFF. In this article, we propose an automatic perception method for initializing UFF that includes the identification of the system's connectivity, the assignment of bond orders as well as UFF atom types. This perception scheme is proposed as a self-contained UFF implementation integrated in a new module for the SAMSON software platform for computational nanoscience (http://www.samson-connect.net). We validate both the automatic perception method and the UFF implementation on a series of benchmarks. PMID:26927616
Atomistic force field for alumina fit to density functional theory
Sarsam, Joanne; Finnis, Michael W.; Tangney, Paul
2013-11-28
We present a force field for bulk alumina (Al{sub 2}O{sub 3}), which has been parametrized by fitting the energies, forces, and stresses of a large database of reference configurations to those calculated with density functional theory (DFT). We use a functional form that is simpler and computationally more efficient than some existing models of alumina parametrized by a similar technique. Nevertheless, we demonstrate an accuracy of our potential that is comparable to those existing models and to DFT. We present calculations of crystal structures and energies, elastic constants, phonon spectra, thermal expansion, and point defect formation energies.
NASA Astrophysics Data System (ADS)
Brilliantov, N. V.; Seidel, C.
2012-01-01
The response of grafted polyelectrolytes to electrostatic field that favors adsorption is studied theoretically and by means of molecular dynamics (MD). Applying to the free chain end a constant force which counteracts adsorption, we analyze the size of the desorbed part as a function of force and electric field. Simulations with different loads, corresponding to linear, nonlinear and Hertzian springs, applied to the free end have been also performed to explore the generation of mechanical force by electric field. The MD results are in good agreement with the predictions of theory.
Particle energization in a chaotic force-free magnetic field
NASA Astrophysics Data System (ADS)
Li, Xiaocan; Li, Gang; Dasgupta, Brahmananda
2015-04-01
A force-free field (FFF) is believed to be a reasonable description of the solar corona and in general a good approximation for low-beta plasma. The equations describing the magnetic field of FFF is similar to the ABC fluid equations which has been demonstrated to be chaotic. This implies that charged particles will experience chaotic magnetic field in the corona. Here, we study particle energization in a time-dependent FFF using a test particle approach. An inductive electric field is introduced by turbulent motions of plasma parcels. We find efficient particle acceleration with power-law like particle energy spectra. The power-law indices depend on the amplitude of plasma parcel velocity field and the spatial scales of the magnetic field fluctuation. The spectra are similar for different particle species. This model provide a possible mechanism for seed population generation for particle acceleration by, e.g., CME-driven shocks. Generalization of our results to certain non-force-free-field (NFFF) is straightforward as the sum of two or multiple FFFs naturally yield NFFF.
Parametrization of an Orbital-Based Linear-Scaling Quantum Force Field for Noncovalent Interactions
2015-01-01
We parametrize a linear-scaling quantum mechanical force field called mDC for the accurate reproduction of nonbonded interactions. We provide a new benchmark database of accurate ab initio interactions between sulfur-containing molecules. A variety of nonbond databases are used to compare the new mDC method with other semiempirical, molecular mechanical, ab initio, and combined semiempirical quantum mechanical/molecular mechanical methods. It is shown that the molecular mechanical force field significantly and consistently reproduces the benchmark results with greater accuracy than the semiempirical models and our mDC model produces errors twice as small as the molecular mechanical force field. The comparisons between the methods are extended to the docking of drug candidates to the Cyclin-Dependent Kinase 2 protein receptor. We correlate the protein–ligand binding energies to their experimental inhibition constants and find that the mDC produces the best correlation. Condensed phase simulation of mDC water is performed and shown to produce O–O radial distribution functions similar to TIP4P-EW. PMID:24803856
NASA Astrophysics Data System (ADS)
Leuty, Gary; Berry, Rajiv; Muratore, Christopher; Varshney, Vikas; Turner, Heath
Transition metal dichalcogenides (TMDs) such as molybdenum disulfide (MoS2) have garnered significant interest in recent years. With a layered structure similar to graphene, TMDs also have an intrinsic band gap. This band gap makes them an attractive alternative to graphene in many applications. MoS2 in particular has received attention due to the placement and tenability of its band gap, via functionalization, mechanical manipulation or physisorption. The latter of these is of interest in biosensor devices. Such applications are dependent on understanding physisorption on the MoS2 surface at the molecular level. This can be difficult experimentally but is possible via computer simulation techniques such as molecular dynamics (MD) simulations. MD simulations, however, require a force field accurate to the process modeled. Such a force field must correctly describe non-bonded interactions between substrate layers and between the surface and adsorbates. The force fields we are aware of have focused on intra-layer covalent bonding for structural and vibrational analysis. This work seeks to develop, through DFT and MD simulations with experimental characterization of surface adsorption, a more accurate parameterization for non-bonded interactions for MoS2.
Exploring Adsorption of Water and Ions on Carbon Surfaces using a Polarizable Force Field
Schyman, Patric; Jorgensen, William L.
2013-01-01
Graphene, carbon nanotubes, and fullerenes are of great interest due to their unique properties and diverse applications in biology, molecular electronics, and materials science. Therefore, there is demand for methods that can accurately model the interface between carbon surfaces and their environment. In this letter we compare results for complexes of water, potassium ion, and chloride ion with graphene, carbon nanotube, and fullerene surfaces using a standard non-polarizable force field (OPLS-AA), a polarizable force field (OPLS-AAP), DFT, and ab initio theory. For interactions with water, OPLS-AA with the TIP3P or TIP4P water models describes the interactions with benzene (C6H6) and coronene (C24H12) well; however, for acenes larger than circumcoronene (C54H18) and especially for C60, the interaction energies are somewhat too weak and polarization is needed. For ions interacting with carbon surfaces, inclusion of polarization is essential, and OPLS-AAP is found to perform well in comparison to the highest-level quantum mechanical methods. Overall, OPLS-AAP provides an accurate and computationally efficient force field for modeling condensed-phase systems featuring carbon surfaces. PMID:23440601
Field measurement of basal forces generated by erosive debris flows
McCoy, S.W.; Tucker, G.E.; Kean, J.W.; Coe, J.A.
2013-01-01
It has been proposed that debris flows cut bedrock valleys in steeplands worldwide, but field measurements needed to constrain mechanistic models of this process remain sparse due to the difficulty of instrumenting natural flows. Here we present and analyze measurements made using an automated sensor network, erosion bolts, and a 15.24 cm by 15.24 cm force plate installed in the bedrock channel floor of a steep catchment. These measurements allow us to quantify the distribution of basal forces from natural debris‒flow events that incised bedrock. Over the 4 year monitoring period, 11 debris‒flow events scoured the bedrock channel floor. No clear water flows were observed. Measurements of erosion bolts at the beginning and end of the study indicated that the bedrock channel floor was lowered by 36 to 64 mm. The basal force during these erosive debris‒flow events had a large‒magnitude (up to 21 kN, which was approximately 50 times larger than the concurrent time‒averaged mean force), high‒frequency (greater than 1 Hz) fluctuating component. We interpret these fluctuations as flow particles impacting the bed. The resulting variability in force magnitude increased linearly with the time‒averaged mean basal force. Probability density functions of basal normal forces were consistent with a generalized Pareto distribution, rather than the exponential distribution that is commonly found in experimental and simulated monodispersed granular flows and which has a lower probability of large forces. When the bed sediment thickness covering the force plate was greater than ~ 20 times the median bed sediment grain size, no significant fluctuations about the time‒averaged mean force were measured, indicating that a thin layer of sediment (~ 5 cm in the monitored cases) can effectively shield the subjacent bed from erosive impacts. Coarse‒grained granular surges and water‒rich, intersurge flow had very similar basal force distributions despite
2013-01-01
The accelerated molecular dynamics (aMD) method has recently been shown to enhance the sampling of biomolecules in molecular dynamics (MD) simulations, often by several orders of magnitude. Here, we describe an implementation of the aMD method for the OpenMM application layer that takes full advantage of graphics processing units (GPUs) computing. The aMD method is shown to work in combination with the AMOEBA polarizable force field (AMOEBA-aMD), allowing the simulation of long time-scale events with a polarizable force field. Benchmarks are provided to show that the AMOEBA-aMD method is efficiently implemented and produces accurate results in its standard parametrization. For the BPTI protein, we demonstrate that the protein structure described with AMOEBA remains stable even on the extended time scales accessed at high levels of accelerations. For the DNA repair metalloenzyme endonuclease IV, we show that the use of the AMOEBA force field is a significant improvement over fixed charged models for describing the enzyme active-site. The new AMOEBA-aMD method is publicly available (http://wiki.simtk.org/openmm/VirtualRepository) and promises to be interesting for studying complex systems that can benefit from both the use of a polarizable force field and enhanced sampling. PMID:24634618
A coarse-grain force field for RDX: Density dependent and energy conserving
NASA Astrophysics Data System (ADS)
Moore, Joshua D.; Barnes, Brian C.; Izvekov, Sergei; Lísal, Martin; Sellers, Michael S.; Taylor, DeCarlos E.; Brennan, John K.
2016-03-01
We describe the development of a density-dependent transferable coarse-grain model of crystalline hexahydro-1,3,5-trinitro-s-triazine (RDX) that can be used with the energy conserving dissipative particle dynamics method. The model is an extension of a recently reported one-site model of RDX that was developed by using a force-matching method. The density-dependent forces in that original model are provided through an interpolation scheme that poorly conserves energy. The development of the new model presented in this work first involved a multi-objective procedure to improve the structural and thermodynamic properties of the previous model, followed by the inclusion of the density dependency via a conservative form of the force field that conserves energy. The new model accurately predicts the density, structure, pressure-volume isotherm, bulk modulus, and elastic constants of the RDX crystal at ambient pressure and exhibits transferability to a liquid phase at melt conditions.
Optical field and attractive force at the subwavelength slit.
Shapiro, David; Nies, Daniel; Belai, Oleg; Wurm, Matthias; Nesterov, Vladimir
2016-07-11
In recent works, a novel light-induced attractive force was predicted between two metal plates. This force arises by the interaction of surface plasmons which are excited at the metal when a transverse magnetic mode propagates through a subwavelength slit between two metal bodies. In this paper, the analytical and numerical calculations of this magnetic field are presented for the perfect metal and for gold. The amplitude and the phase transient curves between the known limiting cases of narrow and wide slits compared to the wavelength are found. The curve is shown to oscillate due to the emergence of new waveguide modes. The analytic solution for the perfect metal is in agreement with the computation for gold by means of the finite element method. The simple asymptotic formula for the light-induced attractive force is found in the limit of a narrow slit. PMID:27410865
Development of Field Excavator with Embedded Force Measurement
NASA Technical Reports Server (NTRS)
Johnson, K.; Creager, C.; Izadnegahdar, A.; Bauman, S.; Gallo, C.; Abel, P.
2012-01-01
A semi-intelligent excavation mechanism was developed for use with the NASA-built Centaur 2 rover prototype. The excavator features a continuously rotatable large bucket supported between two parallel arms, both of which share a single pivot axis near the excavator base attached to the rover. The excavator is designed to simulate the collection of regolith, such as on the Moon, and to dump the collected soil into a hopper up to one meter tall for processing to extract oxygen. Because the vehicle can be autonomous and the terrain is generally unknown, there is risk of damaging equipment or using excessive power when attempting to extract soil from dense or rocky terrain. To minimize these risks, it is critical for the rover to sense the digging forces and adjust accordingly. It is also important to understand the digging capabilities and limitations of the excavator. This paper discusses the implementation of multiple strain gages as an embedded force measurement system in the excavator's arms. These strain gages can accurately measure and resolve multi-axial forces on the excavator. In order to validate these sensors and characterize the load capabilities, a series of controlled excavation tests were performed at Glenn Research Center with the excavator at various depths and cut angles while supported by a six axis load cell. The results of these tests are both compared to a force estimation model and used for calibration of the embedded strain gages. In addition, excavation forces generated using two different types of bucket edge (straight vs. with teeth) were compared.
An assessment of the anharmonic force fields of the halogen cyanides
NASA Astrophysics Data System (ADS)
Lacy, M.
Anharmonic force fields in the ϱ representation (ϱ = Я/ r) have been calculated using the matrix diagonalization technique. A simple model force field and comparison method are used to assess the force fields currently available for these molecules. From this analysis, the potential surfaces defined by the present force fields are estimated to be within 2-5% of the actual molecular potential.
NASA Astrophysics Data System (ADS)
Christophe, B.; Lebat, V.; Foulon, B.; Liorzou, F.; Perrot, E.; Boulanger, D.; Hardy, E.
2014-12-01
ONERA has developed since several years the most accurate accelerometers for the geodesy mission. The accelerometers are still operational in the GRACE mission. Their successors for the GRACE-FO mission are under manufacturing and will fly in 2017. Finally, the GOCE mission has proved the interest of gradiometer for a direct measurement of the gravity field.Now, ONERA proposes a new design of accelerometer, MicroSTAR, for interplanetary mission. It inherits of the same technology but with reduced mass and consumption. It has been proposed in several missions towards outer planets in order to test the deviation to the relativity general over large distance to the sun (with the addition of a bias rejection system). But the same instrument could be interesting to improve our knowledge of the planetary gravitational potential field, allowing a better understanding of the planet interior composition. The success of using accelerometer for geodesy mission could be imported in the planetary science.The paper will present the accuracy achievable on the gravity potential field according to different accelerometer configurations (one accelerometer, one gradiometer arm or a complete 3-axis gradiometer). Then, the instrument will be described and the integration of the instrument inside an interplanetary probe will be evoked.
NASA Astrophysics Data System (ADS)
Jenkins, Micah; Gaylord, Thomas K.
2015-03-01
Most quantitative phase microscopy methods require the use of custom-built or modified microscopic configurations which are not typically available to most bio/pathologists. There are, however, phase retrieval algorithms which utilize defocused bright-field images as input data and are therefore implementable in existing laboratory environments. Among these, deterministic methods such as those based on inverting the transport-of-intensity equation (TIE) or a phase contrast transfer function (PCTF) are particularly attractive due to their compatibility with Köhler illuminated systems and numerical simplicity. Recently, a new method has been proposed, called multi-filter phase imaging with partially coherent light (MFPI-PC), which alleviates the inherent noise/resolution trade-off in solving the TIE by utilizing a large number of defocused bright-field images spaced equally about the focal plane. Despite greatly improving the state-ofthe- art, the method has many shortcomings including the impracticality of high-speed acquisition, inefficient sampling, and attenuated response at high frequencies due to aperture effects. In this report, we present a new method, called bright-field quantitative phase microscopy (BFQPM), which efficiently utilizes a small number of defocused bright-field images and recovers frequencies out to the partially coherent diffraction limit. The method is based on a noiseminimized inversion of a PCTF derived for each finite defocus distance. We present simulation results which indicate nanoscale optical path length sensitivity and improved performance over MFPI-PC. We also provide experimental results imaging live bovine mesenchymal stem cells at sub-second temporal resolution. In all, BFQPM enables fast and accurate phase imaging with unprecedented spatial resolution using widely available bright-field microscopy hardware.
Mapping the force field of a hydrogen-bonded assembly
Sweetman, A. M.; Jarvis, S. P.; Sang, Hongqian; Lekkas, I.; Rahe, P.; Wang, Yu; Wang, Jianbo; Champness, N.R.; Kantorovich, L.; Moriarty, P.
2014-01-01
Hydrogen bonding underpins the properties of a vast array of systems spanning a wide variety of scientific fields. From the elegance of base pair interactions in DNA to the symmetry of extended supramolecular assemblies, hydrogen bonds play an essential role in directing intermolecular forces. Yet fundamental aspects of the hydrogen bond continue to be vigorously debated. Here we use dynamic force microscopy (DFM) to quantitatively map the tip-sample force field for naphthalene tetracarboxylic diimide molecules hydrogen-bonded in two-dimensional assemblies. A comparison of experimental images and force spectra with their simulated counterparts shows that intermolecular contrast arises from repulsive tip-sample interactions whose interpretation can be aided via an examination of charge density depletion across the molecular system. Interpreting DFM images of hydrogen-bonded systems therefore necessitates detailed consideration of the coupled tip-molecule system: analyses based on intermolecular charge density in the absence of the tip fail to capture the essential physical chemistry underpinning the imaging mechanism. PMID:24875276
NASA Astrophysics Data System (ADS)
Sabatino, D. R.; Praisner, T. J.; Smith, C. R.
1998-11-01
The color change of thermochromic liquid crystals with temperature can be effectively utilized as full-field surface temperature sensors to investigate the fundamental structure of wall turbulence. In order to accurately quantify turbulent heat transfer behavior, a new technique has been developed for the calibration of wide-band micro-encapsulated thermochromic liquid crystals. Lighting/viewing arrangements are described and evaluated for ease of implementation and accuracy of the displayed color. This new technique employs images recorded in-situ with the test surface systematically exposed to a series of uniform temperature conditions spanning the bandwidth of the liquid crystals. This sequence of images is used to generate point-wise color/temperature calibration curves for the entire surface. Experimental results will be presented illustrating the application of the technique for assessment of spatial/temporal surface heat transfer behavior due to selected turbulent flows in a water channel
Efficient parametrization of complex molecule-surface force fields.
Gao, David Z; Federici Canova, Filippo; Watkins, Matthew B; Shluger, Alexander L
2015-06-15
We present an efficient scheme for parametrizing complex molecule-surface force fields from ab initio data. The cost of producing a sufficient fitting library is mitigated using a 2D periodic embedded slab model made possible by the quantum mechanics/molecular mechanics scheme in CP2K. These results were then used in conjunction with genetic algorithm (GA) methods to optimize the large parameter sets needed to describe such systems. The derived potentials are able to well reproduce adsorption geometries and adsorption energies calculated using density functional theory. Finally, we discuss the challenges in creating a sufficient fitting library, determining whether or not the GA optimization has completed, and the transferability of such force fields to similar molecules. PMID:25891018
A Maximum-Likelihood Approach to Force-Field Calibration.
Zaborowski, Bartłomiej; Jagieła, Dawid; Czaplewski, Cezary; Hałabis, Anna; Lewandowska, Agnieszka; Żmudzińska, Wioletta; Ołdziej, Stanisław; Karczyńska, Agnieszka; Omieczynski, Christian; Wirecki, Tomasz; Liwo, Adam
2015-09-28
A new approach to the calibration of the force fields is proposed, in which the force-field parameters are obtained by maximum-likelihood fitting of the calculated conformational ensembles to the experimental ensembles of training system(s). The maximum-likelihood function is composed of logarithms of the Boltzmann probabilities of the experimental conformations, calculated with the current energy function. Because the theoretical distribution is given in the form of the simulated conformations only, the contributions from all of the simulated conformations, with Gaussian weights in the distances from a given experimental conformation, are added to give the contribution to the target function from this conformation. In contrast to earlier methods for force-field calibration, the approach does not suffer from the arbitrariness of dividing the decoy set into native-like and non-native structures; however, if such a division is made instead of using Gaussian weights, application of the maximum-likelihood method results in the well-known energy-gap maximization. The computational procedure consists of cycles of decoy generation and maximum-likelihood-function optimization, which are iterated until convergence is reached. The method was tested with Gaussian distributions and then applied to the physics-based coarse-grained UNRES force field for proteins. The NMR structures of the tryptophan cage, a small α-helical protein, determined at three temperatures (T = 280, 305, and 313 K) by Hałabis et al. ( J. Phys. Chem. B 2012 , 116 , 6898 - 6907 ), were used. Multiplexed replica-exchange molecular dynamics was used to generate the decoys. The iterative procedure exhibited steady convergence. Three variants of optimization were tried: optimization of the energy-term weights alone and use of the experimental ensemble of the folded protein only at T = 280 K (run 1); optimization of the energy-term weights and use of experimental ensembles at all three temperatures (run 2
Force-field parameters for beryllium complexes in amorphous layers.
Emelyanova, Svetlana; Chashchikhin, Vladimir; Bagaturyants, Alexander
2016-09-01
Unknown force-field parameters for metal organic beryllium complexes used in emitting and electron transporting layers of OLED structures are determined. These parameters can be used for the predictive atomistic simulations of the structure and properties of amorphous organic layers containing beryllium complexes. The parameters are found for the AMBER force field using a relaxed scan procedure and quantum-mechanical DFT calculations of potential energy curves for specific internal (angular) coordinates in a series of three Be complexes (Bebq2; Be(4-mpp)2; Bepp2). The obtained parameters are verified in calculations of some molecular and crystal structures available from either quantum-mechanical DFT calculations or experimental data. Graphical Abstract Beryllium complexes in amorphous layersᅟ. PMID:27550375
Tuning the mass of chameleon fields in Casimir force experiments.
Brax, Ph; van de Bruck, C; Davis, A C; Shaw, D J; Iannuzzi, D
2010-06-18
We have calculated the chameleon pressure between two parallel plates in the presence of an intervening medium that affects the mass of the chameleon field. As intuitively expected, the gas in the gap weakens the chameleon interaction mechanism with a screening effect that increases with the plate separation and with the density of the intervening medium. This phenomenon might open up new directions in the search of chameleon particles with future long-range Casimir force experiments. PMID:20867290
Quantum mechanical force field for water with explicit electronic polarization
Han, Jaebeom; Mazack, Michael J. M.; Zhang, Peng; Truhlar, Donald G.; Gao, Jiali
2013-08-07
A quantum mechanical force field (QMFF) for water is described. Unlike traditional approaches that use quantum mechanical results and experimental data to parameterize empirical potential energy functions, the present QMFF uses a quantum mechanical framework to represent intramolecular and intermolecular interactions in an entire condensed-phase system. In particular, the internal energy terms used in molecular mechanics are replaced by a quantum mechanical formalism that naturally includes electronic polarization due to intermolecular interactions and its effects on the force constants of the intramolecular force field. As a quantum mechanical force field, both intermolecular interactions and the Hamiltonian describing the individual molecular fragments can be parameterized to strive for accuracy and computational efficiency. In this work, we introduce a polarizable molecular orbital model Hamiltonian for water and for oxygen- and hydrogen-containing compounds, whereas the electrostatic potential responsible for intermolecular interactions in the liquid and in solution is modeled by a three-point charge representation that realistically reproduces the total molecular dipole moment and the local hybridization contributions. The present QMFF for water, which is called the XP3P (explicit polarization with three-point-charge potential) model, is suitable for modeling both gas-phase clusters and liquid water. The paper demonstrates the performance of the XP3P model for water and proton clusters and the properties of the pure liquid from about 900 × 10{sup 6} self-consistent-field calculations on a periodic system consisting of 267 water molecules. The unusual dipole derivative behavior of water, which is incorrectly modeled in molecular mechanics, is naturally reproduced as a result of an electronic structural treatment of chemical bonding by XP3P. We anticipate that the XP3P model will be useful for studying proton transport in solution and solid phases as well as
Quantum mechanical force field for water with explicit electronic polarization
Han, Jaebeom; Mazack, Michael J. M.; Zhang, Peng; Truhlar, Donald G.; Gao, Jiali
2013-01-01
A quantum mechanical force field (QMFF) for water is described. Unlike traditional approaches that use quantum mechanical results and experimental data to parameterize empirical potential energy functions, the present QMFF uses a quantum mechanical framework to represent intramolecular and intermolecular interactions in an entire condensed-phase system. In particular, the internal energy terms used in molecular mechanics are replaced by a quantum mechanical formalism that naturally includes electronic polarization due to intermolecular interactions and its effects on the force constants of the intramolecular force field. As a quantum mechanical force field, both intermolecular interactions and the Hamiltonian describing the individual molecular fragments can be parameterized to strive for accuracy and computational efficiency. In this work, we introduce a polarizable molecular orbital model Hamiltonian for water and for oxygen- and hydrogen-containing compounds, whereas the electrostatic potential responsible for intermolecular interactions in the liquid and in solution is modeled by a three-point charge representation that realistically reproduces the total molecular dipole moment and the local hybridization contributions. The present QMFF for water, which is called the XP3P (explicit polarization with three-point-charge potential) model, is suitable for modeling both gas-phase clusters and liquid water. The paper demonstrates the performance of the XP3P model for water and proton clusters and the properties of the pure liquid from about 900 × 106 self-consistent-field calculations on a periodic system consisting of 267 water molecules. The unusual dipole derivative behavior of water, which is incorrectly modeled in molecular mechanics, is naturally reproduced as a result of an electronic structural treatment of chemical bonding by XP3P. We anticipate that the XP3P model will be useful for studying proton transport in solution and solid phases as well as across
Force-free magnetic fields - Generating functions and footpoint displacements
NASA Technical Reports Server (NTRS)
Wolfson, Richard; Verma, Ritu
1991-01-01
This paper presents analytic and numerical calculations that explore equilibrium sequences of bipolar force-free magnetic fields in relation to displacments of their magnetic footpoints. It is shown that the appearance of magnetic islands - sometimes interpreted as marking the loss of equilibrium in models of the solar atmosphere - is likely associated only with physically unrealistic footpoint displacements such as infinite separation or 'tearing' of the model photosphere. The work suggests that the loss of equilibrium in bipolar configurations, sometimes proposed as a mechanism for eruptive solar events, probably requires either fully three-dimensional field configurations or nonzero plasma pressure. The results apply only to fields that are strictly bipolar, and do not rule out equilibrium loss in more complex structures such as quadrupolar fields.
Design and optimization of force-reduced high field magnets
NASA Astrophysics Data System (ADS)
Rembeczki, Szabolcs
High field magnets have many important applications in different areas of research, in the power industry and also for military purposes. For example, high field magnets are particularly useful in: material sciences, high energy physics, plasma physics (as fusion magnets), high power applications (as energy storage devices), and space applications (in propulsion systems). One of the main issues with high-field magnets is the presence of very large electromagnetic stresses that must be counteracted and therefore require heavy support structures. In superconducting magnets, the problems caused by Lorentz forces are further complicated by the fact that superconductors for high field applications are pressure sensitive. The current carrying capacity is greatly reduced under stress and strain (especially in the case of Nb 3Sn and the new high temperature superconductors) so the reduction of the acting forces is of even greater importance. Different force-reduced magnet concepts have been studied in the past, both numerical and analytical methods have been used to solve this problem. The developed concepts are based on such complex winding geometries that the realization and manufacturing of such coils is extremely difficult and these concepts are mainly of theoretical interest. In the presented research, a novel concept for force-reduced magnets has been developed and analyzed which is easy to realize and therefore is of practical interest. The analysis has been performed with a new methodology, which does not require the time consuming finite element calculations. The developed computer models describe the 3-dimensional winding configuration by sets of filaments (filamentary approximation). This approach is much faster than finite element analysis and therefore allows rapid optimization of concepts. The method has been extensively tested on geometries of force-reduced solenoids where even analytical solutions exist. As a further cross check, the developed computer
Representation of Ion–Protein Interactions Using the Drude Polarizable Force-Field
2016-01-01
Small metal ions play critical roles in numerous biological processes. Of particular interest is how metalloenzymes are allosterically regulated by the binding of specific ions. Understanding how ion binding affects these biological processes requires atomic models that accurately treat the microscopic interactions with the protein ligands. Theoretical approaches at different levels of sophistication can contribute to a deeper understanding of these systems, although computational models must strike a balance between accuracy and efficiency in order to enable long molecular dynamics simulations. In this study, we present a systematic effort to optimize the parameters of a polarizable force field based on classical Drude oscillators to accurately represent the interactions between ions (K+, Na+, Ca2+, and Cl–) and coordinating amino-acid residues for a set of 30 biologically important proteins. By combining ab initio calculations and experimental thermodynamic data, we derive a polarizable force field that is consistent with a wide range of properties, including the geometries and interaction energies of gas-phase ion/protein-like model compound clusters, and the experimental solvation free-energies of the cations in liquids. The resulting models display significant improvements relative to the fixed-atomic-charge additive CHARMM C36 force field, particularly in their ability to reproduce the many-body electrostatic nonadditivity effects estimated from ab initio calculations. The analysis clarifies the fundamental limitations of the pairwise additivity assumption inherent in classical fixed-charge force fields, and shows its dramatic failures in the case of Ca2+ binding sites. These optimized polarizable models, amenable to computationally efficient large-scale MD simulations, set a firm foundation and offer a powerful avenue to study the roles of the ions in soluble and membrane transport proteins. PMID:25578354
Representation of Ion-Protein Interactions Using the Drude Polarizable Force-Field.
Li, Hui; Ngo, Van; Da Silva, Mauricio Chagas; Salahub, Dennis R; Callahan, Karen; Roux, Benoît; Noskov, Sergei Yu
2015-07-23
Small metal ions play critical roles in numerous biological processes. Of particular interest is how metalloenzymes are allosterically regulated by the binding of specific ions. Understanding how ion binding affects these biological processes requires atomic models that accurately treat the microscopic interactions with the protein ligands. Theoretical approaches at different levels of sophistication can contribute to a deeper understanding of these systems, although computational models must strike a balance between accuracy and efficiency in order to enable long molecular dynamics simulations. In this study, we present a systematic effort to optimize the parameters of a polarizable force field based on classical Drude oscillators to accurately represent the interactions between ions (K(+), Na(+), Ca(2+), and Cl(-)) and coordinating amino-acid residues for a set of 30 biologically important proteins. By combining ab initio calculations and experimental thermodynamic data, we derive a polarizable force field that is consistent with a wide range of properties, including the geometries and interaction energies of gas-phase ion/protein-like model compound clusters, and the experimental solvation free-energies of the cations in liquids. The resulting models display significant improvements relative to the fixed-atomic-charge additive CHARMM C36 force field, particularly in their ability to reproduce the many-body electrostatic nonadditivity effects estimated from ab initio calculations. The analysis clarifies the fundamental limitations of the pairwise additivity assumption inherent in classical fixed-charge force fields, and shows its dramatic failures in the case of Ca(2+) binding sites. These optimized polarizable models, amenable to computationally efficient large-scale MD simulations, set a firm foundation and offer a powerful avenue to study the roles of the ions in soluble and membrane transport proteins. PMID:25578354
Zhu, Xiao; MacKerell, Alexander D.
2010-01-01
Condensed-phase computational studies of molecules using molecular mechanics approaches require the use of force fields to describe the energetics of the systems as a function of structure. The advantage of polarizable force fields over non-polarizable (or additive) models lies in their ability to vary their electronic distribution as a function of the environment. Towards development of a polarizable force field for biological molecules, parameters for a series of sulfur-containing molecules are presented. Parameter optimization was performed to reproduce quantum mechanical and experimental data for gas phase properties including geometries, conformational energies, vibrational spectra, and dipole moments as well as for condensed phase properties such as heats of vaporization, molecular volumes, and free energies of hydration. Compounds in the training set include methanethiol, ethanethiol, propanethiol, ethyl methyl sulfide, and dimethyl disulfide. The molecular volumes and heats of vaporization are in good accordance with experimental values, with the polarizable model performing better than the CHARMM22 non-polarizable force field. Improvements with the polarizable model were also obtained for molecular dipole moments and in the treatment of intermolecular interactions as a function of orientation, in part due to the presence of lone pairs and anisotropic atomic polarizability on the sulfur atoms. Significant advantage of the polarizable model was reflected in calculation of the dielectric constants, a property that CHARMM22 systematically underestimates. The ability of this polarizable model to accurately describe a range of gas and condensed phase properties paves the way for more accurate simulation studies of sulfur-containing molecules including cysteine and methionine residues in proteins. PMID:20575015
Levitation forces of a bulk YBCO superconductor in gradient varying magnetic fields
NASA Astrophysics Data System (ADS)
Jiang, J.; Gong, Y. M.; Wang, G.; Zhou, D. J.; Zhao, L. F.; Zhang, Y.; Zhao, Y.
2015-09-01
The levitation forces of a bulk YBCO superconductor in gradient varying high and low magnetic fields generated from a superconducting magnet were investigated. The magnetic field intensity of the superconducting magnet was measured when the exciting current was 90 A. The magnetic field gradient and magnetic force field were both calculated. The YBCO bulk was cooled by liquid nitrogen in field-cooling (FC) and zero-field-cooling (ZFC) condition. The results showed that the levitation forces increased with increasing the magnetic field intensity. Moreover, the levitation forces were more dependent on magnetic field gradient and magnetic force field than magnetic field intensity.
The Amber ff99 Force Field Predicts Relative Free Energy Changes for RNA Helix Formation
Spasic, Aleksandar; Serafini, John; Mathews, David H.
2012-01-01
The ability of the Amber ff99 force field to predict relative free energies of RNA helix formation was investigated. The test systems were three hexaloop RNA hairpins with identical loops and varying stems. The potential of mean force of stretching the hairpins from the native state to an extended conformation was calculated with umbrella sampling. Because the hairpins have identical loop sequence, the differences in free energy changes are only from the stem composition. The Amber ff99 force field was able to correctly predict the order of stabilities of the hairpins, although the magnitude of the free energy change is larger than that determined by optical melting experiments. The two measurements cannot be compared directly because the unfolded state in the optical melting experiments is a random coil, while the end state in the umbrella sampling simulations was an elongated chain. The calculations can be compared to reference data by using a thermodynamic cycle. By applying the thermodynamic cycle to the transitions between the hairpins using simulations and nearest neighbor data, agreement was found to be within the sampling error of simulations, thus demonstrating that ff99 force field is able to accurately predict relative free energies of RNA helix formation. PMID:23112748
Optical forces on a Rayleigh dielectric particle in a patterned near-field landscape
NASA Astrophysics Data System (ADS)
Zelenina, Anna; Quidant, Romain; Petrov, Dmitri; Badenes, Goncal
2005-07-01
The trapping of micro-objects by optical radiation forces, so-called optical tweezers, has become widely used in physical, chemical and biological experiments where accurate and non-invasive manipulation is required. Recent advances in beam shaping render it possible for instance to rotate or to dynamically manipulate independently several elements. Today, one of the remaining challenges of conventional optical tweezers is the direct manipulation of systems with sizes belonging to the sub-wavelength or Rayleigh regime. Indeed, the diffraction limit prevents in that case from achieving a commensurable trapping volume and thus does not allow for minimizing the fluctuations in position of the trapped object due to its strong Brownian motion. In order to overcome this limitation, it has been proposed to use evanescent fields instead of the usual propagating fields. Recent advances in optics of noble metal nano-structures have recently provided new configurations to achieve nano-optical tweezers. Especially, tightly localized modes resulting from the coupling between resonant noble metal nanostructures may offer the gradient forces able to trap and manipulate Rayleigh objects. In this work, we calculate the radiation forces exerted on a nanometric dielectric sphere when exposed to a patterned optical near-field landscape at an interface decorated with resonant gold nanostructures. By comparing their magnitude with other forces that affect the movement of the particle, we discuss the practical ability of our configuration for multiple parallel optical manipulation.
Secondary Structure of Rat and Human Amylin across Force Fields
Hoffmann, Kyle Quynn; McGovern, Michael; Chiu, Chi-cheng; de Pablo, Juan J.
2015-01-01
The aggregation of human amylin has been strongly implicated in the progression of Type II diabetes. This 37-residue peptide forms a variety of secondary structures, including random coils, α-helices, and β-hairpins. The balance between these structures depends on the chemical environment, making amylin an ideal candidate to examine inherent biases in force fields. Rat amylin differs from human amylin by only 6 residues; however, it does not form fibrils. Therefore it provides a useful complement to human amylin in studies of the key events along the aggregation pathway. In this work, the free energy of rat and human amylin was determined as a function of α-helix and β-hairpin content for the Gromos96 53a6, OPLS-AA/L, CHARMM22/CMAP, CHARMM22*, Amberff99sb*-ILDN, and Amberff03w force fields using advanced sampling techniques, specifically bias exchange metadynamics. This work represents a first systematic attempt to evaluate the conformations and the corresponding free energy of a large, clinically relevant disordered peptide in solution across force fields. The NMR chemical shifts of rIAPP were calculated for each of the force fields using their respective free energy maps, allowing us to quantitatively assess their predictions. We show that the predicted distribution of secondary structures is sensitive to the choice of force-field: Gromos53a6 is biased towards β-hairpins, while CHARMM22/CMAP predicts structures that are overly α-helical. OPLS-AA/L favors disordered structures. Amberff99sb*-ILDN, AmberFF03w and CHARMM22* provide the balance between secondary structures that is most consistent with available experimental data. In contrast to previous reports, our findings suggest that the equilibrium conformations of human and rat amylin are remarkably similar, but that subtle differences arise in transient alpha-helical and beta-strand containing structures that the human peptide can more readily adopt. We hypothesize that these transient states enable
Secondary structure of rat and human amylin across force fields
Hoffmann, Kyle Quynn; McGovern, Michael; Chiu, Chi -cheng; de Pablo, Juan J.; Paci, Emanuele
2015-07-29
The aggregation of human amylin has been strongly implicated in the progression of Type II diabetes. This 37-residue peptide forms a variety of secondary structures, including random coils, α-helices, and β-hairpins. The balance between these structures depends on the chemical environment, making amylin an ideal candidate to examine inherent biases in force fields. Rat amylin differs from human amylin by only 6 residues; however, it does not form fibrils. Therefore it provides a useful complement to human amylin in studies of the key events along the aggregation pathway. In this work, the free energy of rat and human amylin wasmore » determined as a function of α-helix and β-hairpin content for the Gromos96 53a6, OPLS-AA/L, CHARMM22/CMAP, CHARMM22*, Amberff99sb*-ILDN, and Amberff03w force fields using advanced sampling techniques, specifically bias exchange metadynamics. This work represents a first systematic attempt to evaluate the conformations and the corresponding free energy of a large, clinically relevant disordered peptide in solution across force fields. The NMR chemical shifts of rIAPP were calculated for each of the force fields using their respective free energy maps, allowing us to quantitatively assess their predictions. We show that the predicted distribution of secondary structures is sensitive to the choice of force-field: Gromos53a6 is biased towards β-hairpins, while CHARMM22/CMAP predicts structures that are overly α-helical. OPLS-AA/L favors disordered structures. Amberff99sb*-ILDN, AmberFF03w and CHARMM22* provide the balance between secondary structures that is most consistent with available experimental data. In contrast to previous reports, our findings suggest that the equilibrium conformations of human and rat amylin are remarkably similar, but that subtle differences arise in transient alpha-helical and beta-strand containing structures that the human peptide can more readily adopt. We hypothesize that these transient states
Secondary structure of rat and human amylin across force fields
Hoffmann, Kyle Quynn; McGovern, Michael; Chiu, Chi -cheng; de Pablo, Juan J.; Paci, Emanuele
2015-07-29
The aggregation of human amylin has been strongly implicated in the progression of Type II diabetes. This 37-residue peptide forms a variety of secondary structures, including random coils, α-helices, and β-hairpins. The balance between these structures depends on the chemical environment, making amylin an ideal candidate to examine inherent biases in force fields. Rat amylin differs from human amylin by only 6 residues; however, it does not form fibrils. Therefore it provides a useful complement to human amylin in studies of the key events along the aggregation pathway. In this work, the free energy of rat and human amylin was determined as a function of α-helix and β-hairpin content for the Gromos96 53a6, OPLS-AA/L, CHARMM22/CMAP, CHARMM22*, Amberff99sb*-ILDN, and Amberff03w force fields using advanced sampling techniques, specifically bias exchange metadynamics. This work represents a first systematic attempt to evaluate the conformations and the corresponding free energy of a large, clinically relevant disordered peptide in solution across force fields. The NMR chemical shifts of rIAPP were calculated for each of the force fields using their respective free energy maps, allowing us to quantitatively assess their predictions. We show that the predicted distribution of secondary structures is sensitive to the choice of force-field: Gromos53a6 is biased towards β-hairpins, while CHARMM22/CMAP predicts structures that are overly α-helical. OPLS-AA/L favors disordered structures. Amberff99sb*-ILDN, AmberFF03w and CHARMM22* provide the balance between secondary structures that is most consistent with available experimental data. In contrast to previous reports, our findings suggest that the equilibrium conformations of human and rat amylin are remarkably similar, but that subtle differences arise in transient alpha-helical and beta-strand containing structures that the human peptide can more readily adopt. We hypothesize that these transient states enable
Secondary Structure of Rat and Human Amylin across Force Fields.
Hoffmann, Kyle Quynn; McGovern, Michael; Chiu, Chi-Cheng; de Pablo, Juan J
2015-01-01
The aggregation of human amylin has been strongly implicated in the progression of Type II diabetes. This 37-residue peptide forms a variety of secondary structures, including random coils, α-helices, and β-hairpins. The balance between these structures depends on the chemical environment, making amylin an ideal candidate to examine inherent biases in force fields. Rat amylin differs from human amylin by only 6 residues; however, it does not form fibrils. Therefore it provides a useful complement to human amylin in studies of the key events along the aggregation pathway. In this work, the free energy of rat and human amylin was determined as a function of α-helix and β-hairpin content for the Gromos96 53a6, OPLS-AA/L, CHARMM22/CMAP, CHARMM22*, Amberff99sb*-ILDN, and Amberff03w force fields using advanced sampling techniques, specifically bias exchange metadynamics. This work represents a first systematic attempt to evaluate the conformations and the corresponding free energy of a large, clinically relevant disordered peptide in solution across force fields. The NMR chemical shifts of rIAPP were calculated for each of the force fields using their respective free energy maps, allowing us to quantitatively assess their predictions. We show that the predicted distribution of secondary structures is sensitive to the choice of force-field: Gromos53a6 is biased towards β-hairpins, while CHARMM22/CMAP predicts structures that are overly α-helical. OPLS-AA/L favors disordered structures. Amberff99sb*-ILDN, AmberFF03w and CHARMM22* provide the balance between secondary structures that is most consistent with available experimental data. In contrast to previous reports, our findings suggest that the equilibrium conformations of human and rat amylin are remarkably similar, but that subtle differences arise in transient alpha-helical and beta-strand containing structures that the human peptide can more readily adopt. We hypothesize that these transient states enable
Heinz, Hendrik; Lin, Tzu-Jen; Mishra, Ratan Kishore; Emami, Fateme S
2013-02-12
The complexity of the molecular recognition and assembly of biotic-abiotic interfaces on a scale of 1 to 1000 nm can be understood more effectively using simulation tools along with laboratory instrumentation. We discuss the current capabilities and limitations of atomistic force fields and explain a strategy to obtain dependable parameters for inorganic compounds that has been developed and tested over the past decade. Parameter developments include several silicates, aluminates, metals, oxides, sulfates, and apatites that are summarized in what we call the INTERFACE force field. The INTERFACE force field operates as an extension of common harmonic force fields (PCFF, COMPASS, CHARMM, AMBER, GROMACS, and OPLS-AA) by employing the same functional form and combination rules to enable simulations of inorganic-organic and inorganic-biomolecular interfaces. The parametrization builds on an in-depth understanding of physical-chemical properties on the atomic scale to assign each parameter, especially atomic charges and van der Waals constants, as well as on the validation of macroscale physical-chemical properties for each compound in comparison to measurements. The approach eliminates large discrepancies between computed and measured bulk and surface properties of up to 2 orders of magnitude using other parametrization protocols and increases the transferability of the parameters by introducing thermodynamic consistency. As a result, a wide range of properties can be computed in quantitative agreement with experiment, including densities, surface energies, solid-water interface tensions, anisotropies of interfacial energies of different crystal facets, adsorption energies of biomolecules, and thermal and mechanical properties. Applications include insight into the assembly of inorganic-organic multiphase materials, the recognition of inorganic facets by biomolecules, growth and shape preferences of nanocrystals and nanoparticles, as well as thermal transitions and
Force field measurements within the exclusion zone of water.
Chen, Chi-Shuo; Chung, Wei-Ju; Hsu, Ian C; Wu, Chien-Ming; Chin, Wei-Chun
2012-01-01
Water molecules play critical roles in many biological functions, such as protein dynamics, enzymatic activities, and cellular responses. Previous nuclear magnetic resonance and neutron scattering studies have shown that water molecules bind to specific sites on surfaces and form localized clusters. However, most current experimental techniques cannot measure dynamic behaviors of ordered water molecules on cell-size (10 μm) scale. Recently, the long-distance effect of structured water has been demonstrated by Pollack and his colleagues. Namely, there is a structured water layer near the hydrophilic surface that can exclude solutes (Zheng et al, Adv Colloid Interface Sci 127:19-27, 2006; Pollack 2006, Adv Colloid Interface Sci 103:173-196, 2003). The repelling forces of water clusters inside this exclusion region are investigated in this study. With a laser tweezers system, we found the existence of an unexpected force fields inside the solute-free exclusion zone near a Nafion surface. Our results suggest that the water clusters could transduce mechanical signals on the micrometer range within the exclusion zone. This unexpected inhomogeneous force field near the hydrophilic surface would provide a new insight into cellular activities, leading to a potential new physical chemistry mechanism for cell biology. PMID:23277674
Coarse-grained force field; general folding theory
Liwo, Adam; He, Yi; Scheraga, Harold A.
2012-01-01
We review the coarse-grained UNited RESidue (UNRES) force field for the simulations of protein structure and dynamics, which is being developed in our laboratory over the last several years. UNRES is a physics-based force field, the prototype of which is defined as a potential of mean force of polypeptide chains in water, where all the degrees of freedom except the coordinates of α-carbon atoms and side-chain centers have been integrated out. We describe the initial implementation of UNRES to protein-structure prediction formulated as a search for the global minimum of the potential-energy function and its subsequent molecular dynamics and extensions of molecular-dynamics implementation, which enabled us to study protein-folding pathways and thermodynamics, as well as to reformulate the protein-structure prediction problem as a search for the conformational ensemble with the lowest free energy at temperatures below the folding-transition temperature. Applications of UNRES to study biological problems are also described. PMID:21643583
Interfacial Force Field Characterization in a Constrained Vapor Bubble Thermosyphon
NASA Technical Reports Server (NTRS)
DasGupta, Sunando; Plawsky, Joel L.; Wayner, Peter C., Jr.
1995-01-01
Isothermal profiles of the extended meniscus in a quartz cuvette were measured in the earth's gravitational field using an image-analyzing interferometer that is based on computer-enhanced video microscopy of the naturally occurring interference fringes. These profiles are a function of the stress field. Experimentally, the augmented Young-Laplace equation is an excellent model for the force field at the solid-liquid-vapor interfaces for heptane and pentane menisci on quartz and tetradecane on SFL6. The effects of refractive indices of the solid and liquid on the measurement techniques were demonstrated. Experimentally obtained values of the disjoining pressure and dispersion constants were compared to those predicted from the Dzyaloshinskii - Lifshitz - Pilaevskii theory for an ideal surface and reasonable agreements were obtained. A parameter introduced gives a quantitative measurement of the closeness of the system to equilibrium. The nonequilibrium behavior of this parameter is also presented
Sodium Chloride, NaCl/ϵ: New Force Field.
Fuentes-Azcatl, Raúl; Barbosa, Marcia C
2016-03-10
A new computational model for sodium chloride, the NaCl/ϵ, is proposed. The force field employed for the description of the NaCl is based on a set of radial particle-particle pair potentials involving Lennard-Jones (LJ) and Coulombic forces. The parametrization is obtained by fitting the density of the crystal and the density and the dielectric constant of the mixture of the salt with water at a diluted solution. Our model shows good agreement with the experimental values for the density and for the surface tension of the pure system, and for the density, the viscosity, the diffusion, and the dielectric constant for the mixture with water at various molal concentrations. The NaCl/ϵ together with the water TIP4P/ϵ models provide a good approximation for studying electrolyte solutions. PMID:26890321
Deformation field of the soft substrate induced by capillary force
NASA Astrophysics Data System (ADS)
Liu, J. L.; Nie, Z. X.; Jiang, W. G.
2009-05-01
Prediction on the deformation of a soft substrate induced by capillary force has been widely paid attention in the broad range of applications, such as metallurgy, material science, astronavigation, micro/nano-technology, etc., which is also a supplementary result to the classical Young's equation. We quantitatively analyzed the deformation of an elastic substrate under capillary force by means of the energy principle and the continuum mechanics method. The actual drop's morphology was investigated and was compared with that calculated based on the classical spherical shape assumption of the droplet. The displacement field of the substrate was obtained, especially, its singularity at the droplet edge was also discussed. The results are beneficial to engineering application and micro/nano-measurement.
Pairwise-additive force fields for selected aqueous monovalent ions from adaptive force matching
NASA Astrophysics Data System (ADS)
Li, Jicun; Wang, Feng
2015-11-01
Simple non-polarizable potentials were developed for Na+, K+, Cl-, and Br- using the adaptive force matching (AFM) method with ab initio MP2 method as reference. Our MP2-AFM force field predicts the solvation free energies of the four salts formed by the ions with an error of no more than 5%. Other properties such as the ion-water radial distribution functions, first solvation shell water tilt angle distributions, ion diffusion constants, concentration dependent diffusion constant of water, and concentration dependent surface tension of the solutions were calculated with this potential. Very good agreement was achieved for these properties. In particular, the diffusion constants of the ions are within 6% of experimental measurements. The model predicts bromide to be enriched at the interface in the 1.6M KBr solution but predicts the ion to be repelled for the surface at lower concentration.
NASA Astrophysics Data System (ADS)
MacCready, P.; Bryan, F.; Tseng, Y. H.; Whitney, M. M.
2014-12-01
The coastal ocean accounts for about half of the global fish harvest, but is poorly resolved in global climate models (a one-degree grid barely sees the continental shelf). Moreover, coastal ocean circulation is strongly modified by river freshwater sources, often coming from estuarine systems that are completely unresolved in the coarse grid. River freshwater input in CESM is added in a practical but ad hoc way, by imposing a surface salinity sink over a region of the ocean approximating the plume area of a given river. Here we present results from a series of model experiments using a high-resolution (1.5 km) ROMS model of the NE Pacific, including the Columbia River and the inland waters of Puget Sound. The base model does multi-year hindcasts using the best available sources of atmospheric (MM5/WRF), ocean (NCOM), river (USGS), and tidal forcing. It has been heavily validated against observations of all sorts, and performs well, so it is an ideal test bed for downscaling experiments. The model framework also does biogeochemistry, including oxygen, and carbon chemistry is being added to make forecasts of Ocean Acidification.This high-resolution ROMS model is systematically run in downscaling experiments for the year 2005 with combinations of CESM forcing (CAM, POP, and rivers) swapped in. Skill is calculated using observations. It is found that the runs with CESM forcing generally retain much of the skill of the base model. A compact metric of response to freshwater forcing is used, which is the mechanical energy required to destratify a shallow coastal volume. This, along with the average temperature and salinity of the volume, are used to characterize and compare runs, including the original CESM-POP fields. Finally the model is run with projected CESM simulation forcing at the end of 21st century based on a set of RCP scenarios, and the compact metrics are used to quantify differences from 2005.
A new cation-exchange method for accurate field speciation of hexavalent chromium
Ball, J.W.; McCleskey, R.B.
2003-01-01
A new method for field speciation of Cr(VI) has been developed to meet present stringent regulatory standards and to overcome the limitations of existing methods. The method consists of passing a water sample through strong acid cation-exchange resin at the field site, where Cr(III) is retained while Cr(VI) passes into the effluent and is preserved for later determination. The method is simple, rapid, portable, and accurate, and makes use of readily available, inexpensive materials. Cr(VI) concentrations are determined later in the laboratory using any elemental analysis instrument sufficiently sensitive to measure the Cr(VI) concentrations of interest. The new method allows measurement of Cr(VI) concentrations as low as 0.05 ??g 1-1, storage of samples for at least several weeks prior to analysis, and use of readily available analytical instrumentation. Cr(VI) can be separated from Cr(III) between pH 2 and 11 at Cr(III)/Cr(VI) concentration ratios as high as 1000. The new method has demonstrated excellent comparability with two commonly used methods, the Hach Company direct colorimetric method and USEPA method 218.6. The new method is superior to the Hach direct colorimetric method owing to its relative sensitivity and simplicity. The new method is superior to USEPA method 218.6 in the presence of Fe(II) concentrations up to 1 mg 1-1 and Fe(III) concentrations up to 10 mg 1-1. Time stability of preserved samples is a significant advantage over the 24-h time constraint specified for USEPA method 218.6.
Continuum Polarizable Force Field within the Poisson-Boltzmann Framework
Tan, Yu-Hong; Tan, Chunhu; Wang, Junmei; Luo, Ray
2008-01-01
We have developed and tested a complete set of nonbonded parameters for a continuum polarizable force field. Our analysis shows that the new continuum polarizable model is consistent with B3LYP/cc-pVTZ in modeling electronic response upon variation of dielectric environment. Comparison with experiment also shows that the new continuum polarizable model is reasonable, with similar accuracy as B3LYP/cc-pVTZ in reproduction of dipole moments of selected organic molecules in the gas phase. We have further tested the validity to interchange the Amber van der Waals parameters between the explicit and continuum polarizable force fields with a series of dimers. It can be found that the continuum polarizable model agrees well with MP2/cc-pVTZ, with deviations in dimer binding energies less than 0.9 kcal/mol in the aqueous dielectric environment. Finally we have optimized atomic cavity radii with respect to experimental solvation free energies of 177 training molecules. To validate the optimized cavity radii, we have tested these parameters against 176 test molecules. It is found that the optimized PB atomic cavity radii transfer well from the training set to the test set, with an overall root-mean-squared deviation of 1.30 kcal/mol, unsigned average error of 1.07 kacl/mol, and correlation coefficient of 92% for all 353 molecules in both the training and test sets. Given the development documented here, the next natural step is the construction of a full protein/nucleic acid force field within the new continuum polarization framework. PMID:18507452
Beekhuizen, Johan; Kromhout, Hans; Bürgi, Alfred; Huss, Anke; Vermeulen, Roel
2015-01-01
The increase in mobile communication technology has led to concern about potential health effects of radio frequency electromagnetic fields (RF-EMFs) from mobile phone base stations. Different RF-EMF prediction models have been applied to assess population exposure to RF-EMF. Our study examines what input data are needed to accurately model RF-EMF, as detailed data are not always available for epidemiological studies. We used NISMap, a 3D radio wave propagation model, to test models with various levels of detail in building and antenna input data. The model outcomes were compared with outdoor measurements taken in Amsterdam, the Netherlands. Results showed good agreement between modelled and measured RF-EMF when 3D building data and basic antenna information (location, height, frequency and direction) were used: Spearman correlations were >0.6. Model performance was not sensitive to changes in building damping parameters. Antenna-specific information about down-tilt, type and output power did not significantly improve model performance compared with using average down-tilt and power values, or assuming one standard antenna type. We conclude that 3D radio wave propagation modelling is a feasible approach to predict outdoor RF-EMF levels for ranking exposure levels in epidemiological studies, when 3D building data and information on the antenna height, frequency, location and direction are available. PMID:24472756
Candel, A.; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Prudencio, E.; Schussman, G.; Uplenchwar, R.; Ko, K.; /SLAC
2009-06-19
Over the past years, SLAC's Advanced Computations Department (ACD), under SciDAC sponsorship, has developed a suite of 3D (2D) parallel higher-order finite element (FE) codes, T3P (T2P) and Pic3P (Pic2P), aimed at accurate, large-scale simulation of wakefields and particle-field interactions in radio-frequency (RF) cavities of complex shape. The codes are built on the FE infrastructure that supports SLAC's frequency domain codes, Omega3P and S3P, to utilize conformal tetrahedral (triangular)meshes, higher-order basis functions and quadratic geometry approximation. For time integration, they adopt an unconditionally stable implicit scheme. Pic3P (Pic2P) extends T3P (T2P) to treat charged-particle dynamics self-consistently using the PIC (particle-in-cell) approach, the first such implementation on a conformal, unstructured grid using Whitney basis functions. Examples from applications to the International Linear Collider (ILC), Positron Electron Project-II (PEP-II), Linac Coherent Light Source (LCLS) and other accelerators will be presented to compare the accuracy and computational efficiency of these codes versus their counterparts using structured grids.
The Anharmonic Force Field of BeH2 Revisited
NASA Technical Reports Server (NTRS)
Martin, Jan M. L.; Lee, Timothy J.
2003-01-01
The anharmonic force field of BeH2 has been calculated near the basis set and n-particle space limits. The computed antisymmetric stretch frequencies of BeH2 and BeD2 are in excellent agreement with recent high-resolution gas-phase measurements. The agreement between theory and experiment for the other spectroscopic constants is also excellent, except for omega(sub 3) and X(sub 33) for BeH2 and G(sub 22) for BeD2. It is concluded that further experimental work is needed in order to resolve these discrepancies.
Micro-gravity: current distributions creating a uniform force field
NASA Astrophysics Data System (ADS)
Vincent-Viry, O.; Mailfert, A.; Colteu, A.; Dael, A.; Gourdin, C.; Quettier, L.
2001-02-01
This paper presents two structures of superconducting coils able to give satisfactory solutions to the problem of generation of uniform field of high magnetic forces. The first structure is modeled by the use of purely surface current densities, whereas the second one can be described with volume current densities. Both of these structures proceed from the study of a particular expression of the complex magnetic potential introduced for structures with two-dimensional geometry. This work is carried out in a research collaboration between the GREEN and the DSM-DAPNIA department of the CEA Saclay.
On the use of quartic force fields in variational calculations
NASA Astrophysics Data System (ADS)
Fortenberry, Ryan C.; Huang, Xinchuan; Yachmenev, Andrey; Thiel, Walter; Lee, Timothy J.
2013-06-01
Quartic force fields (QFFs) have been shown to be one of the most effective ways to efficiently compute vibrational frequencies for small molecules. In this letter we discuss how the simple-internal or bond-length bond-angle (BLBA) coordinates can be transformed into Morse-cosine (-sine) coordinates which produce potential energy surfaces from QFFs that possess proper limiting behavior and can describe the vibrational (or rovibrational) energy levels of an arbitrary molecular system to 5 cm-1 or better compared to experiment. We investigate parameter scaling in the Morse coordinate, symmetry considerations, and examples of transformed QFFs making use of the MULTIMODE, TROVE, and VTET variational vibrational methods.
Validating lipid force fields against experimental data: Progress, challenges and perspectives.
Poger, David; Caron, Bertrand; Mark, Alan E
2016-07-01
Biological membranes display a great diversity in lipid composition and lateral structure that is crucial in a variety of cellular functions. Simulations of membranes have contributed significantly to the understanding of the properties, functions and behaviour of membranes and membrane-protein assemblies. This success relies on the ability of the force field used to describe lipid-lipid and lipid-environment interactions accurately, reproducibly and realistically. In this review, we present some recent progress in lipid force-field development and validation strategies. In particular, we highlight how a range of properties obtained from various experimental techniques on lipid bilayers and membranes, can be used to assess the quality of a force field. We discuss the limitations and assumptions that are inherent to both computational and experimental approaches and how these can influence the comparison between simulations and experimental data. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov. PMID:26850737
A Polarizable Multipole-based Force Field for Dimethyl and Trimethyl Phosphate
Zhang, Changsheng; Lu, Chao; Wang, Qiantao; Ponder, Jay W.; Ren, Pengyu
2016-01-01
Using quantum mechanical studies and liquid phase simulations, the AMOEBA force field for dimethylphosphate (DMP) ion and trimethylphosphate (TMP) has been developed. Based on ab initio calculations, it was found that ion binding and the solution environment significantly impact both the molecular geometry and the energy differences between conformations. Atomic multipole moments are derived from MP2/cc-pVQZ calculations of methyl phosphates at several conformations with accounting of chemical environments. Many-body polarization is handled via a Thole-style induction model using distributed atomic polarizabilities. Van der Waals parameters of phosphate and oxygen atoms are determined by fitting to the quantum mechanical interaction energy curves for water with DMP or TMP. Additional stretch-torsion and angle-torsion coupling terms were introduced in order to capture asymmetry in P-O bond lengths and angles due to the generalized anomeric effect. The resulting force field for DMP and TMP is able to accurately describe both the molecular structure and conformational energy surface, including bond and angle variations with conformation, as well as interaction of both species with water and metal ions. The force field was further validated for liquid TMP by comparing simulated density and heat of vaporization values with experimental data. Structural insight obtained from MD simulations indicates liquid TMP is stabilized by both nonpolar-nonpolar contacts and hydrogen bonding. The current study is an important step towards developing the AMOEBA model for nucleic acids. PMID:26574325
What Do Effective Field Theories Tell Us About the Nuclear Force?
NASA Astrophysics Data System (ADS)
Friar, J. L.
2003-10-01
Potentials are the tools that traditional nuclear physics uses to calculate bound states and reaction rates. First-generation potentials (developed more than a decade ago) used many different mechanisms to motivate shapes and strengths, while their fits to nucleon-nucleon scattering data were indifferent. Second-generation potentials, whose fits to data range from good to excellent, were constructed within the past decade using the same motivations. Accurate third-generation potentials are now being developed using the language and techniques of chiral perturbation theory, which is an effective field theory based on the symmetries of QCD. The philosophy of such effective field theories has made a large impact on the way I view nuclear potentials. Power counting, for example, uses the scales of strongly interacting systems to provide a systematic and concise organizational scheme that subsumes both two-nucleon and three-nucleon potentials. Although this new development has led to a rapid advance in our understanding of nuclear forces (especially in isospin violation), much of this work mirrors traditional approaches. Many of the techniques usually associated with field theory (such as regularization and renormalization, for example) also have direct analogues in more traditional approaches to constructing potentials. Examples of nuclear-force mechanisms and their associated scales will be discussed throughout the talk, together with my appreciation of how well we understand the nuclear force.
2016-01-01
Classical molecular mechanics force fields typically model interatomic electrostatic interactions with point charges or multipole expansions, which can fail for atoms in close contact due to the lack of a description of penetration effects between their electron clouds. These short-range penetration effects can be significant and are essential for accurate modeling of intermolecular interactions. In this work we report parametrization of an empirical charge–charge function previously reported (PiquemalJ.-P.; J. Phys. Chem. A2003, 107, 1035326313624) to correct for the missing penetration term in standard molecular mechanics force fields. For this purpose, we have developed a database (S101×7) of 101 unique molecular dimers, each at 7 different intermolecular distances. Electrostatic, induction/polarization, repulsion, and dispersion energies, as well as the total interaction energy for each complex in the database are calculated using the SAPT2+ method (ParkerT. M.; J. Chem. Phys.2014, 140, 09410624606352). This empirical penetration model significantly improves agreement between point multipole and quantum mechanical electrostatic energies across the set of dimers and distances, while using only a limited set of parameters for each chemical element. Given the simplicity and effectiveness of the model, we expect the electrostatic penetration correction will become a standard component of future molecular mechanics force fields. PMID:26413036
Coarse graining of force fields for metal-organic frameworks.
Dürholt, Johannes P; Galvelis, Raimondas; Schmid, Rochus
2016-03-14
We have adapted our genetic algorithm based optimization approach, originally developed to generate force field parameters from quantum mechanic reference data, to derive a first coarse grained force field for a MOF, taking the atomistic MOF-FF as a reference. On the example of the copper paddle-wheel based HKUST-1, a maximally coarse grained model, using a single bead for each three and four coordinated vertex, was developed as a proof of concept. By adding non-bonded interactions with a modified Buckingham potential, the resulting MOF-FF-CGNB is able to predict local deformation energies of the building blocks as well as bulk properties like the tbovs.pto energy difference or elastic constants in a semi-quantitative way. As expected, the negative thermal expansion of HKUST-1 is not reproduced by the maximally coarse grained model. At the expense of atomic resolution, substantially larger systems (up to tens of nanometers in size) can be simulated with respect to structural and mechanical properties, bridging the gap to the mesoscale. As an example the deformation of the [111] surface of HKUST-1 by a "tip" could be computed without artifacts from periodic images. PMID:26732756
Automated conformational energy fitting for force-field development
Guvench, Olgun; MacKerell, Alexander D.
2010-01-01
We present a general conformational-energy fitting procedure based on Monte Carlo simulated annealing (MCSA) for application in the development of molecular mechanics force fields. Starting with a target potential energy surface and an unparameterized molecular mechanics potential energy surface, an optimized set of either dihedral or grid-based correction map (CMAP) parameters is produced that minimizes the root mean squared error (RMSE) between the parameterized and targeted energies. The fitting is done using an MCSA search in parameter space and consistently converges to the same RMSE irrespective of the randomized parameters used to seed the search. Any number of dihedral parameters can be simultaneously parameterized, allowing for fitting to multi-dimensional potential energy scans. Fitting options for dihedral parameters include non-uniform weighting of the target data, constraining multiple optimized parameters to the same value, constraining parameters to be no greater than a user-specified maximum value, including all or only a subset of multiplicities defining the dihedral Fourier series, and optimization of phase angles in addition to force constants. The dihedral parameter fitting algorithm’s performance is characterized through multi-dimensional fitting of cyclohexane, tetrahydropyran, and hexopyranose monosaccharide energetics, with the latter case having a 30-dimensional parameter space. The CMAP fitting is applied in the context of polypeptides, and is used to develop a parameterization that simultaneously captures the φ, ψ energetics of the alanine dipeptide and the alanine tetrapeptide. Because the dihedral energy term is common to many force fields, we have implemented the dihedral-fitting algorithm in the portable Python scripting language and have made it freely available as Supplementary Material. PMID:18458967
Grid-based backbone correction to the ff12SB protein force field for implicit-solvent simulations.
Perez, Alberto; MacCallum, Justin L; Brini, Emiliano; Simmerling, Carlos; Dill, Ken A
2015-10-13
Force fields, such as Amber's ff12SB, can be fairly accurate models of the physical forces in proteins and other biomolecules. When coupled with accurate solvation models, force fields are able to bring insight into the conformational preferences, transitions, pathways, and free energies for these biomolecules. When computational speed/cost matters, implicit solvent is often used but at the cost of accuracy. We present an empirical grid-like correction term, in the spirit of cMAPs, to the combination of the ff12SB protein force field and the GBneck2 implicit-solvent model. Ff12SB-cMAP is parametrized on experimental helicity data. We provide validation on a set of peptides and proteins. Ff12SB-cMAP successfully improves the secondary structure biases observed in ff12SB + Gbneck2. Ff12SB-cMAP can be downloaded ( https://github.com/laufercenter/Amap.git ) and used within the Amber package. It can improve the agreement of force fields + implicit solvent with experiments. PMID:26574266
Beyond Born-Mayer: Improved Models for Short-Range Repulsion in ab Initio Force Fields.
Van Vleet, Mary J; Misquitta, Alston J; Stone, Anthony J; Schmidt, J R
2016-08-01
Short-range repulsion within intermolecular force fields is conventionally described by either Lennard-Jones (A/r(12)) or Born-Mayer (A exp(-Br)) forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of intermolecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, and robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Finally, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach. PMID:27337546
Beyond Born-Mayer: Improved models for short-range repulsion in ab initio force fields
Van Vleet, Mary J.; Misquitta, Alston J.; Stone, Anthony J.; Schmidt, Jordan R.
2016-06-23
Short-range repulsion within inter-molecular force fields is conventionally described by either Lennard-Jones or Born-Mayer forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of inter-molecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, andmore » robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Lastly, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.« less
Wildman, Jack; Repiščák, Peter; Paterson, Martin J; Galbraith, Ian
2016-08-01
We describe a general scheme to obtain force-field parameters for classical molecular dynamics simulations of conjugated polymers. We identify a computationally inexpensive methodology for calculation of accurate intermonomer dihedral potentials and partial charges. Our findings indicate that the use of a two-step methodology of geometry optimization and single-point energy calculations using DFT methods produces potentials which compare favorably to high level theory calculation. We also report the effects of varying the conjugated backbone length and alkyl side-chain lengths on the dihedral profiles and partial charge distributions and determine the existence of converged lengths above which convergence is achieved in the force-field parameter sets. We thus determine which calculations are required for accurate parametrization and the scope of a given parameter set for variations to a given molecule. We perform simulations of long oligomers of dioctylfluorene and hexylthiophene in explicit solvent and find peristence lengths and end-length distributions consistent with experimental values. PMID:27397762