Molecular mechanisms of Ca$\sp{2+}$ signal transduction from cardiac troponin C to troponin I

$\rm Ca\sp{2+}$-dependent exposure of an N-terminal hydrophobic region in troponin C (TnC) is thought to be important for the regulation of contraction in striated muscle. To study these conformational changes in cardiac troponin (cTnC), the $\varepsilon$C and $\varepsilon$H chemical shifts for all 10 Met residues in cTnC were sequence-specific assigned on NMR spectra using a combination of two dimensional NMR techniques and site-directed mutagenesis. The assigned methyl-Met chemical shifts were used as structural markers to monitor conformational changes induced by $\rm Ca\sp{2+}.$ The results showed that binding of $\rm Ca\sp{2+}$ to the regulatory site in the N-domain induced large changes in the $\varepsilon$H and $\varepsilon$C chemical shifts of Met 45, Met 80, Met 81 in the predicted N-terminal hydrophobic region, but had no effect on the chemical shifts of Met residues located in the C-domain. These results suggest that the $\rm Ca\sp{2+}$-dependent functions of cTnC are mainly through N-terminal domain of cTnC.^ To further define the molecular mechanism by which TnC regulates muscle contraction, single Cys residues were engineered at positions 45, 81, 84 or 85 in the N-terminal hydrophobic region of cTnC to provide sites for attachment of specific blocking groups. Blocking groups were coupled to these Cys residues in cTnC mutants and the covalent adducts were tested for activity in TnC-extracted myofibrils. Covalent modification of cTnC(C45) had no effect on maximal myofibril ATPase activity. Greatly decreased myofibril ATPase activity resulted when the peptide or biotin was conjugated to residue 81 in cTnC(C81), while less inhibition resulted from covalent modification of cTnC(C84) or cTnC(C85). The results suggest that limited sites of the N-terminal hydrophobic region in cTnC are important for transducing the $\rm Ca\sp{2+}$ signal to troponin I (TnI) and are sensitive to modification, while other regions are less important or can adapt to steric hindrances introduced by bulky blocking groups.^ Although the exposed TnI interaction site in the N-terminal hydrophobic region of TnC is crucial for function of TnC, other regions in the N-domain of TnC may also participate in transducing the $\rm Ca\sp{2+}$ signal and conferring the maximal activation of actomyosin ATPase. The interactions between the B-/C-helices of cTnC and cTnI were characterized using a combination of site-directed mutagenesis, fluorescence and covalent modification. The results suggest that the $\rm Ca\sp{2+}$-dependent interactions of the B-/C-helices of cTnC with TnI may be required for the maximal activation of muscle contraction. ^

Bacteriological diagnosis and molecular strain typing of Mycobacterium bovis and Mycobacterium caprae

The primary isolation of a Mycobacterium sp. of the Mycobacterium tuberculosis complex from an infected animal provides a definitive diagnosis of tuberculosis. However, as Mycobacterium bovis and Mycobacterium caprae are difficult to isolate, particularly for animals in the early stages of disease, success is dependent on the optimal performance of all aspects of the bacteriological process, from the initial choice of tissue samples at post-mortem examination or clinical samples, to the type of media and conditions used to cultivate the microorganism. Each step has its own performance characteristics, which can contribute to sensitivity and specificity of the procedure, and may need to be optimized in order to achieve the gold standard diagnosis. Having isolated the slow-growing mycobacteria, species identification and fine resolution strain typing are keys to understanding the epidemiology of the disease and to devise strategies to limit transmission of infection. New technologies have emerged that can now even discriminate different isolates from the same animal. In this review we highlight the key factors that contribute to the accuracy of bacteriological diagnosis of M. bovis and M. caprae, and describe the development of advanced genotyping techniques that are increasingly used in diagnostic laboratories for the purpose of supporting detailed epidemiological investigations.

Gas-phase Lifetimes of Nucleobase Analogues by Picosecond Pumpionization and Streak Techniques

The picosecond (ps) timescale is relevant for the investigation of many molecular dynamical processes such as fluorescence, nonradiative relaxation, intramolecular vibrational relaxation, molecular rotation and intermolecular energy transfer, to name a few. While investigations of ultrafast (femtosecond) processes of biological molecules, e.g. nucleobases and their analogues in the gas phase are available, there are few investigations on the ps time scale. We have constructed a ps pump-ionization setup and a ps streak camera fluorescence apparatus for the determination of lifetimes of supersonic jet-cooled and isolated molecules and clusters. The ps pump-ionization setup was used to determine the lifetimes of the nucleobase analogue 2-aminopurine (2AP) and of two 2AP˙(H2O)n water cluster isomers with n=1 and 2. Their lifetimes lie between 150 ps and 3 ns and are strongly cluster-size dependent. The ps streak camera setup was used to determine accurate fluorescence lifetimes of the uracil analogue 2-pyridone (2PY), its self-dimer (2PY)2, two isomers of its trimer (2PY)3 and its tetramer (2PY)4, which lie in the 7–12 ns range.

