A non-perturbative renormalization group study of the stochastic NavierStokes equation

We study the renormalization group flow of the average action of the stochastic Navier-Stokes equation with power-law forcing. Using Galilean invariance, we introduce a nonperturbative approximation adapted to the zero-frequency sector of the theory in the parametric range of the Hölder exponent 4−2 ɛ of the forcing where real-space local interactions are relevant. In any spatial dimension d, we observe the convergence of the resulting renormalization group flow to a unique fixed point which yields a kinetic energy spectrum scaling in agreement with canonical dimension analysis. Kolmogorov's −5/3 law is, thus, recovered for ɛ=2 as also predicted by perturbative renormalization. At variance with the perturbative prediction, the −5/3 law emerges in the presence of a saturation in the ɛ dependence of the scaling dimension of the eddy diffusivity at ɛ=3/2 when, according to perturbative renormalization, the velocity field becomes infrared relevant.

Stability of vacancy clusters in FeCr alloys: A study of the Cr concentration dependence

Fe–Cr based alloys are the leading structural material candidates in the design of next generation reactors due to their high resistance to swelling and corrosion. Despite these good properties there are others, such as embrittlement, which require a higher level of understanding in order to improve aspects such as safety or lifetime of the reactors. The addition of Cr improves the behavior of the steels under irradiation, but not in a monotonic way. Therefore, understanding the changes in the Fe–Cr based alloys microstructure induced by irradiation and the role played by the alloying element (Cr) is needed in order to predict the response of these materials under the extreme conditions they are going to support. In this work we perform a study of the effect of Cr concentration in a bcc Fe–Cr matrix on formation and binding energies of vacancy clusters up to 5 units. The dependence of the calculated formation and binding energy is investigated with two empirical interatomic potentials specially developed to study radiation damage in Fe–Cr alloys. Results are very similar for both potentials showing an increase of the defect stability with the cluster size and no real dependence on Cr concentration for the binding energy.

Multiscale modeling of the plastic behaviour in single crystal tungsten: from atomistic to crystal plasticity simulations

