Modern Computational Chemistry Methods for Prediction of Ground- and Excited-State Properties in Open-Shell Systems

The main goal of the research presented in this work is to provide some important insights about computational modeling of open-shell species. Such projects are: the investigation of the size-extensivity error in Equation-of-Motion Coupled Cluster methods, the analysis of the Long-Range corrected scheme in predicting UV-Vis spectra of Cu(II) complexes with the 4-imidazole acetate and its ethylated derivative, and the exploration of the importance of choosing a proper basis set for the description of systems such as the lithium monoxide anion. The most significant findings of this research are: (i) The contribution of the left operator to the size-extensivity error of the CR-EOMCC(2,3) approach, (ii) The cause of d-d shifts when varying the range-separation parameter and the amount of the exact exchange arising from the imbalanced treatment of localized vs. delocalized orbitals via the "tuned" CAM-B3LYP* functional, (iii) The proper acidity trend of the first-row hydrides and their lithiated analogs that may be reversed if the basis sets are not correctly selected.

Investigation of phase separation in bulk heterojunction solar cells via supramolecular chemistry

In this work, we have prepared two donor-acceptor-donor (D-A-D) pi-conjugated oligomers to investigate the effect of phase separation on the performance of bulk heterojunction (BHJ) solar cells. These charge transfer low band gap pi-conjugated oligomers (TTB and NMeTTB) were synthesized by Knoevenagel condensation of terthiophenecarbaldehyde and barbiturate appended pyran derivative. The thin film morphology of both the oligomers and along with electron acceptor 6,6]-phenyl-C60-butyric acid methyl ester (PC61BM) was investigated by atomic force microscopy (AFM) and transmission electron microscopy (TEM). The blend of NMeTTB and PC61BM thin film yield highly ordered thin film, whereas there was clear phase separation between TTB and PC61BM in thin film. The BHJ solar cell was fabricated using a blend of NMeTTB and TTB with PC61BM acceptor in 1:1 ratio as active layer, and a power conversion efficiency of 1.8% was obtained. This device characteristic was compared with device having TTB:PC61BM as active layer, and large difference is observed in photocurrents. This poor performance of TTB in BHJ devices was attributed to the difference in the nanoscale morphology of the corresponding derivatives. We rationalize our findings based on the low charge carrier mobility in organic field-effect transistors and miscibility/phase separation parameter of binary components (oligomers and PC61BM) in the active layer of bulk heterojunction solar cells.

Comparing methods for modelling spreading cell fronts

Spreading cell fronts play an essential role in many physiological processes. Classically, models of this process are based on the Fisher-Kolmogorov equation; however, such continuum representations are not always suitable as they do not explicitly represent behaviour at the level of individual cells. Additionally, many models examine only the large time asymptotic behaviour, where a travelling wave front with a constant speed has been established. Many experiments, such as a scratch assay, never display this asymptotic behaviour, and in these cases the transient behaviour must be taken into account. We examine the transient and asymptotic behaviour of moving cell fronts using techniques that go beyond the continuum approximation via a volume-excluding birth-migration process on a regular one-dimensional lattice. We approximate the averaged discrete results using three methods: (i) mean-field, (ii) pair-wise, and (iii) one-hole approximations. We discuss the performace of these methods, in comparison to the averaged discrete results, for a range of parameter space, examining both the transient and asymptotic behaviours. The one-hole approximation, based on techniques from statistical physics, is not capable of predicting transient behaviour but provides excellent agreement with the asymptotic behaviour of the averaged discrete results, provided that cells are proliferating fast enough relative to their rate of migration. The mean-field and pair-wise approximations give indistinguishable asymptotic results, which agree with the averaged discrete results when cells are migrating much more rapidly than they are proliferating. The pair-wise approximation performs better in the transient region than does the mean-field, despite having the same asymptotic behaviour. Our results show that each approximation only works in specific situations, thus we must be careful to use a suitable approximation for a given system, otherwise inaccurate predictions could be made.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Outbreak detection algorithms for seasonal disease data : a case study using Ross River virus disease

Background Detection of outbreaks is an important part of disease surveillance. Although many algorithms have been designed for detecting outbreaks, few have been specifically assessed against diseases that have distinct seasonal incidence patterns, such as those caused by vector-borne pathogens. Methods We applied five previously reported outbreak detection algorithms to Ross River virus (RRV) disease data (1991-2007) for the four local government areas (LGAs) of Brisbane, Emerald, Redland and Townsville in Queensland, Australia. The methods used were the Early Aberration Reporting System (EARS) C1, C2 and C3 methods, negative binomial cusum (NBC), historical limits method (HLM), Poisson outbreak detection (POD) method and the purely temporal SaTScan analysis. Seasonally-adjusted variants of the NBC and SaTScan methods were developed. Some of the algorithms were applied using a range of parameter values, resulting in 17 variants of the five algorithms. Results The 9,188 RRV disease notifications that occurred in the four selected regions over the study period showed marked seasonality, which adversely affected the performance of some of the outbreak detection algorithms. Most of the methods examined were able to detect the same major events. The exception was the seasonally-adjusted NBC methods that detected an excess of short signals. The NBC, POD and temporal SaTScan algorithms were the only methods that consistently had high true positive rates and low false positive and false negative rates across the four study areas. The timeliness of outbreak signals generated by each method was also compared but there was no consistency across outbreaks and LGAs. Conclusions This study has highlighted several issues associated with applying outbreak detection algorithms to seasonal disease data. In lieu of a true gold standard, a quantitative comparison is difficult and caution should be taken when interpreting the true positives, false positives, sensitivity and specificity.

