Optimization problems of the stress-deformed state of thermoelastic shells

We discuss general questions about the choice of optimization criteria for the stressed state of thermoelastic shells taking account of the functional measures of the space-time variability of its characteristics. The efficiency of the proposed approach is illustrated using the example of the solution of a problem for an axisymmetrically heated cylindrical shell.Translated fromMatematicheskie Metody i Fiziko-Mekhanicheskie Polya, Issue 27, 1988, pp. 11–18.     

Toward separation of nuclear spin isomers with coherent light.

We propose an approach for separating nuclear spin isomers with coherent light and illustrate it by numerical calculations using fulvene as a model system. The scheme employs the equivalence of torsion and interchange of equivalent H-atoms in a class of molecules of which fulvene is a simple example. The exchange symmetry couples with the rotational symmetry to produce a spatial distinction between the two photo-excited nuclear spin isomers, and wavepacket interferometry is applied to separate the species.     

A recipe for evaluating and differentiating cos ϕ expressions

We present a simple recipe for calculating and differentiating cosine of bond angle and dihedral angle expressions. The resulting formulas can be incorporated in a straightforward manner into the bond angle and dihedral angle components of potential energy functions. These formulas rely only on expressions and derivatives of dot products, and, in particular, they avoid cross products as well as excessive Fortran function references. Consequently, the expressions derived in this article can be written compactly and evaluated rapidly.     

Hierarchical structure of the energy landscape of proteins revisited by time series analysis. I. Mimicking protein dynamics in different time scales

Time series models, which are constructed from the projections of the molecular-dynamics (MD) runs on principal components (modes), are used to mimic the dynamics of two proteins: tendamistat and immunity protein of colicin E7 (ImmE7). Four independent MD runs of tendamistat and three independent runs of ImmE7 protein in vacuum are used to investigate the energy landscapes of these proteins. It is found that mean-square displacements of residues along the modes in different time scales can be mimicked by time series models, which are utilized in dividing protein dynamics into different regimes with respect to the dominating motion type. The first two regimes constitute the dominance of intraminimum motions during the first 5 ps and the random walk motion in a hierarchically higher-level energy minimum, which comprise the initial time period of the trajectories up to 20-40 ps for tendamistat and 80-120 ps for ImmE7. These are also the time ranges within which the linear nonstationary time series are completely satisfactory in explaining protein dynamics. Encountering energy barriers enclosing higher-level energy minima constrains the random walk motion of the proteins, and pseudorelaxation processes at different levels of minima are detected in tendamistat, depending on the sampling window size. Correlation (relaxation) times of 30-40 ps and 150-200 ps are detected for two energy envelopes of successive levels for tendamistat, which gives an overall idea about the hierarchical structure of the energy landscape. However, it should be stressed that correlation times of the modes are highly variable with respect to conformational subspaces and sampling window sizes, indicating the absence of an actual relaxation. The random-walk step sizes and the time length of the second regime are used to illuminate an important difference between the dynamics of the two proteins, which cannot be clarified by the investigation of relaxation times alone: ImmE7 has lower-energy barriers enclosing the higher-level energy minimum, preventing the protein to relax and letting it move in a random-walk fashion for a longer period of time.     

The Structure of Silylated Amides: N-Methyl-N-Trimethylsilyltrifluoroacetamide,a Reassignment of structure

The ²⁹Si and ¹⁴N nmr spectra of some bissilylamides and monosilylamides are reported. The results are compared with those for model silylesters, amines and imines. The compound of formula CF₃CONCH₃Si(CH₃)₃ (MSTFA) is shown to exist as a tautemeric mixture of a silylamide and a silylimidate. This assignment is confirmed by ¹³C nmr. A rationale is proposed for the structure directing effect of substituents on the amide carbon atom.     

Using spin dynamics of covalently linked radical ion pairs to probe the impact of structural and energetic changes on charge recombination

We have synthesized a series of structurally related, covalently linked electron donor-acceptor triads having highly restricted conformations to study the effects of radical ion pair (RP) structure, energetics, and solvation on charge recombination. The chromophoric electron acceptor in these triads is a 4-aminonaphthalene-1,8-dicarboximide (6ANI), in which the 4-amine nitrogen atom is part of a piperazine ring. The second nitrogen atom of the piperazine ring is part of a para-substituted aniline donor, where the para substituents are X = H, OMe, and NMe(2). The imide group of 6ANI is linked to a naphthalene-1,8:4,5-bis(dicarboximide) (NI) electron acceptor across a phenyl spacer in a meta relationship. The triads undergo two-step photoinduced electron transfer to yield their respective XAn(*)(+)-6ANI-Ph-NI(*)(-) RP states, which undergo radical pair intersystem crossing followed by charge recombination to yield (3)NI. Time-resolved electron paramagnetic resonance experiments on the spin-polarized RPs and triplet states carried out in toluene and in E-7, a mixture of nematic liquid crystals (LCs), show that for all three triads, the XAn(*)(+)-6ANI-Ph-NI(*)(-) RPs are correlated radical pairs and directly yield values of the spin-spin exchange interaction, J, and the dipolar interaction, D. The values of J are all about -1 mT and show that the LC environment most likely enforces the chair conformation at the piperazine ring, for which the RP distance is larger than that for the corresponding boat conformation. The values of D yield effective RP distances that agree well with those calculated earlier from the spin distributions of the radical ions. Within the LC, changing the temperature shows that the CR mechanism can be changed significantly as the energy levels of the RPs change relative to that of the recombination triplet.     