The Exploitation of Versatile Building Blocks for the Self-Assembly of Novel Molecular Magnets

Using molecular building blocks to self-assemble lattices supporting long-range magnetic order is currently an active area of solid-state chemistry. Consequently, it is the realm of supramolecular chemistry that synthetic chemists are turning to in order to develop techniques for the synthesis of structurally well-defined supramolecular materials. In recent years we have investigated the versatility and usefulness of two classes of molecular building blocks, namely, tris-oxalato transition-metal (M. Pilkington and S. Decurtins, in “Magnetoscience—From Molecules to Materials,” Wiley–VCH, 2000), and octacyanometalate complexes (Pilkington and Decurtins, Chimia 54, 593 (2001)), for applications in the field of molecule-based magnets. Anionic, tris-chelated oxalato building blocks are able to build up two-dimensional honeycomb-layered structural motifs as well as three-dimensional decagon frameworks. The discrimination between the crystallization of the two- or three-dimensional structures relies on the choice of the templating counterions (Decurtins, Chimia 52, 539 (1998); Decurtins et al. Mol. Cryst. Liq. Cryst. 273, 167 (1995); New J. Chem. 117 (1998)). These structural types display a range of ferro, ferri, and antiferromagnetic properties (Pilkington and Decurtins, in “Magnetoscience—From Molecules to Materials”). Octacyanometalate building blocks self-assemble to afford two new classes of cyano-bridged compounds namely, molecular clusters and extended three dimensional networks (J. Larionova et al., Angew. Chem. Int. Ed. 39, 1605 (2000); Pilkington et al., in preparation). The molecular cluster with a MnII9MoV6 core has the highest ground state spin value, S=51/2, reported to-date (Larionova et al., Angew. Chem. Int. Ed. 39, 1605 (2000)). In the high-temperature regime, the magnetic properties are characterized by ferromagnetic intracluster coupling. In the magnetic range below 44 K, the magnetic cluster signature is lost as possibly a bulk behavior starts to emerge. The three-dimensional networks exhibit both paramagnetic and ferromagnetic behavior, since the magnetic properties of these materials directly reflect the electronic configuration of the metal ion incorporated into the octacyanometalate building blocks (Pilkington et al., in preparation). For both the oxalate- and cyanide-bridged materials, we are able to manipulate the magnetic properties of the supramolecular assemblies by tuning the electronic configurations of the metal ions incorporated into the appropriate molecular building blocks (Pilkington and Decurtins, in “Magnetoscience—From Molecules to Materials,” Chimia 54, 593 (2000)).

Seasonal variation of high-molecular-weight compounds in the water-soluble fraction of organic urban aerosols

Aerosol samples were collected in Zurich, Switzerland, at an urban background site and were analyzed with size exclusion chromatography (SEC) and laser/desorption ionization mass spectrometry (LDI-MS) for water-soluble organic compounds with high molecular weight. Daily samples were collected during two campaigns in winter and summer, for 1 month each. The concentration of high-molecular-weight compounds (humic-like substances (HULIS)) was between 0.4 and 4 μg/m3 in winter and summer. The most intense signals in the LDI-MS mass spectra were measured between m/z150 and 500, comparing well with the mode of the two main high mass peaks determined with SEC corresponding to masses between 200 and 600 Da. For the maximum molecular weight, however, different results were obtained by the two techniques: whereas a maximum molecular weight between 1300 and 3300 Da was found with SEC, hardly any peaks above m/z700 were measured with LDI-MS. During summer the maximum molecular weight of HULIS (determined with SEC) correlates positively with several parameters such as ozone and increased temperature indicative of enhanced atmospheric photo-oxidation. The HULIS concentration also correlates positively with the oxalic acid concentration in the particles. This suggests that HULIS are generated by secondary processes in summer. The lack of such correlations during winter suggests that other sources and processes might be important during colder seasons.

Molecular control of cardiac sodium homeostasis in health and disease.