Choosing an appropriate modelling framework for analysing multispecies co-culture cell biology experiments

In vitro cell biology assays play a crucial role in informing our understanding of the migratory, proliferative and invasive properties of many cell types in different biological contexts. While mono-culture assays involve the study of a population of cells composed of a single cell type, co-culture assays study a population of cells composed of multiple cell types (or subpopulations of cells). Such co-culture assays can provide more realistic insights into many biological processes including tissue repair, tissue regeneration and malignant spreading. Typically, system parameters, such as motility and proliferation rates, are estimated by calibrating a mathematical or computational model to the observed experimental data. However, parameter estimates can be highly sensitive to the choice of model and modelling framework. This observation motivates us to consider the fundamental question of how we can best choose a model to facilitate accurate parameter estimation for a particular assay. In this work we describe three mathematical models of mono-culture and co-culture assays that include different levels of spatial detail. We study various spatial summary statistics to explore if they can be used to distinguish between the suitability of each model over a range of parameter space. Our results for mono-culture experiments are promising, in that we suggest two spatial statistics that can be used to direct model choice. However, co-culture experiments are far more challenging: we show that these same spatial statistics which provide useful insight into mono-culture systems are insuffcient for co-culture systems. Therefore, we conclude that great care ought to be exercised when estimating the parameters of co-culture assays.

Modelling the movement of interacting cell populations: A moment dynamics approach

Mathematical models describing the movement of multiple interacting subpopulations are relevant to many biological and ecological processes. Standard mean-field partial differential equation descriptions of these processes suffer from the limitation that they implicitly neglect to incorporate the impact of spatial correlations and clustering. To overcome this, we derive a moment dynamics description of a discrete stochastic process which describes the spreading of distinct interacting subpopulations. In particular, we motivate our model by mimicking the geometry of two typical cell biology experiments. Comparing the performance of the moment dynamics model with a traditional mean-field model confirms that the moment dynamics approach always outperforms the traditional mean-field approach. To provide more general insight we summarise the performance of the moment dynamics model and the traditional mean-field model over a wide range of parameter regimes. These results help distinguish between those situations where spatial correlation effects are sufficiently strong, such that a moment dynamics model is required, from other situations where spatial correlation effects are sufficiently weak, such that a traditional mean-field model is adequate.

Hole pairing within an extended Anderson impurity model applicable to the high-Tc cuprates

We calculate the binding energy of a hole pair within the extended Anderson Hamiltonian for the high-Tc cuprates including a Cu impurity and an oxygen-derived band. The results indicate that stable hole pairs can be formed for intra-atomic and interatomic Coulomb repulsion strengths larger than 6 and 3.5 eV, respectively. It is also shown that the total hybridization strength between the Cu 3d and oxygen p band should be less than 2.5 eV. The hole pairing takes place primarily within the oxygen-derived p band. The range of parameter values for which hole pairing occurs is also consistent with the earlier photoemission results from these cuprates.

Analysis of an Aquifer-water table aquitard system

The problem of pumping an aquifer in an aquifer-water table aquitard system is considered, accounting for the elastic properties of both the aquifer and the aquitard, the gravity drainage in the aquitard and treating the water table as an unknown boundary. The coupled partial differential equations are nondimensionalised, yielding three principal parameters governing the problem. The numerical solution of these equations is obtained for a wide range of parameter values. Type curves are generated and their use is illustrated through a field application.

Transposable element-medicated evolution of sex: A population genetic model

For a population made up of individuals capable of sexual as well as asexual modes of reproduction, conditions for the spread of a transposable element are explored using a one-locus, two-haplotype model. The analysis is then extended to include the possibility that the transposable element can modulate the probability of sexual reproduction, thus casting Hickey’s (1982,Genetics 101: 519–531) suggestion in a population genetics framework. The model explicitly includes the cost of sexual reproduction, fitness disadvantage to the transposable element, probability of transposition, and the predisposition for sexual reproduction in the presence and absence of the transposable element. The model predicts several kinds of outcome, including initial frequency dependence and stable polymorphism. More importantly, it is seen that for a wide range of parameter values, the transposable element can go to fixation. Therefore it is able to convert the population from a predominantly asexual to a predominantly sexual mode of reproduction. Viewed in conjunction with recent results implicating short stretches of apparently non-coding DNA in sex determination (McCoubreyet al. 1988,Science 242: 1146–1151), the model hints at the important role this mechanism could have played in the evolution of sexuality.

