A first look at maximally twisted mass lattice QCD calculations at the physical point

In this contribution, a first look at simulations using maximally twisted mass Wilson fermions at the physical point is presented. A lattice action including clover and twisted mass terms is presented and the Monte Carlo histories of one run with two mass-degenerate flavours at a single lattice spacing are shown. Measurements from the light and heavy-light pseudoscalar sectors are compared to previous Nf = 2 results and their phenomenological values. Finally, the strategy for extending simulations to Nf = 2+1+1 is outlined.

Atomic mutagenesis of the ribosomal Sarcin-Ricin-Loop to study EF-G GTPase activation

Translocation factor EF-G, possesses a low basal GTPase activity, which is stimulated by the ribosome. One potential region of the ribosome that triggers GTPase activity of EF-G is the Sarcin-Ricin-Loop (SRL) (helix 95) in domain VI of the 23S rRNA. Structural data showed that the tip of the SRL closely approaches GTP in the active center of EF-G, structural probing data confirmed that EF-G interacts with nucleotides G2655, A2660, G2661 and A2662.1-3 The exocyclic group of adenine at A2660 is required for stimulation of EF-G GTPase activity by the ribosome as demonstrated using atomic mutagenesis.4 Recent crystal structures of EF-G on the ribosome, gave more insights into the molecular mechanism of EF-G GTPase activity.5 Based on the structure of EF-Tu on the ribosome1, the following mechanism of GTPase activation was proposed: upon binding of EF-G to the ribosome, the conserved His92 (E.coli) changes its position, pointing to the γ-phosphate of GTP. In this activated state, the phosphate of residue A2662 of the SRL positions the catalytic His in its active conformation. It was further proposed that the phosphate oxygen of A2662 is involved in a charge-relay system, enabling GTP hydrolysis. In order to test this mechanism, we use the atomic mutagenesis approach, which allows introducing non-natural modifications in the SRL, in the context of the complete 70S ribosome. Therefore, we replaced one of the non-bridging oxygens of A2662 by a methyl group. A methylphosphonat is not able to position or activate a histidine, as it has no free electrons and therefore no proton acceptor function. These modified ribosomes were then tested for stimulation of EF-G GTPase activity. First experiments show that one of the two stereoisomers incorporated into ribosomes does not stimulate GTPase activity of EF-G, whereas the other is active. From this we conclude that indeed the non-bridging phosphate oxygen of A2662 is involved in EF-G GTPase activation by the ribosome. Ongoing experiments aim at revealing the contribution of this non-bridging oxygen at A2662 to the mechanism of EF-G GTPase activation at the atomic level.

Atomic mutagenesis of the ribosomal Sarcin-Ricin-Loop to study EF-G GTPase activation

Translocation factor EF-G, possesses a low basal GTPase activity, which is stimulated by the ribosome. One potential region of the ribosome that triggers GTPase activity of EF-G is the Sarcin-Ricin-Loop (SRL) (helix 95) in domain VI of the 23S rRNA. Structural data showed that the tip of the SRL closely approaches GTP in the active center of EF-G, structural probing data confirmed that EF-G interacts with nucleotides G2655, A2660, G2661 and A2662.1-3 The exocyclic group of adenine at A2660 is required for stimulation of EF-G GTPase activity by the ribosome as demonstrated using atomic mutagenesis.4 Recent crystal structures of EF-G on the ribosome, gave more insights into the molecular mechanism of EF-G GTPase activity.5 Based on the structure of EF-Tu on the ribosome1, the following mechanism of GTPase activation was proposed: upon binding of EF-G to the ribosome, the conserved His92 (E.coli) changes its position, pointing to the γ-phosphate of GTP. In this activated state, the phosphate of residue A2662 of the SRL positions the catalytic His in its active conformation. It was further proposed that the phosphate oxygen of A2662 is involved in a charge-relay system, enabling GTP hydrolysis. In order to test this mechanism, we use the atomic mutagenesis approach, which allows introducing non-natural modifications in the SRL, in the context of the complete 70S ribosome. Therefore, we replaced one of the non-bridging oxygens of A2662 by a methyl group. A methylphosphonat is not able to position or activate a histidine, as it has no free electrons and therefore no proton acceptor function. These modified ribosomes were then tested for stimulation of EF-G GTPase activity. First experiments show that one of the two stereoisomers incorporated into ribosomes does not stimulate GTPase activity of EF-G, whereas the other is active. From this we conclude that indeed the non-bridging phosphate oxygen of A2662 is involved in EF-G GTPase activation by the ribosome. Ongoing experiments aim at revealing the contribution of this non-bridging oxygen at A2662 to the mechanism of EF-G GTPase activation at the atomic level.

