The intrinsic twinning and enigmatic twisting of aragonite crystals.

It is generally accepted that aragonite crystals of biogenic origin are characterized by significantly higher twin densities compared to samples formed during geological processes. Based on our single crystal X-ray diffraction (SCXRD) and transmission electron microscopy (TEM) study of aragonite crystals from various localities, we show that in geological aragonites, the twin densities are comparable to those of the samples from crossed lamellar zones of molluscs shells. The high twin density is consistent with performed calculations, according to which the Gibbs free energy of twin-free aragonite is close to that of periodically twinned aragonite structure. In some cases, high twin densities result in the appearance of diffuse scattering in SCXRD patterns. The obtained TEM and optical micrographs show that besides the twin boundaries (TBs) of growth origin, there are also TBs and especially stacking faults that were likely formed as the result of local strain compensation. SCXRD patterns of the samples from Tazouta, in addition to diffuse scattering lines, show Debye arcs in the [Formula: see text] plane. These Debye arcs are present only on one side of the Bragg reflections and have an azimuthal extent of nearly 30°, making the whole symmetry of the diffraction pattern distinctly chiral, which has not yet been reported for aragonite. By analogy with biogenic calcite crystals, we associate these arcs with the presence of misoriented subgrains formed as a result of crystal twisting during growth.

Controlling inversion disorder in a stoichiometric spinel magnet.

In the study of frustrated quantum magnets, it is essential to be able to control the nature and degree of site disorder during the growth process, as many measurement techniques are incapable of distinguishing between site disorder and frustration-induced spin disorder. Pyrochlore-structured spinel oxides can serve as model systems of geometrically frustrated three-dimensional quantum magnets; however, the nature of the magnetism in one well-studied spinel, ZnFe2O4, remains unclear. Here, we demonstrate simultaneous control of both stoichiometry and inversion disorder in the growth of ZnFe2O4 single crystals, directly yielding a revised understanding of both the collective spin behavior and lattice symmetry. Crystals grown in the stoichiometric limit with minimal site inversion disorder contravene all the previously suggested exotic spin phases in ZnFe2O4. Furthermore, the structure is confirmed on the [Formula: see text] space group with broken inversion symmetry that induces antiferroelectricity. The effective tuning of magnetic behavior by site disorder in the presence of robust antiferroelectricity makes ZnFe2O4 of special interest to multiferroic devices.

Crystal Chemistry of the Copper Oxalate Biomineral Moolooite: The First Single-Crystal X-ray Diffraction Studies and Thermal Behavior.

Moolooite, Cu(C2O4)·nH2O, is a typical biomineral which forms due to Cu-bearing minerals coming into contact with oxalic acid sources such as bird guano deposits or lichens, and no single crystals of moolooite of either natural or synthetic origin have been found yet. This paper reports, for the first time, on the preparation of single crystals of a synthetic analog of the copper-oxalate biomineral moolooite, and on the refinement of its crystal structure from the single-crystal X-ray diffraction (SCXRD) data. Along with the structural model, the SCXRD experiment showed the significant contribution of diffuse scattering to the overall diffraction data, which comes from the nanostructural disorder caused by stacking faults of Cu oxalate chains as they lengthen. This type of disorder should result in the chains breaking, at which point the H2O molecules may be arranged. The amount of water in the studied samples did not exceed 0.15 H2O molecules per formula unit. Apparently, the mechanism of incorporation of H2O molecules governs the absence of good-quality single crystals in nature and a lack of them in synthetic experiments: the more H2O content in the structure, the stronger the disorder will be. A description of the crystal structure indicates that the ideal structure of the Cu oxalate biomineral moolooite should not contain H2O molecules and should be described by the Cu(C2O4) formula. However, it was shown that natural and synthetic moolooite crystals contain a significant portion of "structural" water, which cannot be ignored. Considering the substantially variable amount of water, which can be incorporated into the crystal structure, the formula Cu(C2O4)·nH2O for moolooite is justified.

