Acta Crystallographica Section A: Foundations and Advances最新文献

Analytical calculations of absorption corrections for X-ray powder diffraction experiments on non-ideal samples with surface roughness, porosity or absorption contrasts from multiple phases require complex mathematical models to represent their material distribution. In a computational approach to this problem, a practicable ray-tracing algorithm is formulated which is capable of simulating angle-dependent absorption corrections in reflection geometry for any given rasterized sample model. Single or multiphase systems with arbitrary surface roughness, porosity and spatial distribution of the phases in any combination can be modeled on a voxel grid by assigning respective values to each voxel. The absorption corrections are calculated by tracing the attenuation of X-rays along their individual paths via a modified shear-warp algorithm. The algorithm is presented in detail and the results of simulated absorption corrections on samples with various surface modulations are discussed in the context of published experimental results.

Fundamental equations describing the X-ray and electron diffraction scattering in imperfect crystals have been derived in the form of the matrix Fredholm-Volterra integral equation of the second kind. A theoretical approach has been developed using the perfect-crystal Green function formalism. In contrast, another approach utilizes the wavefield eigenfunctions related to the diagonalized matrix propagators of the conventional Takagi-Taupin and Howie-Whelan equations. Using the Liouville-Neumann-type series formalism for building up the matrix Fredholm-Volterra integral equation solutions, the general resolvent function solutions of the X-ray and electron diffraction boundary-valued Cauchy problems have been obtained. Based on the resolvent-type solutions, the aim is to reveal the features of the diffraction scattering onto the crystal lattice defects, including the mechanisms of intra- and interbranch wave scattering in the strongly deformed regions in the vicinity of crystal lattice defect cores. Using the two-stage resolvent solution of the second order, this approach has been supported by straightforward calculation of the electron bright- and dark-field contrasts of an edge dislocation in a thick foil. The results obtained for the bright- and dark-field profiles of the edge dislocation are discussed and compared with analogous ones numerically calculated by Howie & Whelan [Proc. R. Soc. A (1962), 267, 206].

This paper proposes a new order parameter model which satisfactorily explains complicated symmetry changes, the temperature-pressure (T-P) phase diagram and elastic anomalies observed experimentally with the improper ferroelastic phase transitions in multiferroic KMnF₃ single crystal. First, it is shown that the order parameter model is transformed according to the four-dimensional reducible representation of the wavevector star channel group. Second, based on the order parameter model and the singularity theory, the sixth-order structurally stable Landau potential model is constructed. Finally, the theoretical T-P phase diagram is plotted and the elastic anomalies possible for each of the phase transitions are discussed.

In ab initio indexing, for a given diffraction/scattering pattern, the unit-cell parameters and the Miller indices assigned to reflections in the pattern are determined simultaneously. `Ab initio' means a process performed without any good prior information on the crystal lattice. Newly developed ab initio indexing software is frequently reported in crystallography. However, it is not widely recognized that use of a Bravais lattice determination method, which is tolerant of experimental errors, can simplify indexing algorithms and increase their success rates. One of the goals of this article is to collect information on the lattice-basis reduction theory and its applications. The main result is a Bravais lattice determination algorithm for 2D lattices, along with a mathematical proof that it works even for parameters containing large observational errors. It uses two lattice-basis reduction methods that seem to be optimal for different symmetries, similarly to the algorithm for 3D lattices implemented in the CONOGRAPH software. In indexing, a method for error-stable unit-cell identification is also required to exclude duplicate solutions. Several methods are introduced to measure the difference in unit cells known in crystallography and mathematics.

The shape of a flat-faceted octahedral crystal can be uniquely defined by the measured distances between pairs of its parallel facets and the length of one of its false edges. In total, only five numerical values are involved in this approach. Some interdependencies of parameters that allow one to control the correctness of measurements were derived. The proposed method is suitable for describing the shape as full-faceted, or as incomplete octahedral crystals (e.g. diamond) with unequally developed facets. This so-called `real crystal form' can be considered as one of the typomorphic features of minerals, connecting the dissymmetry to the anisotropy of the host rock. The measurement results can be used in crystallo-morphological analysis, restoration of the lost crystal shape in the case of man-made damage and in the practice of diamond prospecting.