En una planta de fusión, los materiales en contacto con el plasma así como los materiales de primera pared experimentan condiciones particularmente hostiles al estar expuestos a altos flujos de partículas, neutrones y grandes cargas térmicas. Como consecuencia de estas diferentes y complejas condiciones de trabajo, el estudio, desarrollo y diseño de estos materiales es uno de los más importantes retos que ha surgido en los últimos años para la comunidad científica en el campo de los materiales y la energía. Debido a su baja tasa de erosión, alta resistencia al sputtering, alta conductividad térmica, muy alto punto de fusión y baja retención de tritio, el tungsteno (wolframio) es un importante candidato como material de primera pared y como posible material estructural avanzado en fusión por confinamiento magnético e inercial. Sin embargo, el tiempo de vida del tungsteno viene controlado por diversos factores como son su respuesta termo-mecánica en la superficie, la posibilidad de fusión y el fallo por acumulación de helio. Es por ello que el tiempo de vida limitado por la respuesta mecánica del tungsteno (W), y en particular su fragilidad, sean dos importantes aspectos que tienes que ser investigados. El comportamiento plástico en materiales refractarios con estructura cristalina cúbica centrada en las caras (bcc) como el tungsteno está gobernado por las dislocaciones de tipo tornillo a escala atómica y por conjuntos e interacciones de dislocaciones a escalas más grandes. El modelado de este complejo comportamiento requiere la aplicación de métodos capaces de resolver de forma rigurosa cada una de las escalas. El trabajo que se presenta en esta tesis propone un modelado multiescala que es capaz de dar respuestas ingenieriles a las solicitudes técnicas del tungsteno, y que a su vez está apoyado por la rigurosa física subyacente a extensas simulaciones atomísticas. En primer lugar, las propiedades estáticas y dinámicas de las dislocaciones de tipo tornillo en cinco potenciales interatómicos de tungsteno son comparadas, determinando cuáles de ellos garantizan una mayor fidelidad física y eficiencia computacional. Las grandes tasas de deformación asociadas a las técnicas de dinámica molecular hacen que las funciones de movilidad de las dislocaciones obtenidas no puedan ser utilizadas en los siguientes pasos del modelado multiescala. En este trabajo, proponemos dos métodos alternativos para obtener las funciones de movilidad de las dislocaciones: un modelo Monte Cario cinético y expresiones analíticas. El conjunto de parámetros necesarios para formular el modelo de Monte Cario cinético y la ley de movilidad analítica son calculados atomísticamente. Estos parámetros incluyen, pero no se limitan a: la determinación de las entalpias y energías de formación de las parejas de escalones que forman las dislocaciones, la parametrización de los efectos de no Schmid característicos en materiales bcc,etc. Conociendo la ley de movilidad de las dislocaciones en función del esfuerzo aplicado y la temperatura, se introduce esta relación como ecuación de flujo dentro de un modelo de plasticidad cristalina. La predicción del modelo sobre la dependencia del límite de fluencia con la temperatura es validada experimentalmente con ensayos uniaxiales en tungsteno monocristalino. A continuación, se calcula el límite de fluencia al aplicar ensayos uniaxiales de tensión para un conjunto de orientaciones cristalográticas dentro del triángulo estándar variando la tasa de deformación y la temperatura de los ensayos. Finalmente, y con el objetivo de ser capaces de predecir una respuesta más dúctil del tungsteno para una variedad de estados de carga, se realizan ensayos biaxiales de tensión sobre algunas de las orientaciones cristalográficas ya estudiadas en función de la temperatura.-------------------------------------------------------------------------ABSTRACT ----------------------------------------------------------Tungsten and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield and low longterm disposal radioactive footprint. However, tungsten also presents a very low fracture toughness, mostly associated with inter-granular failure and bulk plasticity, that limits its applications. As a result of these various and complex conditions of work, the study, development and design of these materials is one of the most important challenges that have emerged in recent years to the scientific community in the field of materials for energy applications. The plastic behavior of body-centered cubic (bcc) refractory metals like tungsten is governed by the kink-pair mediated thermally activated motion of h¿ (\1 11)i screw dislocations on the atomistic scale and by ensembles and interactions of dislocations at larger scales. Modeling this complex behavior requires the application of methods capable of resolving rigorously each relevant scale. The work presented in this thesis proposes a multiscale model approach that gives engineering-level responses to the technical specifications required for the use of tungsten in fusion energy reactors, and it is also supported by the rigorous underlying physics of extensive atomistic simulations. First, the static and dynamic properties of screw dislocations in five interatomic potentials for tungsten are compared, determining which of these ensure greater physical fidelity and computational efficiency. The large strain rates associated with molecular dynamics techniques make the dislocation mobility functions obtained not suitable to be used in the next steps of the multiscale model. Therefore, it is necessary to employ mobility laws obtained from a different method. In this work, we suggest two alternative methods to get the dislocation mobility functions: a kinetic Monte Carlo model and analytical expressions. The set of parameters needed to formulate the kinetic Monte Carlo model and the analytical mobility law are calculated atomistically. These parameters include, but are not limited to: enthalpy and energy barriers of kink-pairs as a function of the stress, width of the kink-pairs, non-Schmid effects ( both twinning-antitwinning asymmetry and non-glide stresses), etc. The function relating dislocation velocity with applied stress and temperature is used as the main source of constitutive information into a dislocation-based crystal plasticity framework. We validate the dependence of the yield strength with the temperature predicted by the model against existing experimental data of tensile tests in singlecrystal tungsten, with excellent agreement between the simulations and the measured data. We then extend the model to a number of crystallographic orientations uniformly distributed in the standard triangle and study the effects of temperature and strain rate. Finally, we perform biaxial tensile tests and provide the yield surface as a function of the temperature for some of the crystallographic orientations explored in the uniaxial tensile tests.

Protein folding kinetics exhibit an Arrhenius temperature dependence when corrected for the temperature dependence of protein stability

The anomalous temperature dependence of protein folding has received considerable attention. Here we show that the temperature dependence of the folding of protein L becomes extremely simple when the effects of temperature on protein stability are corrected for; the logarithm of the folding rate is a linear function of 1/T on constant stability contours in the temperature–denaturant plane. This convincingly demonstrates that the anomalous temperature dependence of folding derives from the temperature dependence of the interactions that stabilize proteins, rather than from the super Arrhenius temperature dependence predicted for the configurational diffusion constant on a rough energy landscape. However, because of the limited temperature range accessible to experiment, the results do not rule out models with higher order temperature dependences. The significance of the slope of the stability-corrected Arrhenius plots is discussed.

Interaction of the synprint site of N-type Ca2+ channels with the C2B domain of synaptotagmin I

