Electrodynamic heart model construction and ECG simulation

Objectives: In this paper, we present a unified electrodynamic heart model that permits simulations of the body surface potentials generated by the heart in motion. The inclusion of motion in the heart model significantly improves the accuracy of the simulated body surface potentials and therefore also the 12-lead ECG. Methods: The key step is to construct an electromechanical heart model. The cardiac excitation propagation is simulated by an electrical heart model, and the resulting cardiac active forces are used to calculate the ventricular wall motion based on a mechanical model. The source-field point relative position changes during heart systole and diastole. These can be obtained, and then used to calculate body surface ECG based on the electrical heart-torso model. Results: An electromechanical biventricular heart model is constructed and a standard 12-lead ECG is simulated. Compared with a simulated ECG based on the static electrical heart model, the simulated ECG based on the dynamic heart model is more accordant with a clinically recorded ECG, especially for the ST segment and T wave of a V1-V6 lead ECG. For slight-degree myocardial ischemia ECG simulation, the ST segment and T wave changes can be observed from the simulated ECG based on a dynamic heart model, while the ST segment and T wave of simulated ECG based on a static heart model is almost unchanged when compared with a normal ECG. Conclusions: This study confirms the importance of the mechanical factor in the ECG simulation. The dynamic heart model could provide more accurate ECG simulation, especially for myocardial ischemia or infarction simulation, since the main ECG changes occur at the ST segment and T wave, which correspond with cardiac systole and diastole phases.

Mechanical properties of interfacial films formed by lysozyme self-assembly at the air-water interface

We present the first characterization of the mechanical properties of lysozyme films formed by self-assembly at the air-water interface using the Cambridge interfacial tensiometer (CIT), an apparatus capable of subjecting protein films to a much higher level of extensional strain than traditional dilatational techniques. CIT analysis, which is insensitive to surface pressure, provides a direct measure of the extensional stress-strain behavior of an interfacial film without the need to assume a mechanical model (e.g., viscoelastic), and without requiring difficult-to-test assumptions regarding low-strain material linearity. This testing method has revealed that the bulk solution pH from which assembly of an interfacial lysozyme film occurs influences the mechanical properties of the film more significantly than is suggested by the observed differences in elastic moduli or surface pressure. We have also identified a previously undescribed pH dependency in the effect of solution ionic strength on the mechanical strength of the lysozyme films formed at the air-water interface. Increasing solution ionic strength was found to increase lysozyme film strength when assembly occurred at pH 7, but it caused a decrease in film strength at pH 11, close to the pI of lysozyme. This result is discussed in terms of the significant contribution made to protein film strength by both electrostatic interactions and the hydrophobic effect. Washout experiments to remove protein from the bulk phase have shown that a small percentage of the interfacially adsorbed lysozyme molecules are reversibly adsorbed. Finally, the washout tests have probed the role played by additional adsorption to the fresh interface formed by the application of a large strain to the lysozyme film and have suggested the movement of reversibly bound lysozyme molecules from a subinterfacial layer to the interface.

Computer simulation of radio-frequency methane/hydrogen plasmas and their interaction with GaAs Surfaces

The following thesis describes the computer modelling of radio frequency capacitively coupled methane/hydrogen plasmas and the consequences for the reactive ion etching of (100) GaAs surfaces. In addition a range of etching experiments was undertaken over a matrix of pressure, power and methane concentration. The resulting surfaces were investigated using X-ray photoelectron spectroscopy and the results were discussed in terms of physical and chemical models of particle/surface interactions in addition to the predictions for energies, angles and relative fluxes to the substrate of the various plasma species. The model consisted of a Monte Carlo code which followed electrons and ions through the plasma and sheath potentials whilst taking account of collisions with background neutral gas molecules. The ionisation profile output from the electron module was used as input for the ionic module. Momentum scattering interactions of ions with gas molecules were investigated via different models and compared against results given by quantum mechanical code. The interactions were treated as central potential scattering events and the resulting neutral cascades were followed. The resulting predictions for ion energies at the cathode compared well to experimental ion energy distributions and this verified the particular form of the electrical potentials used and their applicability in the particular geometry plasma cell used in the etching experiments. The final code was used to investigate the effect of external plasma parameters on the mass distribution, energy and angles of all species impingent on the electrodes. Comparisons of electron energies in the plasma also agreed favourably with measurements made using a Langmuir electric probe. The surface analysis showed the surfaces all to be depleted in arsenic due to its preferential removal and the resultant Ga:As ratio in the surface was found to be directly linked to the etch rate. The etch rate was determined by the methane flux which was predicted by the code.

Theoretical study of chloroperoxidase catalyzed chlorination of beta-cyclopentanedione and role of water in the chlorination mechanism

