Sharpening the focus of force field analysis

The purpose is to explore the inherent complexity of Kurt Lewin's force field theory through applied analysis of organizational case examples and related methods. The methodology applies a range of tools from the consultancy research domain, including force field analysis of complex organizational scenarios, and applies bricolage and corroboration to emerging discoveries from semi-structured interviews, author experience, critical reflection and literature survey. Findings are that linear representation of internal and external forces in organizational applications of field theory does not fully explain the paradox of inverse vectors in the forces of change. The force field is not an impermeable thing; instead, it morphs. Examples of the inverse principle and its effects are detailed and extended in this analysis. The implications of the research are that force field analysis and related change processes promoted in organizational change literature run the risk of missing key complexities. The inclusion of the inverse principle can provide enhanced, holistic understanding of the prevailing forces for change. The augmentation of the early work of Kurt Lewin, and extension of previous analyses of his legacy in the Journal of Change Management and elsewhere, provide, in this article, change analysis insights that align well with current organizational environments.

First-principles based force-field for the interaction of proteins with Au(100)(5 x 1): an extension of GolP-CHARMM

Noncovalent recognition between peptides and inorganic materials is an established phenomenon. Key to exploiting these interactions in a wide range of materials self-assembly applications would be to harness the facet-selective control of peptide binding onto these materials. Fundamental understanding of what drives facet-selectivity in peptide binding is developing, but as yet is not sufficient to enable design of predictable facet-specific sequences. Computational simulation of the aqueous peptide-gold interface, commonly used to understand the mechanisms driving adsorption at an atomic level, has thus far neglected the role that surface reconstruction might play in facet specificity. Here the polarizable GolP-CHARMM suite of force fields is extended to include the reconstructed Au(100) surface. The force field, compatible with the bio-organic force field CHARMM, is parametrized using first-principles data. Our extended force field is tailored to reproduce the heterogeneity of weak chemisorbing N and S species to specific locations in the Au(100)(5 × 1) surface identified from the first-principles calculations. We apply our new model to predict and compare the three-dimensional structure of liquid water at Au(111), Au(100)(1 × 1), and Au(100)(5 × 1) interfaces. Using molecular dynamics simulations, we predict an increased likelihood for water-mediated peptide adsorption at the aqueous-Au(100)(1 × 1) interface compared with the Au(100)(5 × 1) interface. Therefore, our findings suggest that peptide binding can discriminate between the native and reconstructed Au(100) interfaces and that the role of reconstruction on binding at the Au(100) interface should not be neglected. © 2013 American Chemical Society.

Molecular dynamics study of response of liquid N,N-dimethylformamide to externally applied electric field using a polarizable force field

The behavior of Liquid N,N-dimethylformamide subjected to a wide range of externally applied electric fields (from 0.001 V/nm to 1 V/nm) has been investigated through molecular dynamics simulation. To approach the objective the AMOEBA polarizable force field was extended to include the interaction of the external electric field with atomic partial charges and the contribution to the atomic polarization. The simulation results were evaluated with quantum mechanical calculations. The results from the present force field for the liquid at normal conditions were compared with the experimental and molecular dynamics results with non-polarizable and other polarizable force fields. The uniform external electric fields of higher than 0.01 V/nm have a significant effect on the structure of the liquid, which exhibits a variation in numerous properties, including molecular polarization, local cluster structure, rotation, alignment, energetics, and bulk thermodynamic and structural properties.

Optimal space manoeuvres in a non-Keplerian force field

In this work the problem of a spacecraft bi-impulsive transfer between two given non coplanar elliptical orbits, with minimum fuel consumption, is solved considering a non-Keplerian force field (the perturbing forces include Earth gravity harmonics and atmospheric drag). The problem is transformed in the Two Point Boundary Value Problem. It is developed and implemented a new algorithm, that uses the analytical expressions developed here. A dynamics that considered a Keplerian force field was used to produce an initial guess to solve the Two Point Boundary Value Problem. Several simulations were performed to observe the spacecraft orbital behaviour by different kind of perturbations and constraints, on a fuel consumption optimization point of view. (C) 2002 COSPAR. Published by Elsevier B.V. Ltd. All rights reserved.

Ramachandran-type Plots for Glycosidic Linkages: Examples from Molecular Dynamics Simulations using the Glycam06 Force Field