Accurate rotational constant and bond lengths of hexafluorobenzene by femtosecond rotational Raman coherence spectroscopy and ab initio calculations

The gas-phase rotational motion of hexafluorobenzene has been measured in real time using femtosecond (fs) time-resolved rotational Raman coherence spectroscopy (RR-RCS) at T = 100 and 295 K. This four-wave mixing method allows to probe the rotation of non-polar gas-phase molecules with fs time resolution over times up to ∼5 ns. The ground state rotational constant of hexafluorobenzene is determined as B 0 = 1029.740(28) MHz (2σ uncertainty) from RR-RCS transients measured in a pulsed seeded supersonic jet, where essentially only the v = 0 state is populated. Using this B 0 value, RR-RCS measurements in a room temperature gas cell give the rotational constants B v of the five lowest-lying thermally populated vibrationally excited states ν7/8, ν9, ν11/12, ν13, and ν14/15. Their B v constants differ from B 0 by between −1.02 MHz and +2.23 MHz. Combining the B 0 with the results of all-electron coupled-cluster CCSD(T) calculations of Demaison et al. [Mol. Phys.111, 1539 (2013)] and of our own allow to determine the C-C and C-F semi-experimental equilibrium bond lengths r e(C-C) = 1.3866(3) Å and r e(C-F) = 1.3244(4) Å. These agree with the CCSD(T)/wCVQZ r e bond lengths calculated by Demaison et al. within ±0.0005 Å. We also calculate the semi-experimental thermally averaged bond lengths r g(C-C)=1.3907(3) Å and r g(C-F)=1.3250(4) Å. These are at least ten times more accurate than two sets of experimental gas-phase electron diffraction r g bond lengths measured in the 1960s.

Specificities of Caenorhabditis elegans and human hairpin binding proteins for the first nucleotide in the histone mRNA hairpin loop.

The 3' ends of animal replication-dependent histone mRNAs are formed by endonucleolytic cleavage of the primary transcripts downstream of a highly conserved RNA hairpin. The hairpin-binding protein (HBP) binds to this RNA element and is involved in histone RNA 3' processing. A minimal RNA-binding domain (RBD) of approximately 73 amino acids that has no similarity with other known RNA-binding motifs was identified in human HBP [Wang Z-F et al., Genes & Dev, 1996, 10:3028-3040]. The primary sequence identity between human and Caenorhabditis elegans RBDs is 55% compared to 38% for the full-length proteins. We analyzed whether differences between C. elegans and human HBP and hairpins are reflected in the specificity of RNA binding. The C. elegans HBP and its RBD recognize only their cognate RNA hairpins, whereas the human HBP or RBD can bind both the mammalian and the C. elegans hairpins. This selectivity of C. elegans HBP is mostly mediated by the first nucleotide in the loop, which is C in C. elegans and U in all other metazoans. By converting amino acids in the human RBD to the corresponding C. elegans residues at places where the latter deviates from the consensus, we could identify two amino acid segments that contribute to selectivity for the first nucleotide of the hairpin loop.

All-loop anomalous dimensions in integrable λ-deformed σ-models

We calculate the all-loop anomalous dimensions of current operators in λ-deformed σ-models. For the isotropic integrable deformation and for a semi-simple group G we compute the anomalous dimensions using two different methods. In the first we use the all-loop effective action and in the second we employ perturbation theory along with the Callan–Symanzik equation and in conjunction with a duality-type symmetry shared by these models. Furthermore, using CFT techniques we compute the all-loop anomalous dimension of bilinear currents for the isotropic deformation case and a general G . Finally we work out the anomalous dimension matrix for the cases of anisotropic SU(2) and the two couplings, corresponding to the symmetric coset G/H and a subgroup H, splitting of a group G.

Excitonic splittings in molecular dimers: why static ab initio calculations cannot match them