Optical Differentiation of Brain Tumors Based on Raman Spectroscopy and Cluster Analysis Methods.

In the present study, various combinations of dimensionality reduction methods with data clustering methods for the analysis of biopsy samples of intracranial tumors were investigated. Fresh biopsies of intracranial tumors were studied in the Laboratory of Neurosurgical Anatomy and Preservation of Biological Materials of N.N. Burdenko Neurosurgery Medical Center no later than 4 h after surgery. The spectra of Protoporphyrin IX (Pp IX) fluorescence, diffuse reflectance (DR) and Raman scattering (RS) of biopsy samples were recorded. Diffuse reflectance studies were carried out using a white light source in the visible region. Raman scattering spectra were obtained using a 785 nm laser. Patients diagnosed with meningioma, glioblastoma, oligodendroglioma, and astrocytoma were studied. We used the cluster analysis method to detect natural clusters in the data sample presented in the feature space formed based on the spectrum analysis. For data analysis, four clustering algorithms with eight dimensionality reduction algorithms were considered.

Towards performing high-resolution inelastic X-ray scattering measurements at hard X-ray free-electron lasers coupled with energetic laser drivers.

High-resolution inelastic X-ray scattering is an established technique in the synchrotron community, used to investigate collective low-frequency responses of materials. When fielded at hard X-ray free-electron lasers (XFELs) and combined with high-intensity laser drivers, it becomes a promising technique for investigating matter at high temperatures and high pressures. This technique gives access to important thermodynamic properties of matter at extreme conditions, such as temperature, material sound speed, and viscosity. The successful realization of this method requires the acquisition of many identical laser-pump/X-ray-probe shots, allowing the collection of a sufficient number of photons necessary to perform quantitative analyses. Here, a 2.5-fold improvement in the energy resolution of the instrument relative to previous works at the Matter in Extreme Conditions (MEC) endstation, Linac Coherent Light Source (LCLS), and the High Energy Density (HED) instrument, European XFEL, is presented. Some aspects of the experimental design that are essential for improving the number of photons detected in each X-ray shot, making such measurements feasible, are discussed. A careful choice of the energy resolution, the X-ray beam mode provided by the XFEL, and the position of the analysers used in such experiments can provide a more than ten-fold improvement in the photometrics. The discussion is supported by experimental data on 10 µm-thick iron and 50 nm-thick gold samples collected at the MEC endstation at the LCLS, and by complementary ray-tracing simulations coupled with thermal diffuse scattering calculations.

The ePix10k 2-megapixel hard X-ray detector at LCLS.

The ePix10ka2M (ePix10k) is a new large area detector specifically developed for X-ray free-electron laser (XFEL) applications. The hybrid pixel detector was developed at SLAC to provide a hard X-ray area detector with a high dynamic range, running at the 120 Hz repetition rate of the Linac Coherent Light Source (LCLS). The ePix10k consists of 16 modules, each with 352 × 384 pixels of 100 µm × 100 µm distributed on four ASICs, resulting in a 2.16 megapixel detector, with a 16.5 cm × 16.5 cm active area and ∼80% coverage. The high dynamic range is achieved with three distinct gain settings (low, medium, high) as well as two auto-ranging modes (high-to-low and medium-to-low). Here the three fixed gain modes are evaluated. The resulting dynamic range (from single photon counting to 10000 photons pixel-1 pulse-1 at 8 keV) makes it suitable for a large number of different XFEL experiments. The ePix10k replaces the large CSPAD in operation since 2011. The dimensions of the two detectors are similar, making the upgrade from CSPAD to ePix10k straightforward for most setups, with the ePix10k improving on experimental performance. The SLAC-developed ePix cameras all utilize a similar platform, are tailored to target different experimental conditions and are designed to provide an upgrade path for future high-repetition-rate XFELs. Here the first measurements on this new ePix10k detector are presented and the performance under typical XFEL conditions evaluated during an LCLS X-ray diffuse scattering experiment measuring the 9.5 keV X-ray photons scattered from a thin liquid jet.