The goals of this article are to (1) provide further validation of the Glycam06 force field, specifically for its use in implicit solvent molecular dynamic (MD) simulations, and (2) to present the extension of G.N. Ramachandran's idea of plotting amino acid phi and psi angles to the glycosidic phi, psi, and omega angles formed between carbohydrates. As in traditional Ramachandran plots, these carbohydrate Ramachandran-type (carb-Rama) plots reveal the coupling between the glycosidic angles by displaying the allowed and disallowed conformational space. Considering two-bond glycosidic linkages, there are 18 possible conformational regions that can be defined by (α, ϕ, ψ) and (β, ϕ, ψ), whereas for three-bond linkages, there are 54 possible regions that can be defined by (α, ϕ, ψ, ω) and (β, ϕ, ψ, ω). Illustrating these ideas are molecular dynamic simulations on an implicitly hydrated oligosaccharide (700 ns) and its eight constituent disaccharides (50 ns/disaccharide). For each linkage, we compare and contrast the oligosaccharide and respective disaccharide carb-Rama plots, validate the simulations and the Glycam06 force field through comparison to experimental data, and discuss the general trends observed in the plots.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Cost-Effective Force Field Tailored for Solid-Phase Simulations of OLED Materials

A united atom force field is empirically derived by minimizing the difference between experimental and simulated crystal cells and melting temperatures for eight compounds representative of organic electronic materials used in OLEDs and other devices: biphenyl, carbazole, fluorene, 9,9′-(1,3-phenylene)bis(9H-carbazole)-1,3-bis(N-carbazolyl)benzene (mCP), 4,4′-bis(N-carbazolyl)-1,1′-biphenyl (pCBP), phenazine, phenylcarbazole, and triphenylamine. The force field is verified against dispersion-corrected DFT calculations and shown to also successfully reproduce the crystal structure for two larger compounds employed as hosts in phosphorescent and thermally activated delayed fluorescence OLEDs: N,N′-di(1-naphthyl)-N,N′-diphenyl-(1,1′-biphenyl)-4,4′-diamine (NPD), and 1,3,5-tri(1-phenyl-1H-benzo[d]imidazol-2-yl)phenyl (TPBI). The good performances of the force field coupled to the large computational savings granted by the united atom approximation make it an ideal choice for the simulation of the morphology of emissive layers for OLED materials in crystalline or glassy phases.

Optical force field mapping in microdevices

We present a method for characterizing microscopic optical force fields. Two dimensional vector force maps are generated by measuring the optical force applied to a probe particle for a grid of particle positions. The method is used to map Out the force field created by the beam from a lensed fiber inside a liquid filled microdevice. We find transverse gradient forces and axial scattering forces on the order of 2 pN per 10 mW laser power which are constant over a considerable axial range (> 35 mu m). These findings suggest Future useful applications of lensed fibers for particle guiding/sorting. The propulsion of a small particle at a constant velocity of 200 mu m s(-1) is shown.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Field emission properties of carbon nanotube arrays with defects and impurities

Field emission from carbon nanotubes (CNTs) in the form of arrays or thin films give rise to several strongly correlated process of electromechanical interaction and degradation. Such processes are mainly due to (1) electron-phonon interaction (2) electromechanical force field leading to stretching of CNTs (3) ballistic transport induced thermal spikes, coupled with high dynamic stress, leading to degradation of emission performance at the device scale. Fairly detailed physics based models of CNTs considering the aspects (1) and (2) above have already been developed by these authors, and numerical results indicate good agreement with experimental results. What is missing in such a system level modeling approach is the incorporation of structural defects and vacancies or charge impurities. This is a practical and important problem due to the fact that degradation of field emission performance is indeed observed in experimental I-V curves. What is not clear from these experiments is whether such degradation in the I-V response is due to dynamic reorientation of the CNTs or due to the defects or due to both of these effects combined. Non-equilibrium Green’s function based simulations using a tight-binding Hamiltonian for single CNT segment show up the localization of carrier density at various locations of the CNTs. About 11% decrease in the drive current with steady difference in the drain current in the range of 0.2-0.4V of the gate voltage was reported in literature when negative charge impurity was introduced at various locations of the CNT over a length of ~20nm. In the context of field emission from CNT tips, a simplistic estimate of defects have been introduced by a correction factor in the Fowler-Nordheim formulae. However, a more detailed physics based treatment is required, while at the same time the device-scale simulation is necessary. The novelty of our present approach is the following. We employ a concept of effective stiffness degradation for segments of CNTs, which is due to structural defects, and subsequently, we incorporate the vacancy defects and charge impurity effects in the Green’s function based approach. Field emission induced current-voltage characteristics of a vertically aligned CNT array on a Cu-Cr substrate is then simulated using a detailed nonlinear mechanistic model of CNTs coupled with quantum hydrodynamics. An array of 10 vertically aligned and each 12 m long CNTs is considered for the device scale analysis. Defect regions are introduced randomly over the CNT length. The result shows the decrease in the longitudinal strain due to defects. Contrary to the expected influence of purely mechanical degradation, this result indicates that the charge impurity and hence weaker transport can lead to a different electromechanical force field, which ultimately can reduce the strain. However, there could be significant fluctuation in such strain field due to electron-phonon coupling. The effect of such fluctuations (with defects) is clearly evident in the field emission current history. The average current also decreases significantly due to such defects.

Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