N-type Ca2+ channels mediate Ca2+ influx, which initiates fast exocytosis of neurotransmitters at synapses, and they interact directly with the SNARE proteins syntaxin and SNAP-25 (synaptosome-associated protein of 25 kDa) through a synaptic protein interaction (synprint) site in the intracellular loop connecting domains II and III of their α1B subunits. Introduction of peptides containing the synprint site into presynaptic neurons reversibly inhibits synaptic transmission, confirming the importance of interactions with this site in synaptic transmission. Here we report a direct interaction of the synprint peptide from N-type Ca2+ channels with synaptotagmin I, an important Ca2+ sensor for exocytosis, as measured by an affinity-chromatography binding assay and a solid-phase immunoassay. This interaction is mediated by the second C2 domain (C2B) of synaptotagmin I, but is not regulated by Ca2+. Using both immobilized recombinant proteins and native presynaptic membrane proteins, we found that the synprint peptide and synaptotagmin competitively interact with syntaxin. This interaction is Ca2+-dependent because of the Ca2+ dependence of the interactions between syntaxin and these two proteins. These results provide a molecular basis for a physical link between Ca2+ channels and synaptotagmin, and suggest that N-type Ca2+ channels may undergo a complex series of Ca2+-dependent interactions with multiple presynaptic proteins during neurotransmission.

The DNA Helicase Activity of BLM Is Necessary for the Correction of the Genomic Instability of Bloom Syndrome Cells

Bloom syndrome (BS) is a rare autosomal recessive disorder characterized by growth deficiency, immunodeficiency, genomic instability, and the early development of cancers of many types. BLM, the protein encoded by BLM, the gene mutated in BS, is localized in nuclear foci and absent from BS cells. BLM encodes a DNA helicase, and proteins from three missense alleles lack displacement activity. BLM transfected into BS cells reduces the frequency of sister chromatid exchanges and restores BLM in the nucleus. Missense alleles fail to reduce the sister chromatid exchanges in transfected BS cells or restore the normal nuclear pattern. BLM complements a phenotype of a Saccharomyces cerevisiae sgs1 top3 strain, and the missense alleles do not. This work demonstrates the importance of the enzymatic activity of BLM for its function and nuclear localization pattern.

Analysis of the regulatory phosphorylation site in Acanthamoeba myosin IC by using site-directed mutagenesis

The actin-activated ATPase activity of Acanthamoeba myosin IC is stimulated 15- to 20-fold by phosphorylation of Ser-329 in the heavy chain. In most myosins, either glutamate or aspartate occupies this position, which lies within a surface loop that forms part of the actomyosin interface. To investigate the apparent need for a negative charge at this site, we mutated Ser-329 to alanine, asparagine, aspartate, or glutamate and coexpressed the Flag-tagged wild-type or mutant heavy chain and light chain in baculovirus-infected insect cells. Recombinant wild-type myosin IC was indistinguishable from myosin IC purified from Acanthamoeba as determined by (i) the dependence of its actin-activated ATPase activity on heavy-chain phosphorylation, (ii) the unusual triphasic dependence of its ATPase activity on the concentration of F-actin, (iii) its Km for ATP, and (iv) its ability to translocate actin filaments. The Ala and Asn mutants had the same low actin-activated ATPase activity as unphosphorylated wild-type myosin IC. The Glu mutant, like the phosphorylated wild-type protein, was 16-fold more active than unphosphorylated wild type, and the Asp mutant was 8-fold more active. The wild-type and mutant proteins had the same Km for ATP. Unphosphorylated wild-type protein and the Ala and Asn mutants were unable to translocate actin filaments, whereas the Glu mutant translocated filaments at the same velocity, and the Asp mutant at 50% the velocity, as phosphorylated wild-type proteins. These results demonstrate that an acidic amino acid can supply the negative charge in the surface loop required for the actin-dependent activities of Acanthamoeba myosin IC in vitro and indicate that the length of the side chain that delivers this charge is important.

Direct observation of the enhancement of noncooperative protein self-assembly by macromolecular crowding: Indefinite linear self-association of bacterial cell division protein FtsZ

Recent measurements of sedimentation equilibrium and sedimentation velocity have shown that the bacterial cell division protein FtsZ self-associates to form indefinitely long rod-like linear aggregates in the presence of GDP and Mg2+. In the present study, the newly developed technique of non-ideal tracer sedimentation equilibrium was used to measure the effect of high concentrations—up to 150 g/liter—of each of two inert “crowder” proteins, cyanmethemoglobin or BSA, on the thermodynamic activity and state of association of dilute FtsZ under conditions inhibiting (−Mg2+) and promoting (+Mg2+) FtsZ self-association. Analysis of equilibrium gradients of both FtsZ and crowder proteins indicates that, under the conditions of the present experiment, FtsZ interacts with each of the two crowder proteins essentially entirely via steric repulsion, which may be accounted for quantitatively by a simple model in which hemoglobin, albumin, and monomeric FtsZ are modeled as effective spherical hard particles, and each oligomeric species of FtsZ is modeled as an effective hard spherocylinder. The functional dependence of the sedimentation of FtsZ on the concentrations of FtsZ and either crowder indicates that, in the presence of high concentrations of crowder, both the weight-average degree of FtsZ self-association and the range of FtsZ oligomer sizes present in significant abundance are increased substantially.