Expanding the Dimensions of a Small, Two-Dimensional Diffraction Detector.

We report an approach to expand the effective number of pixels available to small, two-dimensional electron detectors. To do so, we acquire subsections of a diffraction pattern that are then accurately stitched together in post-processing. Using an electron microscopy pixel array detector (EMPAD) that has only 128 × 128 pixels, we show that the field of view can be expanded while achieving high reciprocal-space sampling. Further, we highlight the need to properly account for the detector position (rotation) and the non-orthonormal diffraction shift axes to achieve an accurate reconstruction. Applying the method, we provide examples of spot and convergent beam diffraction patterns acquired with a pixelated detector.

A new diffractometer for diffuse scattering studies on the ID28 beamline at the ESRF.

A new diffractometer is now available to the general user community at the ESRF. The new diffractometer is a side station of the high-resolution inelastic X-ray scattering spectrometer on beamline ID28 and is located in the same experimental hutch. Both instruments can be operated simultaneously. The new diffractometer combines a fast and low-noise hybrid pixel detector with a variable diffraction geometry. The beam spot on the sample is 50 µm × 50 µm, where focusing is achieved by a combination of Be lenses and a KB mirror. Wavelengths from 0.5 to 0.8 Šcan be used for the diffraction experiments. The setup is compatible with a variety of sample environments, allowing studies under non-ambient conditions. The diffractometer is optimized to allow a rapid survey of reciprocal space and diffuse scattering for the identification of regions of interest for subsequent inelastic scattering studies, but can also be employed as a fully independent station for structural studies from both powder and single-crystal diffraction experiments. Several software packages for the transformation and visualization of diffraction data are available. An analysis of data collected with the new diffractometer shows that the ID28 side station is a state-of-the-art instrument for structural investigations using diffraction and diffuse scattering experiments.

Measuring and modeling diffuse scattering in protein X-ray crystallography.

X-ray diffraction has the potential to provide rich information about the structural dynamics of macromolecules. To realize this potential, both Bragg scattering, which is currently used to derive macromolecular structures, and diffuse scattering, which reports on correlations in charge density variations, must be measured. Until now, measurement of diffuse scattering from protein crystals has been scarce because of the extra effort of collecting diffuse data. Here, we present 3D measurements of diffuse intensity collected from crystals of the enzymes cyclophilin A and trypsin. The measurements were obtained from the same X-ray diffraction images as the Bragg data, using best practices for standard data collection. To model the underlying dynamics in a practical way that could be used during structure refinement, we tested translation-libration-screw (TLS), liquid-like motions (LLM), and coarse-grained normal-modes (NM) models of protein motions. The LLM model provides a global picture of motions and was refined against the diffuse data, whereas the TLS and NM models provide more detailed and distinct descriptions of atom displacements, and only used information from the Bragg data. Whereas different TLS groupings yielded similar Bragg intensities, they yielded different diffuse intensities, none of which agreed well with the data. In contrast, both the LLM and NM models agreed substantially with the diffuse data. These results demonstrate a realistic path to increase the number of diffuse datasets available to the wider biosciences community and indicate that dynamics-inspired NM structural models can simultaneously agree with both Bragg and diffuse scattering.

Origins of ultralow velocity zones through slab-derived metallic melt.

Understanding the ultralow velocity zones (ULVZs) places constraints on the chemical composition and thermal structure of deep Earth and provides critical information on the dynamics of large-scale mantle convection, but their origin has remained enigmatic for decades. Recent studies suggest that metallic iron and carbon are produced in subducted slabs when they sink beyond a depth of 250 km. Here we show that the eutectic melting curve of the iron-carbon system crosses the current geotherm near Earth's core-mantle boundary, suggesting that dense metallic melt may form in the lowermost mantle. If concentrated into isolated patches, such melt could produce the seismically observed density and velocity features of ULVZs. Depending on the wetting behavior of the metallic melt, the resultant ULVZs may be short-lived domains that are replenished or regenerated through subduction, or long-lasting regions containing both metallic and silicate melts. Slab-derived metallic melt may produce another type of ULVZ that escapes core sequestration by reacting with the mantle to form iron-rich postbridgmanite or ferropericlase. The hypotheses connect peculiar features near Earth's core-mantle boundary to subduction of the oceanic lithosphere through the deep carbon cycle.

