The report of the Executive Committee for 2017 is presented.

This work presents a technique for extracting the detailed shape of peaks with extended, overlapping tails in an X-ray powder diffraction pattern. The application discussed here concerns crystallite size broadening, though the technique can be applied to spectra of any origin and without regard to how the profiles are to be subsequently analyzed. Historically, the extraction of profile shapes has been difficult due to the complexity of determining the background under the peak, resulting in an offset of the low-frequency components of the Fourier transform of the peak known as the `hook' problem. The use of a carefully considered statistical weighting function in a non-linear least-squares fit, followed by summing the residuals from such a fit with the fit itself, allows one to extract the full shape of an isolated peak, without contributions from either the background or adjacent peaks. The extracted shape, consisting of the fit function recombined with the residuals, is not dependent on any specific shape model. The application of this to analysis of microstructure is performed independently of global parametric models, which would reduce the number of refined parameters; therefore the technique requires high-quality data to produce results of interest. The effectiveness of the technique is demonstrated by extraction of Fourier transforms of peaks from two sets of size-broadened materials with two differing pieces of equipment.

The standard method of joint probability distribution functions, so crucial for the development of direct methods, has been revisited and updated. It consists of three steps: identification of the reflections which may contribute to the estimation of a given structure invariant or seminvariant, calculation of the corresponding joint probability distribution, and derivation of the conditional distribution of the invariant or seminvariant phase given the values of some diffracted amplitudes. In this article the conditional distributions are derived directly without passing through the second step. A good feature of direct methods is that they may work in the absence of any prior information: that is also their weakness. Different types of prior information have been taken into consideration: interatomic distances, interatomic vectors, Patterson peaks, structural model. The method of directly deriving the conditional distributions has been applied to those cases. Some new formulas have been obtained estimating two-, three- and four-phase invariants. Special attention has been dedicated to the practical aspects of the new formulas, in order to simplify their possible use in direct phasing procedures.

The development of antiviral strategies requires a clear understanding of the principles that control the protein arrangements in viral shells. Considered here are those capsids that violate the paradigmatic Caspar and Klug (CK) model, and it is shown that the important structural features of such anomalous shells from the Picobirnaviridae, Flaviviridae and Leviviridae families can be revealed by models in the form of spherical icosahedral packings of equivalent rhombic structural units (SUs). These SUs are composed of protein dimers forming the investigated capsids which, as shown here, are based on the rhombic triacontahedron (RT) geometry. How to modify the original CK approach in order to make it compatible with the considered rhombic tessellations of a sphere is also discussed. Analogies between capsids self-assembled from dimers and trimers are demonstrated. This analysis reveals the principles controlling the localization of receptor proteins (which recognize the host cell) on the capsid surface.

An elementary method is described for finding the coordination sequences for a tiling, based on coloring the underlying graph. The first application is to the two kinds of vertices (tetravalent and trivalent) in the Cairo (or dual-32.4.3.4) tiling. The coordination sequence for a tetravalent vertex turns out, surprisingly, to be 1, 4, 8, 12, 16, …, the same as for a vertex in the familiar square (or 44) tiling. The authors thought that such a simple fact should have a simple proof, and this article is the result. The method is also used to obtain coordination sequences for the 32.4.3.4, 3.4.6.4, 4.82, 3.122 and 34.6 uniform tilings, and the snub-632 tiling. In several cases the results provide proofs for previously conjectured formulas.

The unit-cell reduction described by Selling and used by Delone (whose early publications were under the spelling Delaunay) is explained in a simple form. The transformations needed to implement the reduction are listed. The simplicity of this reduction contrasts with the complexity of Niggli reduction.

The structural and nonlinear optical properties of a new anilinium hybrid crystal of chemical formula (C6H7NCl+·NO3−)3 have been investigated. The crystal structure was determined from single-crystal X-ray diffraction measurements performed at a temperature of 100 K which show that the compound crystallizes in a noncentrosymmetric space group (Pna21). The structural analysis was coupled with Hirshfeld surface analysis to evaluate the contribution of the different intermolecular interactions to the formation of supramolecular assemblies in the solid state that exhibit nonlinear optical features. This analysis reveals that the studied compound is characterized by a three-dimensional network of hydrogen bonds and the main contributions are provided by the O...H, C...H, H...H and Cl...H interactions, which alone represent ∼85% of the total contributions to the Hirshfeld surfaces. It is noteworthy that the halogen...H contributions are quite comparable with those of the H...H contacts. The nonlinear optical properties were investigated by nonlinear diffuse femtosecond-pulse reflectometry and the obtained results were compared with those of the reference material LiNbO3. The hybrid crystals exhibit notable second (SHG) and third (THG) harmonic generation which confirms its polarity is generated by the different intermolecular interactions. These measurements also highlight that the THG signal of the new anilinium compound normalized to its SHG counterpart is more pronounced than for LiNbO3.

In this article, a framework is presented that allows the systematic derivation of planar edge-to-edge k-isocoronal tilings from tile-s-transitive tilings, s ≤ k. A tiling {cal T} is k-isocoronal if its vertex coronae form k orbits or k transitivity classes under the action of its symmetry group. The vertex corona of a vertex x of {cal T} is used to refer to the tiles that are incident to x. The k-isocoronal tilings include the vertex-k-transitive tilings (k-isogonal) and k-uniform tilings. In a vertex-k-transitive tiling, the vertices form k transitivity classes under its symmetry group. If this tiling consists of regular polygons then it is k-uniform. This article also presents the classification of isocoronal tilings in the Euclidean plane.

Compared with X-rays, electron diffraction faces a crucial challenge: dynamical electron scattering compromises structure solution and its effects can only be modelled in specific cases. Dynamical scattering can be reduced experimentally by decreasing crystal size but not without a penalty, as it also reduces the overall diffracted intensity. In this article it is shown that nanometre-sized crystals from organic pharmaceuticals allow positional refinement of the hydrogen atoms, even whilst ignoring the effects of dynamical scattering during refinement. To boost the very weak diffraction data, a highly sensitive hybrid pixel detector was employed. A general likelihood-based computational approach was also introduced for further reducing the adverse effects of dynamic scattering, which significantly improved model accuracy, even for protein crystal data at substantially lower resolution.

It is shown how to reconstruct the stacking sequence from the pairwise correlation functions between layers in close-packed structures. First, of theoretical interest, the analytical formulation and solution of the problem are presented when the exact pairwise correlation counts are known. In the second part, the practical problem is approached. A simulated annealing procedure is developed to solve the problem using as initial guess approximate solutions from previous treatments. The robustness of the procedure is tested with synthetic data, followed by an experimental example. The developed approach performs robustly over different synthetic and experimental data, comparing favorably with the reported methods.

The methods for X-ray crystal orientation are rapidly evolving towards versatility, fewer goniometry measurements, automation, high accuracy and precision. One method that attracts a lot of attention is energy-dispersive X-ray diffraction (EDXRD) which is based on detecting reflections from crystallographic planes in a crystal at fixed angles of a parallel polychromatic X-ray incident beam. In theory, an EDXRD peak can move in a diffraction pattern as a function of a crystallographic plane d-spacing and its orientation relative to a fixed direction in space can change also. This is equivalent to the possibility of measuring the orientation of single crystals. The article provides a modeling for the EDXRD method whose main feature is the nonmoving crystal in the sense of traditional goniometry where the angle measurements of diffracting planes are a must. The article defines the equation of orientation for the method and shows the derivation in great detail. It is shown that the exact solutions of the equations can be obtained using the generalized reduced gradient method, a mathematical subroutine that is implemented in Excel software. The significance and scientific impact of the work are discussed along with the validated tested results.

A new aspherical scattering factor formalism has been implemented in the crystallographic least-squares refinement program SHELXL. The formalism relies on Gaussian functions and can optionally complement the independent atom model to take into account the deformation of electron-density distribution due to chemical bonding and lone pairs. Asphericity contributions were derived from the electron density obtained from quantum-chemical density functional theory computations of suitable model compounds that contain particular chemical environments, as defined by the invariom formalism. Thanks to a new algorithm, invariom assignment for refinement in SHELXL is automated. A suitable parameterization for each chemical environment within the new model was achieved by metaheuristics. Figures of merit, precision and accuracy of crystallographic least-squares refinements improve significantly upon using the new model.

One of the essential components of molecular electronic circuits are switching elements that are stable in two different states and can ideally be switched on and off many times. Here, distinct buckminsterfullerenes within a self-assembled monolayer, forming a two-dimensional dodecagonal quasicrystal on a Pt-terminated Pt3Ti(111) surface, are identified to form well separated molecular rotational switching elements. Employing scanning tunneling microscopy, the molecular-orbital appearance of the fullerenes in the quasicrystalline monolayer is resolved. Thus, fullerenes adsorbed on the 36 vertex configuration are identified to exhibit a distinctly increased mobility. In addition, this finding is verified by differential conductance measurements. The rotation of these mobile fullerenes can be triggered frequently by applied voltage pulses, while keeping the neighboring molecules immobile. An extensive analysis reveals that crystallographic and energetic constraints at the molecule/metal interface induce an inequality of the local potentials for the 36 and 32.4.3.4 vertex sites and this accounts for the switching ability of fullerenes on the 36 vertex sites. Consequently, a local area of the 8/3 approximant in the two-dimensional fullerene quasicrystal consists of single rotational switching fullerenes embedded in a matrix of inert molecules. Furthermore, it is deduced that optimization of the intermolecular interactions between neighboring fullerenes hinders the realization of translational periodicity in the fullerene monolayer on the Pt-terminated Pt3Ti(111) surface.

To date X-ray protein crystallography is the most successful technique available for the determination of high-resolution 3D structures of biological molecules and their complexes. In X-ray protein crystallography the structure of a protein is refined against the set of observed Bragg reflections from a protein crystal. The resolution of the refined protein structure is limited by the highest angle at which Bragg reflections can be observed. In addition, the Bragg reflections alone are typically insufficient (by a factor of two) to determine the structure ab initio, and so prior information is required. Crystals formed from an imperfect packing of the protein molecules may also exhibit continuous diffraction between and beyond these Bragg reflections. When this is due to random displacements of the molecules from each crystal lattice site, the continuous diffraction provides the necessary information to determine the protein structure without prior knowledge, to a resolution that is not limited by the angular extent of the observed Bragg reflections but instead by that of the diffraction as a whole. This article presents an iterative projection algorithm that simultaneously uses the continuous diffraction as well as the Bragg reflections for the determination of protein structures. The viability of this method is demonstrated on simulated crystal diffraction.

Diffuse scattering is a rich source of information about disorder in crystalline materials, which can be modelled using atomistic techniques such as Monte Carlo and molecular dynamics simulations. Modern X-ray and neutron scattering instruments can rapidly measure large volumes of diffuse-scattering data. Unfortunately, current algorithms for atomistic diffuse-scattering calculations are too slow to model large data sets completely, because the fast Fourier transform (FFT) algorithm has long been considered unsuitable for such calculations [Butler & Welberry (1992). J. Appl. Cryst. 25, 391–399]. Here, a new approach is presented for ultrafast calculation of atomistic diffuse-scattering patterns. It is shown that the FFT can actually be used to perform such calculations rapidly, and that a fast method based on sampling theory can be used to reduce high-frequency noise in the calculations. These algorithms are benchmarked using realistic examples of compositional, magnetic and displacive disorder. They accelerate the calculations by a factor of at least 102, making refinement of atomistic models to large diffuse-scattering volumes practical.

This work considers the scaling properties characterizing the hyperuniformity (or anti-hyperuniformity) of long-wavelength fluctuations in a broad class of one-dimensional substitution tilings. A simple argument is presented which predicts the exponent α governing the scaling of Fourier intensities at small wavenumbers, tilings with α > 0 being hyperuniform, and numerical computations confirm that the predictions are accurate for quasiperiodic tilings, tilings with singular continuous spectra and limit-periodic tilings. Quasiperiodic or singular continuous cases can be constructed with α arbitrarily close to any given value between −1 and 3. Limit-periodic tilings can be constructed with α between −1 and 1 or with Fourier intensities that approach zero faster than any power law.

In a Model block, if the instance-specific value of a model argument has an anonymous structure type, an update diagram reports an error when there are multiple bus types that match that anonymous structure type.This bug exists in the following release(s):
R2020a

This bug has a workaround

Interested in Upgrading?

Attempting to display the contents of an HDF5 file containing fixed-length UTF-8 encoded strings results in an unexpected error in MATLAB.

For example, the following code

  h5disp('myHDF5FileWithFixedLenUTF8Strings.h5')

returns this error:

Error using h5infoc
UTF-8 encoding is only supported for variable length strings.

Error in h5info (line 108)
hinfo = h5infoc(filename,location, useUtf8);

Error in h5disp>display_hdf5 (line 121)
hinfo = h5info(options.Filename,options.Location);

Error in h5disp (line 99)
display_hdf5(options);

When creating a System Composer architecture model, if the following functions are called with an invalid or missing stereotype, it can cause the model to get corrupted and might lead to a MATLAB crash:

• addComponent
• addPort
• connect

Simulink might produce an incorrect sorted order for a model that meets all of the following conditions:

• The model contains blocks from the SoC Blockset
• The Signal logging option is selected in the model configuration set
• Signals using asynchronous sample time are configured for logging

If you try to open a text file with UTF-16 or UTF-32 encoding, the Import Tool displays a warning message stating that the encoding is not supported.  If you continue to load the file anyways, it is opened with UTF-8 encoding, and the file may not be displayed or imported as expected.This bug exists in the following release(s):
R2020a

This bug has a workaround

Interested in Upgrading?

When the Signal Editor is saving data, the indicator that the save is occurring does not appear.  You might notice a delay when saving large data files.This bug exists in the following release(s):
R2020a

Interested in Upgrading?

In certain scenarios, a Subsystem Reference block changes its position. It can happen during: 

1. Loading of top model which contains the Subsystem Reference block 
2. Synchronization of Subsystem Reference block's contents
3. Converting a subsystem to Subsystem Reference block