FIELD EMISSION EFFICIENCY OF A CARBON NANOTUBE ARRAY UNDER PARASITIC NONLINEARITIES

Carbon Nanotubes (CNTs) grown on substrates are potential electron sources in field emission applications. Several studies have reported the use of CNTs in field emission devices, including field emission displays, X-ray tube, electron microscopes, cathode-ray lamps, etc. Also, in recent years, conventional cold field emission cathodes have been realized in micro-fabricated arrays for medical X-ray imaging. CNTbased field emission cathode devices have potential applications in a variety of industrial and medical applications, including cancer treatment. Field emission performance of a single isolated CNT is found to be remarkable, but the situation becomes complex when an array of CNTs is used. At the same time, use of arrays of CNTs is practical and economical. Indeed, such arrays on cathode substrates can be grown easily and their collective dynamics can be utilized in a statistical sense such that the average emission intensity is high enough and the collective dynamics lead to longer emission life. The authors in their previous publications had proposed a novel approach to obtain stabilized field emission current from a stacked CNT array of pointed height distribution. A mesoscopic modeling technique was employed, which took into account electro-mechanical forces in the CNTs, as well as transport of conduction electron coupled with electron phonon induced heat generation from the CNT tips. The reported analysis of pointed arrangements of the array showed that the current density distribution was greatly localized in the middle of the array, the scatter due to electrodynamic force field was minimized, and the temperature transients were much smaller compared to those in an array with random height distribution. In the present paper we develop a method to compute the emission efficiency of the CNT array in terms of the amount of electrons hitting the anode surface using trajectory calculations. Effects of secondary electron emission and parasitic capacitive nonlinearity on the current-voltage signals are accounted. Field emission efficiency of a stacked CNT array with various pointed height distributions are compared to that of arrays with random and uniform height distributions. Effect of this parasitic nonlinearity on the emission switch-on voltage is estimated by model based simulation and Monte Carlo method.

A body force-based method for prediction of multiple-pure-tone noise: Validation

Emissions, fuel burn, and noise are the main drivers for innovative aircraft design. Embedded propulsion systems, such as for example used in hybrid-wing body aircraft, can offer fuel burn and noise reduction benefits but the impact of inlet flow distortion on the generation and propagation of turbomachinery noise has yet to be assessed. A novel approach is used to quantify the effects of non-uniform flow on the creation and propagation of multiple pure tone (MPT) noise. The ultimate goal is to conduct a parametric study of S-duct inlets to quantify the effects of inlet design parameters on the acoustic signature. The key challenge is that the effects of distortion transfer, noise source generation and propagation through the non-uniform flow field are inherently coupled such that a simultaneous computation of the aerodynamics and acoustics is required to capture the mechanisms at play. The technical approach is based on a body force description of the fan blade row that is able to capture the distortion transfer and the blade-to-blade flow variations that cause the MPT noise while reducing computational cost. A single, 3-D full-wheel CFD simulation, in which the Euler equations are solved to second-order spatial and temporal accuracy, simultaneously computes the MPT noise generation and its propagation in distorted inlet flow. A new method of producing the blade-to-blade variations in the body force field for MPT noise generation has been developed and validated. The numerical dissipation inherent to the solver is quantified and used to correct for non-physical attenuation in the far-field noise spectra. Source generation, acoustic propagation and acoustic energy transfer between modes is examined in detail. The new method is validated on NASA's Source Diagnostic Test fan and inlet, showing good agreement with experimental data for aerodynamic performance, acoustic source generation, and far-field noise spectra. The next steps involve the assessment of MPT noise in serpentine inlet ducts and the development of a reduced order formulation suitable for incorporation into NASA's ANOPP framework. © 2010 by Jeff Defoe, Alex Narkaj & Zoltan Spakovszky.

Field configured assembly: Programmed manipulation and self-assembly at the mesoscale

Field configured assembly is a programmable force field method that permits rapid, "hands-free" manipulation, assembly, and integration of mesoscale objects and devices. In this method, electric fields, configured by specific addressing of receptor and counter electrode sites pre-patterned at a silicon chip substrate, drive the field assisted transport, positioning, and localization of mesoscale devices at selected receptor locations. Using this approach, we demonstrate field configured deterministic and stochastic self-assembly of model mesoscale devices, i.e., 50 mum diameter, 670 nm emitting GaAs-based light emitting diodes, at targeted receptor sites on a silicon chip. The versatility of the field configured assembly method suggests that it is applicable to self-assembly of a wide variety of functionally integrated nanoscale and mesoscale systems.

Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.