X-ray scattering by the fused silica surface etched by low-energy Ar ions.

Anomalously high x-ray scattering at a wavelength of 0.154 nm by super-polished substrates of fused silica, which were etched by the argon ions with the energy of 300 eV, is detected. The scattering intensity increases monotonically with increasing of the etching depth. The effect is explained by the scattering on the volume inhomogeneities with the lateral size greater than 0.5 µm of the subsurface "damaged" layer. The concentration of volume inhomogeneities increases with the increase of the fluence of argon ions, but the concentration of implanted argon atoms in the layer quickly reaches the maximum value and then begins a trend of going down. The thickness of the "damaged" layer is approximately equal to the penetration depth of the Ar atoms and can be directly determined from the x-ray specular reflection. It is shown that the presence of volume inhomogeneities of the subsurface "damaged" layer does not affect the geometric roughness of the surface. The observed effect imposes limitations on the usage of grazing incidence x-ray optics without reflective coatings and of the diffuse x-ray scattering (DXRS) method for studying the substrate roughness. A new method that potentially enables to evaluate the applicability of the DXRS method in practice is proposed.

HIV-1 matrix-31 membrane binding peptide interacts differently with membranes containing PS vs. PI(4,5)P2.

Efficient assembly of HIV-1 at the plasma membrane (PM) of the T-cell specifically requires PI(4,5)P2. It was previously shown that a highly basic region (HBR) of the matrix protein (MA) on the Gag precursor polyprotein Pr55Gag is required for membrane association. MA is N-terminally myristoylated, which enhances its affinity to membranes. In this work we used X-ray scattering and neutron reflectivity to determine how the physical properties and structure of lipid bilayers respond to the addition of binding domain peptides, either in the myristoylated form (MA31myr) or without the myristoyl group (MA31). Neutron reflectivity measurements showed the peptides predominantly located in the hydrocarbon interior. Diffuse X-ray scattering showed differences in membrane properties upon addition of peptides and the direction of the changes depended on lipid composition. The PI(4,5)P2-containing bilayers softened, thinned and became less ordered as peptide concentration increased. In contrast, POPS-containing bilayers with equivalent net charge first stiffened, thickened and became more ordered with increasing peptide concentration. As softening the host cell's PM upon contact with the protein lowers the free energy for membrane restructuring, thereby potentially facilitating budding of viral particles, our results suggest that the role of PI(4,5)P2 in viral assembly goes beyond specific stereochemical membrane binding. These studies reinforce the importance of lipids in virology.

Conformational dynamics of a crystalline protein from microsecond-scale molecular dynamics simulations and diffuse X-ray scattering.

X-ray diffraction from protein crystals includes both sharply peaked Bragg reflections and diffuse intensity between the peaks. The information in Bragg scattering is limited to what is available in the mean electron density. The diffuse scattering arises from correlations in the electron density variations and therefore contains information about collective motions in proteins. Previous studies using molecular-dynamics (MD) simulations to model diffuse scattering have been hindered by insufficient sampling of the conformational ensemble. To overcome this issue, we have performed a 1.1-µs MD simulation of crystalline staphylococcal nuclease, providing 100-fold more sampling than previous studies. This simulation enables reproducible calculations of the diffuse intensity and predicts functionally important motions, including transitions among at least eight metastable states with different active-site geometries. The total diffuse intensity calculated using the MD model is highly correlated with the experimental data. In particular, there is excellent agreement for the isotropic component of the diffuse intensity, and substantial but weaker agreement for the anisotropic component. Decomposition of the MD model into protein and solvent components indicates that protein-solvent interactions contribute substantially to the overall diffuse intensity. We conclude that diffuse scattering can be used to validate predictions from MD simulations and can provide information to improve MD models of protein motions.