Differential Expression of the Arabidopsis Nia1 and Nia2 Genes1 : Cytokinin-Induced Nitrate Reductase Activity Is Correlated With Increased Nia1 Transcription and mRNA Levels

Nitrate reductase (NR) activity increased up to 14-fold in response to treatment of Arabidopsis thaliana seedlings with the cytokinin benzyladenine. NR induction was observed in seedlings germinated directly on cytokinin-containing medium, seedlings transferred to cytokinin medium, and seedlings grown in soil in which cytokinin was applied directly to the leaves. About the same level of induction was seen in both wild-type and Nia2-deletion mutants, indicating that increased NR activity is related to the expression of the minor NR gene, Nia1. The steady-state Nia1 mRNA level was increased severalfold in both wild-type and mutant seedlings after benzyladenine treatment. Transcript levels of the Nia2 gene, which is responsible for 90% of the NR activity in developing wild-type seedlings, did not show any changes upon cytokinin treatment. Nuclear run-on assays demonstrated that Nia1 gene transcription increased dramatically after cytokinin treatment.

The relaxation dynamics of the excited electronic states of retinal in bacteriorhodopsin by two-pump-probe femtosecond studies

We present the results of two-pump and probe femtosecond experiments designed to follow the relaxation dynamics of the lowest excited state (S1) populated by different modes. In the first mode, a direct (S0 → S1) radiative excitation of the ground state is used. In the second mode, an indirect excitation is used where the S1 state is populated by the use of two femtosecond laser pulses with different colors and delay times between them. The first pulse excites the S0 → S1 transition whereas the second pulse excites the S1 → Sn transition. The nonradiative relaxation from the Sn state populates the lowest excited state. Our results suggest that the S1 state relaxes faster when populated nonradiatively from the Sn state than when pumped directly by the S0 → S1 excitation. Additionally, the Sn → S1 nonradiative relaxation time is found to change by varying the delay time between the two pump pulses. The observed dependence of the lowest excited state population as well as its dependence on the delay between the two pump pulses are found to fit a kinetic model in which the Sn state populates a different surface (called S′1) than the one being directly excited (S1). The possible involvement of the Ag type states, the J intermediate, and the conical intersection leading to the S0 or to the isomerization product (K intermediate) are discussed in the framework of the proposed model.

Force and kinetic barriers to unzipping of the DNA double helix

A theory of the unzipping of double-stranded DNA is presented and is compared to recent micromanipulation experiments. It is shown that the interactions that stabilize the double helix and the elastic rigidity of single strands simply determine the sequence-dependent ≈12-pN force threshold for DNA strand separation. Using a semimicroscopic model of the binding between nucleotide strands, we show that the greater rigidity of the strands when formed into double-stranded DNA, relative to that of isolated strands, gives rise to a potential barrier to unzipping. The effects of this barrier are derived analytically. The force to keep the extremities of the molecule at a fixed distance, the kinetic rates for strand unpairing at fixed applied force, and the rupture force as a function of loading rate are calculated. The dependence of the kinetics and of the rupture force on molecule length is also analyzed.

Nonuniformity of the constitutive law parameters for shear rupture and quasistatic nucleation to dynamic rupture: a physical model of earthquake generation processes.

Based on the recent high-resolution laboratory experiments on propagating shear rupture, the constitutive law that governs shear rupture processes is discussed in view of the physical principles and constraints, and a specific constitutive law is proposed for shear rupture. It is demonstrated that nonuniform distributions of the constitutive law parameters on the fault are necessary for creating the nucleation process, which consists of two phases: (i) a stable, quasistatic phase, and (ii) the subsequent accelerating phase. Physical models of the breakdown zone and the nucleation zone are presented for shear rupture in the brittle regime. The constitutive law for shear rupture explicitly includes a scaling parameter Dc that enables one to give a common interpretation to both small scale rupture in the laboratory and large scale rupture as earthquake source in the Earth. Both the breakdown zone size Xc and the nucleation zone size L are prescribed and scaled by Dc, which in turn is prescribed by a characteristic length lambda c representing geometrical irregularities of the fault. The models presented here make it possible to understand the earthquake generation process from nucleation to unstable, dynamic rupture propagation in terms of physics. Since the nucleation process itself is an immediate earthquake precursor, deep understanding of the nucleation process in terms of physics is crucial for the short-term (or immediate) earthquake prediction.