Selective Flotation of Molybdenite and Chalcopyrite

Separation of molybdenite from chalcopyrite from the bulk concentrate of a low grade molybdenite ore from Chintamani, Kolar District, Karnataka, has been attempted in a modified Hallimond tube microflotation set-up. The flotation parameters studied are, sodium sulphide and potassium cyanide as depressants for chalcopyrite, pH, depressant concentration, steaming time and time of flotation. Selectivity index was used as a separation parameter. Steaming for 25 min followed by flotation for 7 min at 80 mg/l concentration of potassium cyanide and pH 9 gave the best results for separation of molybdenite. 19 ref.--AA

Coupled Wavenumbers in an Infinite Flexible Fluid-Filled Circular Cylindrical Shell: Comparison between Different Shell Theories

Analytical expressions are found for the wavenumbers in an infinite flexible in vacuo I fluid-filled circular cylindrical shell based on different shell-theories using asymptotic methods. Donnell-Mushtari theory (the simplest shell theory) and four higher order theories, namely Love-Timoshenko, Goldenveizer-Novozhilov, Flugge and Kennard-simplified are considered. Initially, in vacuo and fluid-coupled wavenumber expressions are presented using the Donnell-Mushtari theory. Subsequently, the wavenumbers using the higher order theories are presented as perturbations on the Donnell-Mushtari wavenumbers. Similarly, expressions for the resonance frequencies in a finite shell are also presented, using each shell theory. The basic differences between the theories being what they are, the analytical expressions obtained from the five theories allow one to see how these differences propagate into the asymptotic expansions. Also, they help to quantify the difference between the theories for a wide range of parameter values such as the frequency range, circumferential order, thickness ratio of the shell, etc.

Bose-Hubbard models in confining potentials: Inhomogeneous mean-field theory

We present an extensive study of Mott insulator (MI) and superfluid (SF) shells in Bose-Hubbard (BH) models for bosons in optical lattices with harmonic traps. For this we apply the inhomogeneous mean-field theory developed by Sheshadri et al. Phys. Rev. Lett. 75, 4075 (1995)]. Our results for the BH model with one type of spinless bosons agree quantitatively with quantum Monte Carlo simulations. Our approach is numerically less intensive than such simulations, so we are able to perform calculations on experimentally realistic, large three-dimensional systems, explore a wide range of parameter values, and make direct contact with a variety of experimental measurements. We also extend our inhomogeneous mean-field theory to study BH models with harmonic traps and (a) two species of bosons or (b) spin-1 bosons. With two species of bosons, we obtain rich phase diagrams with a variety of SF and MI phases and associated shells when we include a quadratic confining potential. For the spin-1 BH model, we show, in a representative case, that the system can display alternating shells of polar SF and MI phases, and we make interesting predictions for experiments in such systems.

Status of supersymmetric type-I seesaw in SO(10) inspired models

We report on the status of supersymmetric seesaw models in the light of recent experimental results on mu -> e + gamma, theta(13) and the light Higgs mass at the LHC. SO(10)-like relations are assumed for neutrino Dirac Yukawa couplings and two cases of mixing, one large, PMNS-like, and another small, CKM-like, are considered. It is shown that for the large mixing case, only a small range of parameter space with moderate tan beta is still allowed. This remaining region can be ruled out by an order of magnitude improvement in the current limit on BR(mu -> e + gamma). We also explore a model with non-universal Higgs mass boundary conditions at the high scale. It is shown that the renormalization group induced flavor violating slepton mass terms are highly sensitive to the Higgs boundary conditions. Depending on the choice of the parameters, they can either lead to strong enhancements or cancellations within the flavor violating terms. Such cancellations might relax the severe constraints imposed by lepton flavor violation compared to mSUGRA. Nevertheless for a large region of parameter space the predicted rates lie within the reach of future experiments once the light Higgs mass constraint is imposed. We also update the potential of the ongoing and future experimental searches for lepton flavor violation in constraining the supersymmetric parameter space.

The role of intra- and inter-site exchange correlations in the extended Falicov-Kimball model on a triangular lattice

Ground state magnetic properties of the spin-dependent Falicov-Kimball model (FKM) are studied by incorporating the intrasite exchange correlation J (between itinerant d- and localized f-electrons) and intersite (superexchange) correlation J (between localized f-electrons) on a triangular lattice for two different fillings. Numerical diagonalization and Monte-Carlo techniques are used to determine the ground state magnetic properties. Transitions from antiferromagnetic to ferromagnetic and again to re-entrant antiferromagnetic phase is observed in a wide range of parameter space. The magnetic moments of d- and f-electrons are observed to depend strongly on the value off, J and also on the total number of d-electrons (N-d). (C) 2015 Elsevier Ltd. All rights reserved.