Crystal pathologies in macromolecular crystallography.

Macromolecules, such as proteins or nucleic acids, form crystals with a large volume fraction of water, ~50% on average. Apart from typical physical defects and rather trivial poor quality problems, macromolecular crystals, as essentially any crystals, can also suffer from several kinds of pathologies, in which everything seems to be perfect, except that from the structural point of view the interpretation may be very difficult, sometimes even impossible. A frequent nuisance is pseudosymmetry, or non-crystallographic symmetry (NCS), which is particularly nasty when it has translational character. Lattice-translocation defects, also called order-disorder twinning (OD-twinning), occur when molecules are packed regularly in layers but the layers are stacked (without rotation) in two (or more) discrete modes, with a unique translocation vector. Crystal twinning arises when twin domains have different orientations, incompatible with the symmetry of the crystal structure. There are also crystals in which the periodic (lattice) order is broken or absent altogether. When the strict short-range translational order from one unit cell to the next is lost but the long-range order is restored by a periodic modulation, we have a modulated crystal structure. In quasicrystals (not observed for macromolecules yet), the periodic order (in 3D space) is lost completely and the diffraction pattern (which is still discrete) cannot be even indexed using three hkl indices. In addition, there are other physical defects and phenomena (such as high mosaicity, diffraction anisotropy, diffuse scattering, etc.) which make diffraction data processing and structure solution difficult or even impossible.

Predicting X-ray diffuse scattering from translation-libration-screw structural ensembles.

Identifying the intramolecular motions of proteins and nucleic acids is a major challenge in macromolecular X-ray crystallography. Because Bragg diffraction describes the average positional distribution of crystalline atoms with imperfect precision, the resulting electron density can be compatible with multiple models of motion. Diffuse X-ray scattering can reduce this degeneracy by reporting on correlated atomic displacements. Although recent technological advances are increasing the potential to accurately measure diffuse scattering, computational modeling and validation tools are still needed to quantify the agreement between experimental data and different parameterizations of crystalline disorder. A new tool, phenix.diffuse, addresses this need by employing Guinier's equation to calculate diffuse scattering from Protein Data Bank (PDB)-formatted structural ensembles. As an example case, phenix.diffuse is applied to translation-libration-screw (TLS) refinement, which models rigid-body displacement for segments of the macromolecule. To enable the calculation of diffuse scattering from TLS-refined structures, phenix.tls_as_xyz builds multi-model PDB files that sample the underlying T, L and S tensors. In the glycerophosphodiesterase GpdQ, alternative TLS-group partitioning and different motional correlations between groups yield markedly dissimilar diffuse scattering maps with distinct implications for molecular mechanism and allostery. These methods demonstrate how, in principle, X-ray diffuse scattering could extend macromolecular structural refinement, validation and analysis.

From deep TLS validation to ensembles of atomic models built from elemental motions.

The translation-libration-screw model first introduced by Cruickshank, Schomaker and Trueblood describes the concerted motions of atomic groups. Using TLS models can improve the agreement between calculated and experimental diffraction data. Because the T, L and S matrices describe a combination of atomic vibrations and librations, TLS models can also potentially shed light on molecular mechanisms involving correlated motions. However, this use of TLS models in mechanistic studies is hampered by the difficulties in translating the results of refinement into molecular movement or a structural ensemble. To convert the matrices into a constituent molecular movement, the matrix elements must satisfy several conditions. Refining the T, L and S matrix elements as independent parameters without taking these conditions into account may result in matrices that do not represent concerted molecular movements. Here, a mathematical framework and the computational tools to analyze TLS matrices, resulting in either explicit decomposition into descriptions of the underlying motions or a report of broken conditions, are described. The description of valid underlying motions can then be output as a structural ensemble. All methods are implemented as part of the PHENIX project.

Correction of complex nonlinear signal response from a pixel array detector.