Interaction with the recombination hot spot chi in vivo converts the RecBCD enzyme of Escherichia coli into a chi-independent recombinase by inactivation of the RecD subunit.

The RecBCD enzyme of Escherichia coli promotes recombination preferentially at chi nucleotide sequences and has in vivo helicase and strong duplex DNA exonuclease (exoV) activities. The enzyme without the RecD subunit, as in a recD null mutant, promotes recombination efficiently but independently of chi and has no nucleolytic activity. Employing phage lambda red gam crosses, phage T4 2- survival measurements, and exoV assays, it is shown that E. coli cells in which RecBCD has extensive opportunity to interact with linear chi-containing DNA (produced by rolling circle replication of a plasmid with chi or by bleomycin-induced fragmentation of the cellular chromosome) acquire the phenotype of a recD mutant and maintain this for approximately 2 h. It is concluded that RecBCD is converted into RecBC during interaction with chi by irreversible inactivation of RecD. After conversion, the enzyme is released and initiates recombination on other DNA molecules in a chi-independent fashion. Overexpression of recD+ (from a plasmid) prevented the phenotypic change and providing RecD after the change restored chi-stimulated recombination. The observed recA+ dependence of the downregulation of exoV could explain the previously noted "reckless" DNA degradation of recA mutants. It is proposed that chi sites are regulatory elements for the RecBCD to RecBC switch and thereby function as cis- and trans-acting stimulators of RecBC-dependent recombination.

Simulation of the secondary electrons energy deposition produced by proton beams in PMMA: influence of the target electronic excitation description

We have studied the radial dependence of the energy deposition of the secondary electron generated by swift proton beams incident with energies T = 50 keV–5 MeV on poly(methylmethacrylate) (PMMA). Two different approaches have been used to model the electronic excitation spectrum of PMMA through its energy loss function (ELF), namely the extended-Drude ELF and the Mermin ELF. The singly differential cross section and the total cross section for ionization, as well as the average energy of the generated secondary electrons, show sizeable differences at T ⩽ 0.1 MeV when evaluated with these two ELF models. In order to know the radial distribution around the proton track of the energy deposited by the cascade of secondary electrons, a simulation has been performed that follows the motion of the electrons through the target taking into account both the inelastic interactions (via electronic ionizations and excitations as well as electron-phonon and electron trapping by polaron creation) and the elastic interactions. The radial distribution of the energy deposited by the secondary electrons around the proton track shows notable differences between the simulations performed with the extended-Drude ELF or the Mermin ELF, being the former more spread out (and, therefore, less peaked) than the latter. The highest intensity and sharpness of the deposited energy distributions takes place for proton beams incident with T ~ 0.1–1 MeV. We have also studied the influence in the radial distribution of deposited energy of using a full energy distribution of secondary electrons generated by proton impact or using a single value (namely, the average value of the distribution); our results show that differences between both simulations become important for proton energies larger than ~0.1 MeV. The results presented in this work have potential applications in materials science, as well as hadron therapy (due to the use of PMMA as a tissue phantom) in order to properly consider the generation of electrons by proton beams and their subsequent transport and energy deposition through the target in nanometric scales.

Stepwise consolidation of glacigenic sediments of the Prydz Bay

A series of samples from the five sites drilled across the continental shelf and upper slope in Prydz Bay during ODP Leg 119 were consolidation tested in an oedometer. Preconsolidation stresses increase downcore at Sites 739 and 742 in a stepwise manner, and the steps are interpreted to represent periods of increased action of grounded glaciers covering the entire shelf. By the use of theoretical ice sheet surface profiles giving the range of possible ice thicknesses, sediment loading and subsequent erosion seem to be the most important factor for increasing the overconsolidation ratios, and a total glacial erosion exceeding 1 km is possible. Four separate steps in consolidation, here termed "load events" have been identified. The lowermost load event, 1, is correlated to the onset of glaciations reaching the shelf edge and an early period of extensive glaciations, starting in early Oligocene or possibly earlier. Glacial activity related to the buildup of ice in West Antarctica in the late Miocene is tentatively correlated to load event 2. Event 3 is the trace of relatively extensive glacial erosion probably in the Pliocene, whereas the upper step in preconsolidation stress, load event 4, results from the last glaciation reaching the shelf edge, possibly during the late Weichselian. Correlations to other data related to Antarctic glacial history are, however, hampered by the poor age control of the cored diamictites. Consolidation tests may provide a tool for finding the position for hiatuses and unconformities formed subglacially and obscured by subglacial reworking.