INTRODUCTION Cardiac myocytes utilize three high-capacity Na transport processes whose precise function can determine myocyte fate and the triggering of arrhythmias in pathological settings. We present recent results on the regulation of all three transporters that may be important for an understanding of cardiac function during ischemia/reperfusion episodes. METHODS AND RESULTS Refined ion selective electrode (ISE) techniques and giant patch methods were used to analyze the function of cardiac Na/K pumps, Na/Ca exchange (NCX1), and Na/H exchange (NHE1) in excised cardiac patches and intact myocytes. To consider results cohesively, simulations were developed that account for electroneutrality of the cytoplasm, ion homeostasis, water homeostasis (i.e., cell volume), and cytoplasmic pH. The Na/K pump determines the average life-time of Na ions (3-10 minutes) as well as K ions (>30 minutes) in the cytoplasm. The long time course of K homeostasis can determine the time course of myocyte volume changes after ion homeostasis is perturbed. In excised patches, cardiac Na/K pumps turn on slowly (-30 seconds) with millimolar ATP dependence, when activated for the first time. In steady state, however, pumps are fully active with <0.2 mM ATP and are nearly unaffected by high ADP (2 mM) and Pi (10 mM) concentrations as may occur in ischemia. NCX1s appear to operate with slippage that contributes to background Na influx and inward current in heart. Thus, myocyte Na levels may be regulated by the inactivation reactions of the exchanger which are both Na- and proton-dependent. NHE1 also undergo strong Na-dependent inactivation, whereby a brief rise of cytoplasmic Na can cause inactivation that persists for many minutes after cytoplasmic Na is removed. This mechanism is blocked by pertussis toxin, suggesting involvement of a Na-dependent G-protein. Given that maximal NCX1- and NHE1-mediated ion fluxes are much greater than maximal Na/K pump-mediated Na extrusion in myocytes, the Na-dependent inactivation mechanisms of NCX1 and NHE1 may be important determinants of cardiac Na homeostasis. CONCLUSIONS Na/K pumps appear to be optimized to continue operation when energy reserves are compromised. Both NCX1 and NHE1 activities are regulated by accumulation of cytoplasmic Na. These principles may importantly control cardiac cytoplasmic Na and promote myocyte survival during ischemia/reperfusion episodes by preventing Ca overload.

First international external quality assessment of molecular diagnostics for Mers-CoV

BACKGROUND: Since the discovery of Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012, diagnostic protocols were quickly published and deployed globally. OBJECTIVES: We set out to assess the quality of MERS-CoV molecular diagnostics worldwide. STUDY DESIGN: Both sensitivity and specificity were assessed using 12 samples containing different viral loads of MERS-CoV or common coronaviruses (OC43, 229E, NL63, HKU1). RESULTS: The panel was sent to more than 106 participants, of which 99 laboratories from 6 continents returned 189 panel results.Scores ranged from 100% (84 laboratories) to 33% (1 laboratory). 15% of respondents reported quantitative results, 61% semi-quantitative (Ct-values or time to positivity) and 24% reported qualitative results. The major specific technique used was real-time RT-PCR using the WHO recommended targets upE, ORF1a and ORF1b. The evaluation confirmed that RT-PCRs targeting the ORF1b are less sensitive, and therefore not advised for primary diagnostics. CONCLUSIONS: The first external quality assessment MERS-CoV panel gives a good insight in molecular diagnostic techniques and their performances for sensitive and specific detection of MERS-CoV RNA globally. Overall, all laboratories were capable of detecting MERS-CoV with some differences in sensitivity. The observation that 8% of laboratories reported false MERS-CoV positive single assay results shows room for improvement, and the importance of using confirmatory targets.

STUDIES OF THE MOLECULAR MECHANISMS OF CHROMOSOME ABERRATION FORMATION (DNA REPAIR, CROSSLINKS, MITOMYCIN C)

