Simplifying the representation of complex free-energy landscapes using sketch-map

A new scheme, sketch-map, for obtaining a low-dimensional representation of the region of phase space explored during an enhanced dynamics simulation is proposed. We show evidence, from an examination of the distribution of pairwise distances between frames, that some features of the free-energy surface are inherently high-dimensional. This makes dimensionality reduction problematic because the data does not satisfy the assumptions made in conventional manifold learning algorithms We therefore propose that when dimensionality reduction is performed on trajectory data one should think of the resultant embedding as a quickly sketched set of directions rather than a road map. In other words, the embedding tells one about the connectivity between states but does not provide the vectors that correspond to the slow degrees of freedom. This realization informs the development of sketch-map, which endeavors to reproduce the proximity information from the high-dimensionality description in a space of lower dimensionality even when a faithful embedding is not possible.

General, PDB-based collective variables for protein folding

New, automated forms of data-analysis are required in order to understand the high-dimensional trajectories that are obtained from molecular dynamics simulations on proteins. Dimensionality reduction algorithms are particularly appealing in this regard as they allow one to construct unbiased, low-dimensional representations of the trajectory using only the information encoded in the trajectory. The downside of this approach is that different sets of coordinates are required for each different chemical systems under study precisely because the coordinates are constructed using information from the trajectory. In this paper we show how one can resolve this problem by using the sketch-map algorithm that we recently proposed to construct a low-dimensional representation of the structures contained in the protein data bank (PDB). We show that the resulting coordinates are as useful for analysing trajectory data as coordinates constructed using landmark configurations taken from the trajectory and that these coordinates can thus be used for understanding protein folding across a range of systems.

The Validation of a Two-Dimensional Transient Catalyst Model for Direct Injection Two-Stroke Applications

This paper describes the detailed validation of a computer model designed to simulate the transient light-off in a two-stroke oxidation catalyst. A plug flow reactor is employed to provide measurements of temperature and gas concentration at various radial and axial locations inside the catalyst. These measurements are recorded at discrete intervals during a transient light-off in which the inlet temperature is increased from ambient to 300oC at rates of up to 6oC/sec. The catalyst formulation used in the flow reactor, and its associated test procedures, are then simulated by the computer and a comparison made between experimental readings and model predictions. The design of the computer model to which this validation exercise relates is described in detail in a separate technical paper. The first section of the paper investigates the warm-up characteristics of the substrate and examines the validity of the heat transfer predictions between the wall and the gas in the absence of chemical reactions. The predictions from a typical single-component CO transient light-off test are discussed in the second section and are compared with experimental data. In particular the effect of the temperature ramp on the light-off curve and reaction zone development is examined. An analysis of the C3H6 conversion is given in the third section while the final section examines the accuracy of the light-off curves which are produced when both CO and C3H6 are present in the feed gas. The analysis shows that the heat and mass transfer calculations provided reliable predictions of the warm-up behaviour and post light-off gas concentration profiles. The self-inhibition and cross-inhibition terms in the global rate expressions were also found to be reasonably reliable although the surface reaction rates required calibration with experimental data.

Modelling the Geometry of the Repeat Unit Cell of Three-Dimensional Weave Architectures,

The three-dimensional (3D) weaving process offers the ability to tailor the mechanical properties via design of the weave architecture. One repeat of the 3D woven fabric is represented by the unit cell. The model accepts basic weaver and material manufacturer data as inputs in order to calculate the geometric characteristics of the 3D woven unit cell. The specific weave architecture manufactured and subsequently modelled had an angle interlock type binding configuration. The modelled result was shown to have a close approximation compared to the experimentally measured values and highlighted the importance of the representation of the binder tow path.

Three dimensional morphology and compressive behaviour of sintered biodegradable composite scaffolds

Porous poly-L-lactide acid (PLA) scaffolds are prepared using polymer sintering and porogen leaching method. Different weight fractions of the Hydroxyapatite (HA) are added to the PLA to control the acidity and degradation rate. The three dimensional morphology and surface porosity are tested using micro CT, optical microscopy and scanning electron microscopy (SEM). Results indicate that the surface porosity does not change by addition of HA. The micro Ct examinations show slight decrease in the pore size and increase in wall thickness accompanied with reduced anisotropy for the scaffolds containing HA. SEM micrographs show detectable interconnected pores for the scaffold with pure PLA. Addition of the HA results in agglomeration of the HA which blocks some of the pores. Compression tests of the scaffold identify three stages in the stress-strain curve. The addition of HA adversely affects the modulus of the scaffold at the first stage, but this was reversed for the second and third stages of the compression. The results of these tests are compared with the cellular material model. The manufactured scaffold have acceptable properties for a scaffold, however improvement to the mixing of the phases of PLA and HA is required to achieve better integrity of the composite scaffolds.

