Elasticity-driven interaction between vortices in type-II superconductors

The contribution to the vortex lattice energy which is due to the vortex-induced strains is calculated covering all the magnetic field range which defines the vortex state. This contribution is compared with previously reported ones what shows that, in the most part of the vortex state, it has been notably underestimated until now. The reason of such underestimation is the assumption that only the vortex cores induce strains. In contrast to what is generally assumed, both core and non-core regions are important sources of strains in high-$\kappa$ superconductors.Comment: 10 pages, 1 figure, revtex

Maximum Flux Transition Paths of Conformational Change

Given two metastable states A and B of a biomolecular system, the problem is to calculate the likely paths of the transition from A to B. Such a calculation is more informative and more manageable if done for a reduced set of collective variables chosen so that paths cluster in collective variable space. The computational task becomes that of computing the "center" of such a cluster. A good way to define the center employs the concept of a committor, whose value at a point in collective variable space is the probability that a trajectory at that point will reach B before A. The committor "foliates" the transition region into a set of isocommittors. The maximum flux transition path is defined as a path that crosses each isocommittor at a point which (locally) has the highest crossing rate of distinct reactive trajectories. (This path is different from that of the MaxFlux method of Huo and Straub.) It is argued that such a path is nearer to an ideal path than others that have been proposed with the possible exception of the finite-temperature string method path. To make the calculation tractable, three approximations are introduced, yielding a path that is the solution of a nonsingular two-point boundary-value problem. For such a problem, one can construct a simple and robust algorithm. One such algorithm and its performance is discussed.Comment: 7 figure

The regulatory subunit of PKA-I remains partially structured and undergoes β-aggregation upon thermal denaturation

Background: The regulatory subunit (R) of cAMP-dependent protein kinase (PKA) is a modular flexible protein that responds with large conformational changes to the binding of the effector cAMP. Considering its highly dynamic nature, the protein is rather stable. We studied the thermal denaturation of full-length RIα and a truncated RIα(92-381) that contains the tandem cyclic nucleotide binding (CNB) domains A and B. Methodology/Principal Findings: As revealed by circular dichroism (CD) and differential scanning calorimetry, both RIα proteins contain significant residual structure in the heat-denatured state. As evidenced by CD, the predominantly α-helical spectrum at 25°C with double negative peaks at 209 and 222 nm changes to a spectrum with a single negative peak at 212-216 nm, characteristic of β-structure. A similar α→β transition occurs at higher temperature in the presence of cAMP. Thioflavin T fluorescence and atomic force microscopy studies support the notion that the structural transition is associated with cross-β-intermolecular aggregation and formation of non-fibrillar oligomers. Conclusions/Significance: Thermal denaturation of RIα leads to partial loss of native packing with exposure of aggregation-prone motifs, such as the B' helices in the phosphate-binding cassettes of both CNB domains. The topology of the β-sandwiches in these domains favors inter-molecular β-aggregation, which is suppressed in the ligand-bound states of RIα under physiological conditions. Moreover, our results reveal that the CNB domains persist as structural cores through heat-denaturation. © 2011 Dao et al

Ground State Vortex Lattice Structures in d-wave Superconductors

We show in a realistic $d_{x^{2}-y^{2}}$ symmetry gap model for a cuprate superconductor that the clean vortex lattice has discontinuous structural transitions (at and near T=0), as a function of the magnetic field $B$ along the c-axis. The transitions arise from the singular nonlocal and anisotropic susceptibility of the $d_{x^{2}-y^{2}}$ superconductor to the perturbation caused by supercurrents associated with vortices. The susceptibility, due to virtual Dirac quasiparticle-hole excitation, is calculated carefully, and leads to a ground state transition for the triangular lattice from an orientation along one of the crystal axis to one at 45$^o$ to them, i.e, along the gap zero direction. The field scale is seen to be 5 Tesla $\sim (\Delta_{0}/ta)^{2}\Phi_{0}$, where $\Delta_{0}$ is the gap maximum, $t$ is the nearest neighbour hopping, $a$ is the lattice constant, and $\Phi_{0}$ is the flux quantum. At much higher fields ($\sim 28T$) there is a discontinuous transition to a centred square structure. The source of the differences from existing calculations, and experimental observability are discussed, the latter especially in view of the very small (a few degrees $K$ per vortex) differences in the ground state energy.Comment: To be published in Phys. Rev.

Exploring the conformational dynamics of alanine dipeptide in solution subjected to an external electric field: A nonequilibrium molecular dynamics simulation

In this paper, we investigate the conformational dynamics of alanine dipeptide under an external electric field by nonequilibrium molecular dynamics simulation. We consider the case of a constant and of an oscillatory field. In this context we propose a procedure to implement the temperature control, which removes the irrelevant thermal effects of the field. For the constant field different time-scales are identified in the conformational, dipole moment, and orientational dynamics. Moreover, we prove that the solvent structure only marginally changes when the external field is switched on. In the case of oscillatory field, the conformational changes are shown to be as strong as in the previous case, and non-trivial nonequilibrium circular paths in the conformation space are revealed by calculating the integrated net probability fluxes.Comment: 23 pages, 12 figure

Volume-energy correlations in the slow degrees of freedom of computer-simulated phospholipid membranes

Constant-pressure molecular-dynamics simulations of phospholipid membranes in the fluid phase reveal strong correlations between equilibrium fluctuations of volume and energy on the nanosecond time-scale. The existence of strong volume-energy correlations was previously deduced indirectly by Heimburg from experiments focusing on the phase transition between the fluid and the ordered gel phases. The correlations, which are reported here for three different membranes (DMPC, DMPS-Na, and DMPSH), have volume-energy correlation coefficients ranging from 0.81 to 0.89. The DMPC membrane was studied at two temperatures showing that the correlation coefficient increases as the phase transition is approached

Experimental investigation of the stochastic early flame propagation after ignition by a low-energy electrical discharge

In the context of explosion protection, very conservative safety factors need to be considered, e.g. in the design of electrical devices. This is due to standards which are mainly based on empirical data as opposed to a detailed knowledge of the underlying physiochemical processes. In this work, the early phase of ignition of burnable gas mixtures close to their respective minimum ignition energy is investigated experimentally by means of high-speed schlieren imaging. Our data quantifies how the ignition process at such low energies becomes less repeatable which is evidenced by a high scattering of the flame propagation. It was found that, depending on the mixture, the flow field induced by the electrical discharge may exhibit a considerable effect on the ignition process. This effect is more pronounced for mixtures which are characterized by a large Lewis number, thus, leading to a more random flame propagation

Catalytic Cycle of Multicopper Oxidases Studied by Combined Quantum- and Molecular-Mechanical Free-Energy Perturbation Methods

We have used combined quantum mechanical and molecular mechanical free-energy perturbation methods in combination with explicit solvent simulations to study the reaction mechanism of the multicopper oxidases, in particular the regeneration of the reduced state from the native intermediate. For 52 putative states of the trinuclear copper cluster, differing in the oxidation states of the copper ions and the protonation states of water- and O2-derived ligands, we have studied redox potentials, acidity constants, isomerisation reactions, as well as water- and O2 binding reactions. Thereby, we can propose a full reaction mechanism of the multicopper oxidases with atomic detail. We also show that the two copper sites in the protein communicate so that redox potentials and acidity constants of one site are affected by up to 0.2 V or 3 pKa units by a change in the oxidation state of the other site