Chloroperoxidase (CPO) is a potential biocatalyst for use in asymmetric synthesis. The mechanisms of CPO catalysis are therefore of interest. The halogenation reaction, one of several chemical reactions that CPO catalyzes, is not fully understood and is the subject of this dissertation. The mechanism by which CPO catalyzes halogenation is disputed. It has been postulated that halogenation of substrates occurs at the active site. Alternatively, it has been proposed that hypochlorous acid, produced at the active site via oxidation of chloride, is released prior to reaction, so that halogenation occurs in solution. The free-solution mechanism is supported by the observation that halogenation of most substrates often occurs non-stereospecifically. On the other hand, the enzyme-bound mechanism is supported by the observation that some large substrates undergo halogenation stereospecifically. The major purpose of this research is to compare chlorination of the substrate β-cyclopentanedione in the two environments. One study was of the reaction with limited hydration because such a level of hydration is typical of the active site. For this work, a purely quantum mechanical approach was used. To model the aqueous environment, the limited hydration environment approach is not appropriate. Instead, reaction precursor conformations were obtained from a solvated molecular dynamics simulation, and reaction of potentially reactive molecular encounters was modeled with a hybrid quantum mechanical/molecular mechanical approach. Extensive work developing parameters for small molecules was pre-requisite for the molecular dynamics simulation. It is observed that a limited and optimized (active-site-like) hydration environment leads to a lower energetic barrier than the fully solvated model representative of the aqueous environment at room temperature, suggesting that the stable water network near the active site is likely to facilitate the chlorination mechanism. The influence of the solvent environment on the reaction barrier is critical. It is observed that stabilization of the catalytic water by other solvent molecules lowers the barrier for keto-enol tautomerization. Placement of water molecules is more important than the number of water molecules in such studies. The fully-solvated model demonstrates that reaction proceeds when the instantaneous dynamical water environment is close to optimal for stabilizing the transition state.

Nonlinear quantum transport in low-dimensional electronic devices

The study of transport processes in low-dimensional semiconductors requires a rigorous quantum mechanical treatment. However, a full-fledged quantum transport theory of electrons (or holes) in semiconductors of small scale, applicable in the presence of external fields of arbitrary strength, is still not available. In the literature, different approaches have been proposed, including: (a) the semiclassical Boltzmann equation, (b) perturbation theory based on Keldysh's Green functions, and (c) the Quantum Boltzmann Equation (QBE), previously derived by Van Vliet and coworkers, applicable in the realm of Kubo's Linear Response Theory (LRT). ^ In the present work, we follow the method originally proposed by Van Wet in LRT. The Hamiltonian in this approach is of the form: H = H 0(E, B) + λV, where H0 contains the externally applied fields, and λV includes many-body interactions. This Hamiltonian differs from the LRT Hamiltonian, H = H0 - AF(t) + λV, which contains the external field in the field-response part, -AF(t). For the nonlinear problem, the eigenfunctions of the system Hamiltonian, H0(E, B), include the external fields without any limitation on strength. ^ In Part A of this dissertation, both the diagonal and nondiagonal Master equations are obtained after applying projection operators to the von Neumann equation for the density operator in the interaction picture, and taking the Van Hove limit, (λ → 0, t → ∞, so that (λ2 t)n remains finite). Similarly, the many-body current operator J is obtained from the Heisenberg equation of motion. ^ In Part B, the Quantum Boltzmann Equation is obtained in the occupation-number representation for an electron gas, interacting with phonons or impurities. On the one-body level, the current operator obtained in Part A leads to the Generalized Calecki current for electric and magnetic fields of arbitrary strength. Furthermore, in this part, the LRT results for the current and conductance are recovered in the limit of small electric fields. ^ In Part C, we apply the above results to the study of both linear and nonlinear longitudinal magneto-conductance in quasi one-dimensional quantum wires (1D QW). We have thus been able to quantitatively explain the experimental results, recently published by C. Brick, et al., on these novel frontier-type devices. ^

Nonlinear quantum transport in low-dimensional electronic devices