Explaining the multiparty systems of Western Europe prior to the adoption of proportional representation

The literature has difficulty explaining why the number of parties in majoritarian electoral systems often exceeds the two-party predictions associated with Duverger’s Law. To understand why this is the case, I examine several party systems in Western Europe before the adoption of proportional representation. Drawing from the social cleavage approach, I argue that the emergence of multiparty systems was because of the development of the class cleavage, which provided a base of voters sizeable enough to support third parties. However, in countries where the class cleavage became the largest cleavage, the class divide displaced other cleavages and the number of parties began to converge on two. The results show that the effect of the class cleavage was nonlinear, producing the greatest party system fragmentation in countries where class cleavages were present – but not dominant – and smaller in countries where class cleavages were either dominant or non-existent.

Optimizing geometric accuracy of four dimensional CT scans acquired using the wall and couch mounted Varian Real-Time Position Management camera systems

Objective:

The aim of this study was to identify sources of anatomical misrepresentation due to the location of camera mounting, tumour motion velocity and image processing artefacts in order to optimise the 4DCT scan protocol and improve geometrical-temporal accuracy.

A phantom with an imaging insert was driven with a sinusoidal superior-inferior motion of varying amplitude and period for 4DCT scanning. The length of a high density cube within the insert was measured using treatment planning software to determine the accuracy of its spatial representation. Scan parameters were varied including the tube rotation period and the cine time between reconstructed images. A CT image quality phantom was used to measure various image quality signatures under the scan parameters tested.

No significant difference in spatial accuracy was found for 4DCT scans carried out using the wall mounted or couch mounted camera for sinusoidal target motion. Greater spatial accuracy was found for 4DCT scans carried out using a tube rotation speed of 0.5s rather than 1.0s. The reduction in image quality when using a faster rotation speed was not enough to require an increase in patient dose.

4DCT accuracy may be increased by optimising scan parameters, including choosing faster tube rotation speeds. Peak misidentification in the recorded breathing trace leads to spatial artefacts and this risk can be reduced by using a couch mounted infrared camera.

This study explicitly shows that 4DCT scan accuracy is improved by scanning with a faster CT tube rotation speed.

Investigation of tanpura string vibrations using a two-dimensional time-domain model incorporating coupling and bridge friction

Tanpura string vibrations have been investigated previously using numerical models based on energy conserving schemes derived from a Hamiltonian description in one-dimensional form. Such time-domain models have the property that, for the lossless case, the numerical Hamiltonian (representing total energy of the system) can be proven to be constant from one time step
to the next, irrespective of any of the system parameters; in practice the Hamiltonian can be shown to be conserved within machine precision. Models of this kind can reproduce a jvari effect, which results from the bridge-string interaction. However the one-dimensional formulation has recently been shown to fail to replicate the jvaris strong dependence on the thread placement. As a first step towards simulations which accurately emulate this sensitivity to the thread placement, a twodimensional model is proposed, incorporating coupling of controllable level between the two string polarisations at the string termination opposite from the barrier. In addition, a friction force acting when the string slides across the bridge in horizontal direction is introduced, thus effecting a further damping mechanism. In this preliminary study, the string is terminated at the position of the thread. As in the one-dimensional model, an implicit scheme has to be used to solve the system, employing Newton's method to calculate the updated positions and momentums of each string segment. The two-dimensional model is proven to be energy conserving when the loss parameters are set to zero, irrespective of the coupling constant. Both frequency-dependent and independent losses are then added to the string, so that the model can be compared to analogous instruments. The influence of coupling and the bridge friction are investigated.

Time-resolved four-wave-mixing spectroscopy for inner-valence transitions

Noncollinear four-wave-mixing (FWM) techniques at near-infrared (NIR), visible, and ultraviolet frequencies have been widely used to map vibrational and electronic couplings, typically in complex molecules. However, correlations between spatially localized inner-valence transitions among different sites of a molecule in the extreme ultraviolet (XUV) spectral range have not been observed yet. As an experimental step toward this goal, we perform time-resolved FWM spectroscopy with femtosecond NIR and attosecond XUV pulses. The first two pulses (XUV-NIR) coincide in time and act as coherent excitation fields, while the third pulse (NIR) acts as a probe. As a first application, we show how coupling dynamics between odd- and even-parity, inner-valence excited states of neon can be revealed using a two-dimensional spectral representation. Experimentally obtained results are found to be in good agreement with ab initio time-dependent R-matrix calculations providing the full description of multielectron interactions, as well as few-level model simulations. Future applications of this method also include site-specific probing of electronic processes in molecules.