The objective of this research has been to study the molecular basis for chromosome aberration formation. Predicated on a variety of data, Mitomycin C (MMC)-induced DNA damage has been postulated to cause the formation of chromatid breaks (and gaps) by preventing the replication of regions of the genome prior to mitosis. The basic protocol for these experiments involved treating synchronized Hela cells in G(,1)-phase with a 1 (mu)g/ml dose of MMC for one hour. After removing the drug, cells were then allowed to progress to mitosis and were harvested for analysis by selective detachment. Utilizing the alkaline elution assay for DNA damage, evidence was obtained to support the conclusion that Hela cells can progress through S-phase into mitosis with intact DNA-DNA interstrand crosslinks. A higher level of crosslinking was observed in those cells remaining in interphase compared to those able to reach mitosis at the time of analysis. Dual radioisotope labeling experiments revealed that, at this dose, these crosslinks were associated to the same extent with both parental and newly replicated DNA. This finding was shown not to be the result of a two-step crosslink formation mechanism in which crosslink levels increase with time after drug treatment. It was also shown not to be an artefact of the double-labeling protocol. Using neutral CsCl density gradient ultracentrifugation of mitotic cells containing BrdU-labeled newly replicated DNA, control cells exhibited one major peak at a heavy/light density. However, MMC-treated cells had this same major peak at the heavy/light density, in addition to another minor peak at a density characteristic for light/light DNA. This was interpreted as indicating either: (1) that some parental DNA had not been replicated in the MMC treated sample or; (2) that a recombination repair mechanism was operational. To distinguish between these two possibilities, flow cytometric DNA fluorescence (i.e., DNA content) measurements of MMC-treated and control cells were made. These studies revealed that the mitotic cells that had been treated with MMC while in G(,1)-phase displayed a 10-20% lower DNA content than untreated control cells when measured under conditions that neutralize chromosome condensation effects (i.e., hypotonic treatment). These measurements were made under conditions in which the binding of the drug, MMC, was shown not to interfere with the stoichiometry of the ethidium bromide-mithramycin stain. At the chromosome level, differential staining techniques were used in an attempt to visualize unreplicated regions of the genome, but staining indicative of large unreplicated regions was not observed. These results are best explained by a recombinogenic mechanism. A model consistent with these results has been proposed.^

Characterization of fecal contamination targeting 16S rRNA bacteroides molecular markers in the Santa Cruz River, Arizona

Understanding the origins, transport and fate of contamination is essential to effective management of water resources and public health. Individuals and organizations with management responsibilities need to understand the risks to ecosystems and to humans from contact with contamination. Managers also need to understand how key contaminants vary over time and space in order to design and prioritize mitigation strategies. Tumacacori National Historic Park (NHP) is responsible for management of its water resources for the benefit of the park and for the health of its visitors. The existence of microbial contaminants in the park poses risks that must be considered in park planning and operations. The water quality laboratory at the Maricopa Agricultural Center (in collaboration with stakeholder groups and individuals located in the ADEQ-targeted watersheds) identified biological changes in surface water quality in impaired reaches rivers to determine the sources of Escherichia coli (E. coli); bacteria utilizing innovative water quality microbial/bacterial source tracking methods. The end goal was to support targeted watershed groups and ADEQ towards E. coli reductions. In the field monitoring was conducted by the selected targeted watershed groups in conjunction with The University of Arizona Maricopa Agricultural Center Water Quality Laboratory. This consisted of collecting samples for Bacteroides testing from multiple locations on select impaired reaches, to determine contamination resulting from cattle, human recreation, and other contributions. Such testing was performed in conjunction with high flow and base flow conditions in order to accurately portray water quality conditions and variations. Microbial monitoring was conducted by The University of Arizona Water Quality Laboratory at the Maricopa Agricultural Center using genetic typing to differentiate among two categories of Bacteroides: human and all (total). Testing used microbial detection methodologies and molecular source tracking techniques.^

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Molecular properties of gases of ODP Leg 204 sites

The recognition of finely disseminated gas hydrate in deep marine sediments heavily depends on various indirect techniques because this mineral quickly decomposes upon recovery from in situ pressure and temperature conditions. Here, we discuss molecular properties of closely spaced gas voids (formed as a result of core recovery) and gas hydrates from an area of relatively low gas flux at the flanks of the southern Hydrate Ridge offshore Oregon (ODP Sites 1244, 1245 and 1247). Within the gas hydrate occurrence zone (GHOZ), the concentration of ethane (C2) and propane (C3) in adjacent gas voids shows large variability. Sampled gas hydrates are enriched in C2 relative to void gases but do not contain C3. We suggest that the observed variations in the composition of void gases is a result of molecular fractionation during crystallization of structure I gas hydrate that contains C2 but excludes C3 from its crystal lattice. This hypothesis is used to identify discrete intervals of finely disseminated gas hydrate in cored sediments. Variations in gas composition help better constrain gas hydrate distribution near the top of the GHOZ along with variations in pore water chemistry and core temperature. Sediments near the base of the gas hydrate stability zone are relatively enriched in C2+ hydrocarbon gases. Complex and poorly understood geological and geochemical processes in these deeper sediments make the identification of gas hydrate based on molecular properties of void gases more ambiguous. The proposed technique appears to be a useful tool to better understand the distribution of gas hydrate in marine sediments and ultimately the role of gas hydrate in the global carbon cycle.

A molecular dynamics study of swift heavy ion irradiation of amorphous silica: the role of thermal effects