After decades of research on molecular excitons, only few molecular dimers are available on which exciton and vibronic coupling theories can be rigorously tested. In centrosymmetric H-bonded dimers consisting of identical (hetero)aromatic chromophores, the monomer electronic transition dipole moment vectors subtract or add, yielding S0 → S1 and S0 → S2 transitions that are symmetry-forbidden or -allowed, respectively. Symmetry breaking by 12C/13C or H/D isotopic substitution renders the forbidden transition weakly allowed. The excitonic coupling (Davydov splitting) can then be measured between the S0 → S1 and S0 → S2 vibrationless bands. We discuss the mass-specific excitonic spectra of five H-bonded dimers that are supersonically cooled to a few K and investigated using two-color resonant two-photon ionization spectroscopy. The excitonic splittings Δcalc predicted by ab initio methods are 5–25 times larger than the experimental excitonic splittings Δexp. The purely electronic ab initio splittings need to be reduced (“quenched”), reflecting the coupling of the electronic transition to the optically active vibrations of the monomers. The so-called quenching factors Γ < 1 can be determined from experiment (Γexp) and/or calculation (Γcalc). The vibronically quenched splittings Γ·Δcalc are found to nicely reproduce the experimental exciton splittings.

Paleo-sea surface temperature calculations in the equatorial east Atlantic

We present two ~270 kyr paleo-sea surface temperature (SST) records from the Equatorial Divergence and the South Equatorial Current derived from Mg/Ca ratios in the planktic foraminifer Globigerinoides sacculifer. The present study suggests that the magnesium signature of G. sacculifer provides a seasonal SST estimate from the upper ~50 m of the water column generated during upwelling in austral low-latitude fall/winter. Common to both down-core records is a glacial-interglacial amplitude of ~3°-3.5°C for the last climatic changes and lower Holocene and glacial oxygen isotope stage 2 temperatures compared with interglacial stage 5.5 and glacial stage 6 temperatures, respectively. The comparison to published SST estimates from alkenones, oxygen isotopes, and foraminiferal transfer function from the same core material pinpoints discrepancies and conformities between methods.

Propagation of nuclear data uncertainties in fuel cycle calculations using MONTE-CARLO Technique