The study of transport processes in low-dimensional semiconductors requires a rigorous quantum mechanical treatment. However, a full-fledged quantum transport theory of electrons (or holes) in semiconductors of small scale, applicable in the presence of external fields of arbitrary strength, is still not available. In the literature, different approaches have been proposed, including: (a) the semiclassical Boltzmann equation, (b) perturbation theory based on Keldysh's Green functions, and (c) the Quantum Boltzmann Equation (QBE), previously derived by Van Vliet and coworkers, applicable in the realm of Kubo's Linear Response Theory (LRT). In the present work, we follow the method originally proposed by Van Vliet in LRT. The Hamiltonian in this approach is of the form: H = H°(E, B) + λV, where H0 contains the externally applied fields, and λV includes many-body interactions. This Hamiltonian differs from the LRT Hamiltonian, H = H° - AF(t) + λV, which contains the external field in the field-response part, -AF(t). For the nonlinear problem, the eigenfunctions of the system Hamiltonian, H°(E, B) , include the external fields without any limitation on strength. In Part A of this dissertation, both the diagonal and nondiagonal Master equations are obtained after applying projection operators to the von Neumann equation for the density operator in the interaction picture, and taking the Van Hove limit, (λ → 0 , t → ∞ , so that (λ2 t)n remains finite). Similarly, the many-body current operator J is obtained from the Heisenberg equation of motion. In Part B, the Quantum Boltzmann Equation is obtained in the occupation-number representation for an electron gas, interacting with phonons or impurities. On the one-body level, the current operator obtained in Part A leads to the Generalized Calecki current for electric and magnetic fields of arbitrary strength. Furthermore, in this part, the LRT results for the current and conductance are recovered in the limit of small electric fields. In Part C, we apply the above results to the study of both linear and nonlinear longitudinal magneto-conductance in quasi one-dimensional quantum wires (1D QW). We have thus been able to quantitatively explain the experimental results, recently published by C. Brick, et al., on these novel frontier-type devices.

Multiple sulfur isotope fractionations in inorganic aqueous systems

New constraints on isotope fractionation factors in inorganic aqueous sulfur systems based on theoretical and experimental techniques relevant to studies of the sulfur cycle in modern environments and the geologic rock record are presented in this dissertation. These include theoretical estimations of equilibrium isotope fractionation factors utilizing quantum mechanical software and a water cluster model approach for aqueous sulfur compounds that span the entire range of oxidation state for sulfur. These theoretical calculations generally reproduce the available experimental determinations from the literature and provide new constraints where no others are available. These theoretical calculations illustrate in detail the relationship between sulfur bonding environment and the mass dependence associated with equilibrium isotope exchange reactions involving all four isotopes of sulfur. I additionally highlight the effect of isomers of protonated compounds (compounds with the same chemical formula but different structure, where protons are bound to either sulfur or oxygen atoms) on isotope partitioning in the sulfite (S4+) and sulfoxylate (S2+) systems, both of which are key intermediates in oxidation-reduction processes in the sulfur cycle. I demonstrate that isomers containing the highest degree of coordination around sulfur (where protonation occurs on the sulfur atom) have a strong influence on isotopic fractionation factors, and argue that isomerization phenomenon should be considered in models of the sulfur cycle. Additionally, experimental results of the reaction rates and isotope fractionations associated with the chemical oxidation of aqueous sulfide are presented. Sulfide oxidation is a major process in the global sulfur cycle due largely to the sulfide-producing activity of anaerobic microorganisms in organic-rich marine sediments. These experiments reveal relationships between isotope fractionations and reaction rate as a function of both temperature and trace metal (ferrous iron) catalysis that I interpret in the context of the complex mechanism of sulfide oxidation. I also demonstrate that sulfide oxidation is a process associated with a mass dependence that can be described as not conforming to the mass dependence typically associated with equilibrium isotope exchange. This observation has implications for the inclusion of oxidative processes in environmental- and global-scale models of the sulfur cycle based on the mass balance of all four isotopes of sulfur. The contents of this dissertation provide key reference information on isotopic fractionation factors in aqueous sulfur systems that will have far-reaching applicability to studies of the sulfur cycle in a wide variety of natural settings.

Circuit Quantum Electrodynamics with Transmon Qubits in the Ultrastrong Coupling Regime

68 pg.

Multi-objective optimization of material model parameters of an adhesive layer by using SPEA2

The usage of multi material structures in industry, especially in the automotive industry are increasing. To overcome the difficulties in joining these structures, adhesives have several benefits over traditional joining methods. Therefore, accurate simulations of the entire process of fracture including the adhesive layer is crucial. In this paper, material parameters of a previously developed meso mechanical finite element (FE) model of a thin adhesive layer are optimized using the Strength Pareto Evolutionary Algorithm (SPEA2). Objective functions are defined as the error between experimental data and simulation data. The experimental data is provided by previously performed experiments where an adhesive layer was loaded in monotonically increasing peel and shear. Two objective functions are dependent on 9 model parameters (decision variables) in total and are evaluated by running two FEsimulations, one is loading the adhesive layer in peel and the other in shear. The original study converted the two objective functions into one function that resulted in one optimal solution. In this study, however, a Pareto frontis obtained by employing the SPEA2 algorithm. Thus, more insight into the material model, objective functions, optimal solutions and decision space is acquired using the Pareto front. We compare the results and show good agreement with the experimental data.