A powerful truncated Newton method for potential energy minimization

With advances in computer architecture and software, Newton methods are becoming not only feasible for large-scale nonlinear optimization problems, but also reliable, fast and efficient. Truncated Newton methods, in particular, are emerging as a versatile subclass. In this article we present a truncated Newton algorithm specifically developed for potential energy minimization. The method is globally convergent with local quadratic convergence. Its key ingredients are: (1) approximation of the Newton direction far away from local minima, (2) solution of the Newton equation iteratively by the linear Conjugate Gradient method, and (3) preconditioning of the Newton equation by the analytic second-derivative components of the “local” chemical interactions: bond length, bond angle and torsional potentials. Relaxation of the required accuracy of the Newton search direction diverts the minimization search away from regions where the function is nonconvex and towards physically interesting regions. The preconditioning strategy significantly accelerates the iterative solution for the Newton search direction, and therefore reduces the computation time for each iteration. With algorithmic variations, the truncated Newton method can be formulated so that storage and computational requirements are comparable to those of the nonlinear Conjugate Gradient method. As the convergence rate of nonlinear Conjugate Gradient methods is linear and performance less predictable, the application of the truncated Newton code to potential energy functions is promising.     

Inherent speedup limitations in multiple time step/particle mesh Ewald algorithms

Multiple time step (MTS) algorithms present an effective integration approach to reduce the computational cost of dynamics simulations. By using force splitting to allow larger time steps for the more slowly varying force components, computational savings can be realized. The Particle-Mesh-Ewald (PME) method has been independently devised to provide an effective and efficient treatment of the long-range electrostatics interactions. Here we examine the performance of a combined MTS/PME algorithm previously developed for AMBER on a large polymerase beta/DNA complex containing 40,673 atoms. Our goal is to carefully combine the robust features of the Langevin/MTS (LN) methodology implemented in CHARMM-which uses position rather than velocity Verlet with stochasticity to make possible outer time steps of 150 fs-with the PME formulation. The developed MTS/PME integrator removes fast terms from the reciprocal-space Ewald component by using switch functions. We analyze the advantages and limitations of the resulting scheme by comparing performance to the single time step leapfrog Verlet integrator currently used in AMBER by evaluating different time-step protocols using three assessors for accuracy, speedup, and stability, all applied to long (i.e., nanosecond) simulations to ensure proper energy conservation. We also examine the performance of the algorithm on a parallel, distributed shared-memory computer (SGI Origin 2000 with 8 300-MHz R12000 processors). Good energy conservation and stability behavior can be demonstrated, for Newtonian protocols with outer time steps of up to 8 fs and Langevin protocols with outer time steps of up to 16 fs. Still, we emphasize the inherent limitations imposed by the incorporation of MTS methods into the PME formulation that may not be widely appreciated. Namely, the limiting factor on the largest outer time-step size, and hence speedup, is an intramolecular cancellation error inherent to PME. This error stems from the excluded-nonbonded correction term contained in the reciprocal-space component. This cancellation error varies in time and introduces artificial frequencies to the governing dynamics motion. Unfortunately, we find that this numerical PME error cannot be easily eliminated by refining the PME parameters (grid resolution and/or order of interpolating polynomial). We suggest that methods other than PME for fast electrostatics may allow users to reap the full advantages from MTS algorithms.     

Substitution-Inert Metal Complex Mobile Phases in Non-Suppressed Cation Chromatography

Substitution-inert metal complexes, [Co(edda)(H₂O)₂]⁺, (Co(edda)(en)]⁺, [Co(edda)(dmen)]⁺, [Co(en)₂-(gly)]²⁺, [Co(en)₂(acac)]²⁺, and [Co(trien)(gly)]²⁺ in their nitrate salt solutions are used as eluents in nonsuppressed cation chromatography (where edda = ethylenediamine-N,N′-diacetic acid, en = ethylenediamine, dmen = N,N′-dime-thylethylenediamine, gly = glycine, acac = acetylacetone, and trien = triethylenetetraamine). It is found that all the mono- and di-valent charged complexes can be used to separate alkali and alkaline earth metal cations, respectively. The separations for monovalent cations are sometimes comparable to those using ultrapure HNO₃ solutions. However, the divalent Ca²⁺ and Sr²⁺ ions cannot be resolved using the metal complex eluents. On the other hand, a selected, direct non-suppressed IC separation of zinc(II) and cadmium(II) ions is demonstrated for the first time using a substitution-inert metal complex eluent and a conductivity detector. Comparisons of these eluents with those reported previously, i. e. HNO₃ and ethylenediammonium salt solution are made and explanations are proposed to account for the different selectivities observed where possible. The future development of this technique is also briefly discussed.