The uncertainty propagation in fuel cycle calculations due to Nuclear Data (ND) is a important important issue for : issue for : • Present fuel cycles (e.g. high burnup fuel programme) • New fuel cycles designs (e.g. fast breeder reactors and ADS) Different error propagation techniques can be used: • Sensitivity analysis • Response Response Surface Method Surface Method • Monte Carlo technique Then, p p , , in this paper, it is assessed the imp y pact of ND uncertainties on the decay heat and radiotoxicity in two applications: • Fission Pulse Decay ( y Heat calculation (FPDH) • Conceptual design of European Facility for Industrial Transmutation (EFIT)

The PSE-3D instability analysis methodology for flows depending strongly on two and weakly on the third spatial dimension

The present contribution discusses the development of a PSE-3D instability analysis algorithm, in which a matrix forming and storing approach is followed. Alternatively to the typically used in stability calculations spectral methods, new stable high-order finitedifference-based numerical schemes for spatial discretization 1 are employed. Attention is paid to the issue of efficiency, which is critical for the success of the overall algorithm. To this end, use is made of a parallelizable sparse matrix linear algebra package which takes advantage of the sparsity offered by the finite-difference scheme and, as expected, is shown to perform substantially more efficiently than when spectral collocation methods are used. The building blocks of the algorithm have been implemented and extensively validated, focusing on classic PSE analysis of instability on the flow-plate boundary layer, temporal and spatial BiGlobal EVP solutions (the latter necessary for the initialization of the PSE-3D), as well as standard PSE in a cylindrical coordinates using the nonparallel Batchelor vortex basic flow model, such that comparisons between PSE and PSE-3D be possible; excellent agreement is shown in all aforementioned comparisons. Finally, the linear PSE-3D instability analysis is applied to a fully three-dimensional flow composed of a counter-rotating pair of nonparallel Batchelor vortices.

Proposal of an Alternative Model for Speed-Flow Relationship in Two-Lane Highways

A research has been carried out in two-lanehighways in the Madrid Region to propose an alternativemodel for the speed-flowrelationship using regular loop data. The model is different in shape and, in some cases, slopes with respect to the contents of Highway Capacity Manual (HCM). A model is proposed for a mountainous area road, something for which the HCM does not provide explicitly a solution. The problem of a mountain road with high flows to access a popular recreational area is discussed, and some solutions are proposed. Up to 7 one-way sections of two-lanehighways have been selected, aiming at covering a significant number of different characteristics, to verify the proposed method the different classes of highways on which the Manual classifies them. In order to enunciate the model and to verify the basic variables of these types of roads a high number of data have been used. The counts were collected in the same way that the Madrid Region Highway Agency performs their counts. A total of 1.471 hours have been collected, in periods of 5 minutes. The models have been verified by means of specific statistical test (R2, T-Student, Durbin-Watson, ANOVA, etc.) and with the diagnostics of the contrast of assumptions (normality, linearity, homoscedasticity and independence). The model proposed for this type of highways with base conditions, can explain the different behaviors as traffic volumes increase, and follows a polynomial multiple regression model of order 3, S shaped. As secondary results of this research, the levels of service and the capacities of this road have been measured with the 2000 HCM methodology, and the results discussed. © 2011 Published by Elsevier Ltd.

Luminescence from two-dimensional electron gases in InAlN/GaN heterostructures with different In content

The luminescence properties of InxAl1−xN/GaN heterostructures are investigated systematically as a function of the In content (x = 0.067 − 0.208). The recombination between electrons confined in the two-dimensional electron gas and free holes in the GaN template is identified and analyzed. We find a systematic shift of the recombination with increasing In content from about 80 meV to only few meV below the GaN exciton emission. These results are compared with model calculations and can be attributed to the changing band profile and originating from the polarization gradient between InAlN and GaN.

Revisiting the method to obtain the mechanical properties of hydrided fuel cladding in the hoop direction

The method reported in the literature to calculate the stress–strain curve of nuclear fuel cladding from ring tensile test is revisited in this paper and a new alternative is presented. In the former method, two universal curves are introduced under the assumption of small strain. In this paper it is shown that these curves are not universal, but material-dependent if geometric nonlinearity is taken into account. The new method is valid beyond small strains, takes geometric nonlinearity into consideration and does not need universal curves. The stress–strain curves in the hoop direction are determined by combining numerical calculations with experimental results in a convergent loop. To this end, ring tensile tests were performed in unirradiated hydrogen-charged samples. The agreement among the simulations and the experimental results is excellent for the range of concentrations tested (up to 2000 wppm hydrogen). The calculated stress–strain curves show that the mechanical properties do not depend strongly on the hydrogen concentration, and that no noticeable strain hardening occurs. However, ductility decreases with the hydrogen concentration, especially beyond 500 wppm hydrogen. The fractographic results indicate that as-received samples fail in a ductile fashion, whereas quasicleavage is bserved in the hydrogen-charged samples.

Instabilities in two-liquid layers subject to a horizontal temperature gradient

We study theoretically the stability of two superposed fluid layers heated laterally. The fluids are supposed to be immiscible, the interface undeformable and of infinite horizontal extension. Combined thermocapillary and buoyancy forces give rise to a basic flow when a temperature difference is applied. The calculations are performed for a melt of GaAs under a layer of molten B2 O3 , a configuration of considerable technological importance. Four dif- ferent flow patterns and five temperature configurations are found for the basic state in this system. A linear stability analysis shows that the basic state may be destabilized by oscilla- tory motions leading to the so-called hydrothermal waves. Depending on the relative height of the two layers these hydrothermal waves propagate parallel or perpendicular to the temperature gradient. This analysis reveals that these perturbations can alter significantly the liquid flow in the liquid-encapsulated crystal growth techniques.

Ab-Initio calculations including Van der Waals interactions: the SnS2 layered material

Tin disulfide SnS2 was recently proposed as a high efficiency solar cell precursor [1]. The aim of this work is a deep study of the structural disposition of the most important polytipes of this layered material, not only describing the electronic correlation but also the interatomic Van der Waals interactions that is present between the layers. The two recent implementations to take Van der Waals interactions into account in the VASP code are the self-consistent Dion et al. [2] functional optimized for solids by Michaelides et al [3] and the Grimme [4] dispersion correction that is applied after each autoconsistent PBE electronic calculation. In this work these two methods are compared with DFT PBE functional. The results we will presented at this Conference, demonstrates the enhancement of the geometric parameters by the use of the Van der Waals interactions in agreement with the experimental values.