Following the work of Day & Hawthorne [Acta Cryst. (2022), A78, 212-233] and Day et al. [Acta Cryst. (2024), A80, 258-281], the program GraphT-T has been developed to embed graphical representations of observed and hypothetical chains of (SiO₄)^4- tetrahedra into 2D and 3D Euclidean space. During embedding, the distance between linked vertices (T-T distances) and the distance between unlinked vertices (T...T separations) in the resultant unit-distance graph are restrained to the average observed distance between linked Si tetrahedra (3.06±0.15 Å) and the minimum separation between unlinked vertices is restrained to be equal to or greater than the minimum distance between unlinked Si tetrahedra (3.713 Å) in silicate minerals. The notional interactions between vertices are described by a 3D spring-force algorithm in which the attractive forces between linked vertices behave according to Hooke's law and the repulsive forces between unlinked vertices behave according to Coulomb's law. Embedding parameters (i.e. spring coefficient, k, and Coulomb's constant, K) are iteratively refined during embedding to determine if it is possible to embed a given graph to produce a unit-distance graph with T-T distances and T...T separations that are compatible with the observed T-T distances and T...T separations in crystal structures. The resultant unit-distance graphs are denoted as compatible and may form crystal structures if and only if all distances between linked vertices (T-T distances) agree with the average observed distance between linked Si tetrahedra (3.06±0.15 Å) and the minimum separation between unlinked vertices is equal to or greater than the minimum distance between unlinked Si tetrahedra (3.713 Å) in silicate minerals. If the unit-distance graph does not satisfy these conditions, it is considered incompatible and the corresponding chain of tetrahedra is unlikely to form crystal structures. Using GraphT-T, Day et al. [Acta Cryst. (2024), A80, 258-281] have shown that several topological properties of chain graphs influence the flexibility (and rigidity) of the corresponding chains of Si tetrahedra and may explain why particular compatible chain arrangements (and the minerals in which they occur) are more common than others and/or why incompatible chain arrangements do not occur in crystals despite being topologically possible.

A crystal structure with N atoms in its unit cell can be solved starting from a model with atoms 1 to j - 1 being located. To locate the next atom j, the method uses a modified definition of the traditional R1 factor where its dependencies on the locations of atoms j + 1 to N are removed. This modified R1 is called the single-atom R1 (sR1), because the locations of atoms 1 to j - 1 in sR1 are the known parameters, and only the location of atom j is unknown. Finding the correct position of atom j translates thus into the optimization of the sR1 function, with respect to its fractional coordinates, x_j, y_j, z_j. Using experimental data, it has been verified that an sR1 has a hole near each missing atom. Further, it has been verified that an algorithm based on sR1, hereby called the sR1 method, can solve crystal structures (with up to 156 non-hydrogen atoms in the unit cell). The strategy to carry out this calculation has also been optimized. The main feature of the sR1 method is that, starting from a single arbitrarily positioned atom, the structure is gradually revealed. With the user's help to delete poorly determined parts of the structure, the sR1 method can build the model to a high final quality. Thus, sR1 is a viable and useful tool for solving crystal structures.

In Part I of this series, all topologically possible 1-periodic infinite graphs (chain graphs) representing chains of tetrahedra with up to 6-8 vertices (tetrahedra) per repeat unit were generated. This paper examines possible restraints on embedding these chain graphs into Euclidean space such that they are compatible with the metrics of chains of tetrahedra in observed crystal structures. Chain-silicate minerals with T = Si⁴⁺ (plus P⁵⁺, V⁵⁺, As⁵⁺, Al³⁺, Fe³⁺, B³⁺, Be²⁺, Zn²⁺ and Mg²⁺) have a grand nearest-neighbour ⟨T-T⟩ distance of 3.06±0.15 Å and a minimum T...T separation of 3.71 Å between non-nearest-neighbour tetrahedra, and in order for embedded chain graphs (called unit-distance graphs) to be possible atomic arrangements in crystals, they must conform to these metrics, a process termed equalization. It is shown that equalization of all acyclic chain graphs is possible in 2D and 3D, and that equalization of most cyclic chain graphs is possible in 3D but not necessarily in 2D. All unique ways in which non-isomorphic vertices may be moved are designated modes of geometric modification. If a mode (m) is applied to an equalized unit-distance graph such that a new geometrically distinct unit-distance graph is produced without changing the lengths of any edges, the mode is designated as valid (m_v); if a new geometrically distinct unit-distance graph cannot be produced, the mode is invalid (m_i). The parameters m_v and m_i are used to define ranges of rigidity of the unit-distance graphs, and are related to the edge-to-vertex ratio, e/n, of the parent chain graph. The program GraphT-T was developed to embed any chain graph into Euclidean space subject to the metric restraints on T-T and T...T. Embedding a selection of chain graphs with differing e/n ratios shows that the principal reason why many topologically possible chains cannot occur in crystal structures is due to violation of the requirement that T...T > 3.71 Å. Such a restraint becomes increasingly restrictive as e/n increases and indicates why chains with stoichiometry TO_<2.5 do not occur in crystal structures.

The kisrhombille tiling is the dual tessellation of one of the semi-regular tessellations. It consists of right-angled triangle tiles with 12 different orientations. An adequate coordinate system for the tiles of the grid has been defined that allows a formal description of the grid. In this paper, two tiles are considered to be neighbors if they share at least one point in their boundary. Paths are sequences of tiles such that any two consecutive tiles are neighbors. The digital distance is defined as the minimum number of steps in a path between the tiles, and the distance formula is proven through constructing minimum paths. In fact, the distance between triangles is almost twice the hexagonal distance of their embedding hexagons.

Recent advances in quantum crystallography have shown that, beyond conventional charge density refinement, a one-electron reduced density matrix (1-RDM) satisfying N-representability conditions can be reconstructed using jointly experimental X-ray structure factors and directional Compton profiles (DCP) through semidefinite programming. So far, such reconstruction methods for 1-RDM, not constrained to idempotency, have been tested only on a toy model system (CO₂). In this work, a new method is assessed on crystalline urea [CO(NH₂)₂] using static (0 K) and dynamic (50 K) artificial experimental data. An improved model, including symmetry constraints and frozen core-electron contribution, is introduced to better handle the increasing system complexity. Reconstructed 1-RDMs, deformation densities and DCP anisotropy are analysed, and it is demonstrated that the changes in the model significantly improve the reconstruction quality, even when there is insufficient information and data corruption. The robustness of the model and the strategy are thus shown to be well adapted to address the reconstruction problem from actual experimental scattering data.