Modelling word activation in semantic networks : three scaled entanglement models compared

Modelling how a word is activated in human memory is an important requirement for determining the probability of recall of a word in an extra-list cueing experiment. Previous research assumed a quantum-like model in which the semantic network was modelled as entangled qubits, however the level of activation was clearly being over-estimated. This paper explores three variations of this model, each of which are distinguished by a scaling factor designed to compensate the overestimation.

Perceptions of document relevance

This article presents a study of how humans perceive and judge the relevance of documents. Humans are adept at making reasonably robust and quick decisions about what information is relevant to them, despite the ever increasing complexity and volume of their surrounding information environment. The literature on document relevance has identified various dimensions of relevance (e.g., topicality, novelty, etc.), however little is understood about how these dimensions may interact. We performed a crowdsourced study of how human subjects judge two relevance dimensions in relation to document snippets retrieved from an internet search engine. The order of the judgment was controlled. For those judgments exhibiting an order effect, a q–test was performed to determine whether the order effects can be explained by a quantum decision model based on incompatible decision perspectives. Some evidence of incompatibility was found which suggests incompatible decision perspectives is appropriate for explaining interacting dimensions of relevance in such instances.

Towards a multi-element simulation of heavily over-consolidated bagasse using a material subroutine

A better understanding of the behaviour of prepared cane and bagasse, and the ability to model the mechanical behaviour of bagasse as it is squeezed in a milling unit to extract juice, would help identify how to improve the current process. There are opportunities to decrease bagasse moisture from a milling unit. The behaviour of bagasse in chutes is poorly understood. Previous investigations have shown that juice flow through bagasse obeys Darcy’s permeability law, that the grip of the rough surface of the grooves on the bagasse can be represented by the Mohr-Coulomb failure criterion for soils, and that the internal mechanical behaviour of the bagasse is critical state behaviour similar to that for sand and clay. Progress has been made in the last 11 years towards implementing a mechanical model for bagasse in finite element software. The objective is to be able to correctly simulate various simple mechanical loading conditions measured in the laboratory. Combining these behaviours together is thought to have a high probability of reproducing the complicated stress conditions in a milling unit. This paper reports on progress made towards modelling the fifth and final (and most challenging) of the simple loading conditions: the shearing of heavily over-consolidated bagasse, using a specific model for bagasse in a multi-element simulation.

Determining the material properties for heavily over-consolidated bagasse through parameter estimation

A better understanding of the behaviour of prepared cane and bagasse, and the ability to model the mechanical behaviour of bagasse as it is squeezed in a milling unit to extract juice, would help identify how to improve the current process. For example, there are opportunities to decrease bagasse moisture from a milling unit. Also, the behaviour of bagasse in chutes is poorly understood. Previous investigations have shown that juice flow through bagasse obeys Darcy’s permeability law, that the grip of the rough surface of the grooves on the bagasse can be represented by the Mohr-Coulomb failure criterion for soils, and that the internal mechanical behaviour of the bagasse is critical state behaviour similar to that for sand and clay. Progress has been made in the last ten years towards implementing a mechanical model for bagasse in finite element software. The objective has been to be able to simulate simple mechanical loading conditions measured in the laboratory, which, when combined together, have a high probability of reproducing the complicated stress conditions in a milling unit. This paper reports on the successful simulation of part of the fifth and final (and most challenging) loading condition, the shearing of heavily over-consolidated bagasse, and determining material property values through the use of powerful and free parameter estimation software.

Modelling attitudes to climate change : an order effect and a test between alternatives

Quantum-like models can be fruitfully used to model attitude change in a social context. Next steps require data, and higher dimensional models. Here, we discuss an exploratory study that demonstrates an order effect when three question sets about Climate Beliefs, Political Affiliation and Attitudes Towards Science are presented in different orders within a larger study of n=533 subjects. A quantum-like model seems possible, and we propose a new experiment which could be used to test between three possible models for this scenario.

Resonance raman intensity analysis of polyatomic molecules

A time-dependent quantum mechanical (TDQM) method of wavepacket propagation in computing resonance Raman intensities for polyatomic systems, has been developed and demonstrated by applying it tocis-stilbene andtrans-azobenzene. In the case of the former, Raman excitation profiles (REPs) for the various vibrational modes have also been computed. It is observed that the calculated absorption spectrum and the REPs compare very well with the experimental results. A comparison of these results with those of the often semiclassical approach reveals that the TDQM method can be used to study polyatomic systems with as much ease as the semiclassical wavepacket method.