The pulsed free-electron laser light sources represent a new challenge to photon area detectors due to the intrinsic spontaneous X-ray photon generation process that makes single-pulse detection necessary. Intensity fluctuations up to 100% between individual pulses lead to high linearity requirements in order to distinguish small signal changes. In real detectors, signal distortions as a function of the intensity distribution on the entire detector can occur. Here a robust method to correct this nonlinear response in an area detector is presented for the case of exposures to similar signals. The method is tested for the case of diffuse scattering from liquids where relevant sub-1% signal changes appear on the same order as artifacts induced by the detector electronics.

Modelling dynamical 3D electron diffraction intensities. II. The role of inelastic scattering.

The strong interaction of high-energy electrons with a crystal results in both dynamical elastic scattering and inelastic events, particularly phonon and plasmon excitation, which have relatively large cross sections. For accurate crystal structure refinement it is therefore important to uncover the impact of inelastic scattering on the Bragg beam intensities. Here a combined Bloch wave-Monte Carlo method is used to simulate phonon and plasmon scattering in crystals. The simulated thermal and plasmon diffuse scattering are consistent with experimental results. The simulations also confirm the empirical observation of a weaker unscattered beam intensity with increasing energy loss in the low-loss regime, while the Bragg-diffracted beam intensities do not change significantly. The beam intensities include the diffuse scattered background and have been normalized to adjust for the inelastic scattering cross section. It is speculated that the random azimuthal scattering angle during inelastic events transfers part of the unscattered beam intensity to the inner Bragg reflections. Inelastic scattering should not significantly influence crystal structure refinement, provided there are no artefacts from any background subtraction, since the relative intensity of the diffracted beams (which includes the diffuse scattering) remains approximately constant in the low energy loss regime.

Scaling and merging macromolecular diffuse scattering with mdx2.

Diffuse scattering is a promising method to gain additional insight into protein dynamics from macromolecular crystallography (MX) experiments. Bragg intensities yield the average electron density, while the diffuse scattering can be processed to obtain a three-dimensional reciprocal space map, that is further analyzed to determine correlated motion. To make diffuse scattering techniques more accessible, we have created software for data processing called mdx2 that is both convenient to use and simple to extend and modify. Mdx2 is written in Python, and it interfaces with DIALS to implement self-contained data reduction workflows. Data are stored in NeXus format for software interchange and convenient visualization. Mdx2 can be run on the command line or imported as a package, for instance to encapsulate a complete workflow in a Jupyter notebook for reproducible computing and education. Here, we describe mdx2 version 1.0, a new release incorporating state-of-the-art techniques for data reduction. We describe the implementation of a complete multi-crystal scaling and merging workflow, and test the methods using a high-redundancy dataset from cubic insulin. We show that redundancy can be leveraged during scaling to correct systematic errors, and obtain accurate and reproducible measurements of weak diffuse signals.

Scaling and merging macromolecular diffuse scattering with mdx2.

Diffuse scattering is a promising method to gain additional insight into protein dynamics from macromolecular crystallography experiments. Bragg intensities yield the average electron density, while the diffuse scattering can be processed to obtain a three-dimensional reciprocal-space map that is further analyzed to determine correlated motion. To make diffuse scattering techniques more accessible, software for data processing called mdx2 has been created that is both convenient to use and simple to extend and modify. mdx2 is written in Python, and it interfaces with DIALS to implement self-contained data-reduction workflows. Data are stored in NeXus format for software interchange and convenient visualization. mdx2 can be run on the command line or imported as a package, for instance to encapsulate a complete workflow in a Jupyter notebook for reproducible computing and education. Here, mdx2 version 1.0 is described, a new release incorporating state-of-the-art techniques for data reduction. The implementation of a complete multi-crystal scaling and merging workflow is described, and the methods are tested using a high-redundancy data set from cubic insulin. It is shown that redundancy can be leveraged during scaling to correct systematic errors and obtain accurate and reproducible measurements of weak diffuse signals.