Irradiation with swift heavy ions (SHI), roughly defined as those having atomic masses larger than 15 and energies exceeding 1 MeV/amu, may lead to significant modification of the irradiated material in a nanometric region around the (straight) ion trajectory (latent tracks). In the case of amorphous silica, SHI irradiation originates nano-tracks of higher density than the virgin material (densification). As a result, the refractive index is increased with respect to that of the surroundings. Moreover, track overlapping leads to continuous amorphous layers that present a significant contrast with respect to the pristine substrate. We have recently demonstrated that SHI irradiation produces a large number of point defects, easily detectable by a number of experimental techniques (work presented in the parallel conference ICDIM). The mechanisms of energy transfer from SHI to the target material have their origin in the high electronic excitation induced in the solid. A number of phenomenological approaches have been employed to describe these mechanisms: coulomb explosion, thermal spike, non-radiative exciton decay, bond weakening. However, a detailed microscopic description is missing due to the difficulty of modeling the time evolution of the electronic excitation. In this work we have employed molecular dynamics (MD) calculations to determine whether the irradiation effects are related to the thermal phenomena described by MD (in the ps domain) or to electronic phenomena (sub-ps domain), e.g., exciton localization. We have carried out simulations of up to 100 ps with large boxes (30x30x8 nm3) using a home-modified version of MDCASK that allows us to define a central hot cylinder (ion track) from which heat flows to the surrounding cold bath (unirradiated sample). We observed that once the cylinder has cooled down, the Si and O coordination numbers are 4 and 2, respectively, as in virgin silica. On the other hand, the density of the (cold) cylinder increases with respect to that of silica and, furthermore, the silica network ring size decreases. Both effects are in agreement with the observed densification. In conclusion, purely thermal effects do not explain the generation of point defects upon irradiation, but they do account for the silica densification.

Preliminary approach to neutron instrument selecction at ESSB-Bilbao based on experience at ISIS molecular spectroscopy group

Collaborative efforts between the Neutronics and Target Design Group at the Instituto de Fusión Nuclear and the Molecular Spectroscopy Group at the ISIS Pulsed Neutron and Muon Source date back to 2012 in the context of the ESS-Bilbao project. The rationale for these joint activities was twofold, namely: to assess the realm of applicability of the low-energy neutron source proposed by ESS-Bilbao - for details; and to explore instrument capabilities for pulsed-neutron techniques in the range 0.05-3 ms, a time range where ESS-Bilbao and ISIS could offer a significant degree of synergy and complementarity. As part of this collaboration, J.P. de Vicente has spent a three-month period within the ISIS Molecular Spectroscopy Group, to gain hands-on experience on the practical aspects of neutron-instrument design and the requisite neutron-transport simulations. To date, these activities have resulted in a joint MEng thesis as well as a number of publications and contributions to national and international conferences. Building upon these previous works, the primary aim of this report is to provide a self-contained discussion of general criteria for instrument selection at ESS-Bilbao, the first accelerator-driven, low-energy neutron source designed in Spain. To this end, Chapter 1 provides a brief overview of the current design parameters of the accelerator and target station. Neutron moderation is covered in Chapter 2, where we take a closer look at two possible target-moderator-reflector configurations and pay special attention to the spectral and temporal characteristics of the resulting neutron pulses. This discussion provides a necessary starting point to assess the operation of ESSB in short- and long-pulse modes. These considerations are further explored in Chapter 3, dealing with the primary characteristics of ESS-Bilbao as a short- or long-pulse facility in terms of accessible dynamic range and spectral resolution. Other practical aspects including background suppression and the use of fast choppers are also discussed. The guiding principles introduced in the first three chapters are put to use in Chapter 4 where we analyse in some detail the capabilities of a small-angle scattering instrument, as well as how specific scientific requirements can be mapped onto the optimal use of ESS-Bilbao for condensed-matter research. Part 2 of the report contains additional supporting documentation, including a description of the ESSB McStas component, a detailed characterisation of moderator response and neutron pulses, and estimates ofparameters associated with the design and operation of neutron choppers. In closing this brief foreword, we wish to thank both ESS-Bilbao and ISIS for their continuing encouragement and support along the way.

Molecular Dynamics Simulations of Couette flow

The first steps towards developing a continuum-molecular coupled simulations techniques are presented, for the purpose of computing macroscopic systems of confined fluids. The idea is to compute the interface wall-fluid by Molecular Dynamics simulations, where Lennard-Jones potential (and others) have been employed for the molecular interactions, so the usual non slip boundary condition is not specified. Instead, a shear rate can be imposed at the wall, which allows to obtain the properties of the wall material by means of an iterative method. The remaining fluid region will be computed by a spectral hp method. We present MD simulations of a Couette flow, and the results of the developed boundary conditions from the wall fluid interaction.