Large scale parallelization in stochastic coupled cluster.

Coupled cluster theory is a vital cornerstone of electronic structure theory and is being applied to ever-larger systems. Stochastic approaches to quantum chemistry have grown in importance and offer compelling advantages over traditional deterministic algorithms in terms of computational demands, theoretical flexibility, or lower scaling with system size. We present a highly parallelizable algorithm of the coupled cluster Monte Carlo method involving sampling of clusters of excitors over multiple time steps. The behavior of the algorithm is investigated on the uniform electron gas and the water dimer at coupled-cluster levels including up to quadruple excitations. We also describe two improvements to the original sampling algorithm, full non-composite, and multi-spawn sampling. A stochastic approach to coupled cluster results in an efficient and scalable implementation at arbitrary truncation levels in the coupled cluster expansion.

Accurate Exchange-Correlation Energies for the Warm Dense Electron Gas.

The density matrix quantum Monte Carlo (DMQMC) method is used to sample exact-on-average N-body density matrices for uniform electron gas systems of up to 10^{124} matrix elements via a stochastic solution of the Bloch equation. The results of these calculations resolve a current debate over the accuracy of the data used to parametrize finite-temperature density functionals. Exchange-correlation energies calculated using the real-space restricted path-integral formalism and the k-space configuration path-integral formalism disagree by up to ∼10% at certain reduced temperatures T/T_{F}≤0.5 and densities r_{s}≤1. Our calculations confirm the accuracy of the configuration path-integral Monte Carlo results available at high density and bridge the gap to lower densities, providing trustworthy data in the regime typical of planetary interiors and solids subject to laser irradiation. We demonstrate that the DMQMC method can calculate free energies directly and present exact free energies for T/T_{F}≥1 and r_{s}≤2.

Understanding and improving the efficiency of full configuration interaction quantum Monte Carlo.

Within full configuration interaction quantum Monte Carlo, we investigate how the statistical error behaves as a function of the parameters which control the stochastic sampling. We define the inefficiency as a measure of the statistical error per particle sampling the space and per time step and show there is a sizeable parameter regime where this is minimised. We find that this inefficiency increases sublinearly with Hilbert space size and can be reduced by localising the canonical Hartree-Fock molecular orbitals, suggesting that the choice of basis impacts the method beyond that of the sign problem.

Linked coupled cluster Monte Carlo.

We consider a new formulation of the stochastic coupled cluster method in terms of the similarity transformed Hamiltonian. We show that improvement in the granularity with which the wavefunction is represented results in a reduction in the critical population required to correctly sample the wavefunction for a range of systems and excitation levels and hence leads to a substantial reduction in the computational cost. This development has the potential to substantially extend the range of the method, enabling it to be used to treat larger systems with excitation levels not easily accessible with conventional deterministic methods.

Interaction picture density matrix quantum Monte Carlo.

The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the N-body thermal density matrix and hence provides access to exact properties of many-particle quantum systems at arbitrary temperatures. We demonstrate that moving to the interaction picture provides substantial benefits when applying DMQMC to interacting fermions. In this first study, we focus on a system of much recent interest: the uniform electron gas in the warm dense regime. The basis set incompleteness error at finite temperature is investigated and extrapolated via a simple Monte Carlo sampling procedure. Finally, we provide benchmark calculations for a four-electron system, comparing our results to previous work where possible.

Semi-stochastic full configuration interaction quantum Monte Carlo: Developments and application.

We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable.

Minimising biases in full configuration interaction quantum Monte Carlo.

We show that Full Configuration Interaction Quantum Monte Carlo (FCIQMC) is a Markov chain in its present form. We construct the Markov matrix of FCIQMC for a two determinant system and hence compute the stationary distribution. These solutions are used to quantify the dependence of the population dynamics on the parameters defining the Markov chain. Despite the simplicity of a system with only two determinants, it still reveals a population control bias inherent to the FCIQMC algorithm. We investigate the effect of simulation parameters on the population control bias for the neon atom and suggest simulation setups to, in general, minimise the bias. We show a reweight ing scheme to remove the bias caused by population control commonly used in diffusion Monte Carlo [Umrigar et al., J. Chem. Phys. 99, 2865 (1993)] is effective and recommend its use as a post processing step.

Linear-scaling time-dependent density-functional theory in the linear response formalism.

We present an implementation of time-dependent density-functional theory (TDDFT) in the linear response formalism enabling the calculation of low energy optical absorption spectra for large molecules and nanostructures. The method avoids any explicit reference to canonical representations of either occupied or virtual Kohn-Sham states and thus achieves linear-scaling computational effort with system size. In contrast to conventional localised orbital formulations, where a single set of localised functions is used to span the occupied and unoccupied state manifold, we make use of two sets of in situ optimised localised orbitals, one for the occupied and one for the unoccupied space. This double representation approach avoids known problems of spanning the space of unoccupied Kohn-Sham states with a minimal set of localised orbitals optimised for the occupied space, while the in situ optimisation procedure allows for efficient calculations with a minimal number of functions. The method is applied to a number of medium sized organic molecules and a good agreement with traditional TDDFT methods is observed. Furthermore, linear scaling of computational cost with system size is demonstrated on (10,0) carbon nanotubes of different lengths.

The effect of quantization on the full configuration interaction quantum Monte Carlo sign problem.

The sign problem in full configuration interaction quantum Monte Carlo (FCIQMC) without annihilation can be understood as an instability of the psi-particle population to the ground state of the matrix obtained by making all off-diagonal elements of the Hamiltonian negative. Such a matrix, and hence the sign problem, is basis dependent. In this paper, we discuss the properties of a physically important basis choice: first versus second quantization. For a given choice of single-particle orbitals, we identify the conditions under which the fermion sign problem in the second quantized basis of antisymmetric Slater determinants is identical to the sign problem in the first quantized basis of unsymmetrized Hartree products. We also show that, when the two differ, the fermion sign problem is always less severe in the second quantized basis. This supports the idea that FCIQMC, even in the absence of annihilation, improves the sign problem relative to first quantized methods. Finally, we point out some theoretically interesting classes of Hamiltonians where first and second quantized sign problems differ, and others where they do not.

The sign problem and population dynamics in the full configuration interaction quantum Monte Carlo method.

The recently proposed full configuration interaction quantum Monte Carlo method allows access to essentially exact ground-state energies of systems of interacting fermions substantially larger than previously tractable without knowledge of the nodal structure of the ground-state wave function. We investigate the nature of the sign problem in this method and how its severity depends on the system studied. We explain how cancellation of the positive and negative particles sampling the wave function ensures convergence to a stochastic representation of the many-fermion ground state and accounts for the characteristic population dynamics observed in simulations.

Expression and characterization of recombinant interferon gamma (IFN-gamma) from the nine-banded armadillo (Dasypus novemcinctus) and its effect on Mycobacterium leprae-infected macrophages.

Armadillos (Dasypus novemcinctus) manifest the full histopathological spectrum of leprosy, and are hosts of choice for in vivo propagation of Mycobacterium leprae. Though potentially useful as a model of leprosy pathogenesis, few armadillo-specific reagents exist. We have identified a region of high homology to the interferon gamma (IFN-gamma) of other mammals within the recently published armadillo whole genomic sequence. cDNA was made from ConA-stimulated armadillo peripheral blood mononuclear cells (PBMC), amplified, and cloned into a pET expression vector for transformation and over-expression in Escherichia coli. The recombinant protein (rDnIFN-gamma) was characterized by western blot and its biological function confirmed with bioassays including intracellular killing of Toxoplasma gondii and induction of indoleamine 2, 3-dioxygenase activity. In using rIFN-gamma to activate macrophages from mice, humans or armadillos, similar to humans, rIFN-gamma-activated armadillo MPhi did not produce nitrite and or inhibit the viability of M. leprae in vitro. Conversely, murine rIFN-gamma-activated mouse MPhi produced high levels of nitrite and killed intracellular M. leprae in vitro. These data indicate that the response of armadillo MPhi to rDnIFN-gamma is similar to that which occurs in humans, and demonstrates a potentially important value of the armadillo as a model in leprosy research.

Immunological characterization of antigens encoded by the RD1 region of the Mycobacterium tuberculosis genome.

Development of immunoassays specific for the diagnosis of tuberculosis requires antigens unique to Mycobacterium tuberculosis. In a search for such antigens we tested six proteins encoded by RD1, a region present in M. tuberculosis and virulent M. bovis genomes but missing from the DNA of all substrains of M. bovis Bacillus Calmette-Guerin (BCG). The six proteins (Rv3871, Rv3872, Rv3873, MTSA-10, ESAT-6 and Rv3878) were purified to near-homogeneity from recombinant Escherichia coli. When tested for the ability to elicit antibody responses and delayed type hypersensitivity in tuberculous guinea pigs, only two of six antigens, ESAT-6 and MTSA-10, elicited strong skin reactions, while vigorous antibody responses were observed to all six proteins. When antibody responses to RD1 antigens were evaluated in sera from patients having pulmonary tuberculosis and from control subjects (patients having mycobacterioses other than tuberculosis, and healthy persons), a sizeable proportion (25%) of tuberculosis patients but none of the control subjects, had antibodies against MTSA-10 and/or ESAT-6. We conclude that MTSA-10 and ESAT-6 are promising candidates for immunodiagnostic assays specific for tuberculosis.

The application of proteomics in defining the T cell antigens of Mycobacterium tuberculosis.

The complete sequencing of the Mycobacterium tuberculosis genome offers a unique opportunity to fully elucidate the biology of this human pathogen. One aspect of significant importance is the definition of T cell antigens. This report describes the development and implementation of a proteomic approach to defining such antigens. Large quantities of subcellular protein fractions of M. tuberculosis were resolved by two-dimensional liquid phase electrophoresis (2-D LPE), resulting in 355 and 299 fractions of culture filtrate and cytosolic proteins, respectively. Analysis of these fractions against splenocytes of C57Bl/6 mice infected with M. tuberculosis resulted in the identification of 37 fractions that stimulated a dominant T cell response, as measured by the production of interferon-gamma. Additionally, when the 2-D LPE fractions were assayed against splenocytes harvested at 10 and 40 days post infection significant changes in the T cell response were observed. Molecular characterization of the proteins contained in each of the 38 immunodominant fractions by liquid chromatography-mass spectrometry and liquid chromatography-tandem mass spectrometry resulted in the identification of 30 individual proteins. Many of these represented previously defined antigens. However 17 of these proteins were novel T cell antigens. The data presented demonstrate that proteomics offers a rapid and facile approach for elucidation of immunodominant T cell antigens of pathogenic bacteria.

Characterization of the secreted MPT53 antigen of Mycobacterium tuberculosis.

MPT53 is a secreted protein of Mycobacterium tuberculosis. Southern transfer and hybridization showed mpt53 to be conserved in the M. tuberculosis complex and to have homology with DNA from Mycobacterium avium and other nontuberculous mycobacteria. However, anti-MPT53 polyclonal antibodies detected no antigen in the culture filtrates of M. avium and other nontuberculous mycobacteria. MPT53 of M. tuberculosis induced strong, tuberculosis-specific antibody responses in guinea pigs but induced no delayed-type hypersensitivity. Involvement in immune responses during human tuberculosis was very modest.

Further biochemical characterization of Mycobacterium leprae laminin-binding proteins.

It has been demonstrated that the alpha2 chain of laminin-2 present on the surface of Schwann cells is involved in the process of attachment of Mycobacterium leprae to these cells. Searching for M. leprae laminin-binding molecules, in a previous study we isolated and characterized the cationic proteins histone-like protein (Hlp) and ribosomal proteins S4 and S5 as potential adhesins involved in M. leprae-Schwann cell interaction. Hlp was shown to bind alpha2-laminins and to greatly enhance the attachment of mycobacteria to ST88-14 Schwann cells. In the present study, we investigated the laminin-binding capacity of the ribosomal proteins S4 and S5. The genes coding for these proteins were PCR amplified and their recombinant products were shown to bind alpha2-laminins in overlay assays. However, when tested in ELISA-based assays and in adhesion assays with ST88-14 cells, in contrast to Hlp, S4 and S5 failed to bind laminin and act as adhesins. The laminin-binding property and adhesin capacity of two basic host-derived proteins were also tested, and only histones, but not cytochrome c, were able to increase bacterial attachment to ST88-14 cells. Our data suggest that the alanine/lysine-rich sequences shared by Hlp and eukaryotic H1 histones might be involved in the binding of these cationic proteins to laminin.

Further biochemical characterization of Mycobacterium leprae laminin-binding proteins

It has been demonstrated that the alpha2 chain of laminin-2 present on the surface of Schwann cells is involved in the process of attachment of Mycobacterium leprae to these cells. Searching for M. leprae laminin-binding molecules, in a previous study we isolated and characterized the cationic proteins histone-like protein (Hlp) and ribosomal proteins S4 and S5 as potential adhesins involved in M. leprae-Schwann cell interaction. Hlp was shown to bind alpha2-laminins and to greatly enhance the attachment of mycobacteria to ST88-14 Schwann cells. In the present study, we investigated the laminin-binding capacity of the ribosomal proteins S4 and S5. The genes coding for these proteins were PCR amplified and their recombinant products were shown to bind alpha2-laminins in overlay assays. However, when tested in ELISA-based assays and in adhesion assays with ST88-14 cells, in contrast to Hlp, S4 and S5 failed to bind laminin and act as adhesins. The laminin-binding property and adhesin capacity of two basic host-derived proteins were also tested, and only histones, but not cytochrome c, were able to increase bacterial attachment to ST88-14 cells. Our data suggest that the alanine/lysine-rich sequences shared by Hlp and eukaryotic H1 histones might be involved in the binding of these cationic proteins to laminin

Cell-mediated immune response to tuberculosis antigens: comparison of skin testing and measurement of in vitro gamma interferon production in whole-blood culture.

Although delayed-type hypersensitivity skin testing with tuberculin purified protein derivative (PPD) is the standard for tuberculosis screening, its variability suggests the need for a more sensitive, noninvasive test. An in vitro whole-blood assay has been proposed as an alternative. Using health care worker volunteers, we confirmed the correlation between PPD skin test (PPD-ST) results (positive, induration of >15 mm) and a standardized gamma interferon (IFN-gamma) assay, QuantiFERON-TB (Q-IFN), manufactured by CSL Biosciences in Australia, and we evaluated Mycobacterium tuberculosis culture subfractions as potential substitutes for PPD. Twenty healthy volunteers with positive PPD-ST results and 20 PPD-ST-negative controls were enrolled. Whole blood was cultured with human PPD antigens (HuPPD), Mycobacterium avium complex (MAC) PPD, phytohemagglutinin (PHA), and four M. tuberculosis culture subfractions: low-molecular-weight culture, filtrate, culture filtrate without lipoarabinomannan, soluble cell wall proteins, and cytosolic proteins, all developed from M. tuberculosis strain H(37)RV. Secretion of IFN-gamma (expressed as international units per milliliter) was measured by an enzyme immunoassay. The PPD or subculture fraction response as a percentage of the PHA response was used to determine positivity. Sixteen of 20 PPD-ST-positive individuals were classified as M. tuberculosis positive by Q-IFN, and 1 was classified as MAC positive. Sixteen of 20 PPD-ST-negative individuals were M. tuberculosis negative by Q-IFN, 2 were MAC positive, and 2 were M. tuberculosis positive. The tuberculosis culture subfractions stimulated IFN-gamma production in PPD-ST-positive volunteers, and significant differences could be seen between the two PPD-ST groups with all subfractions except soluble cell wall protein; however, the response was variable and no better than the Q-IFN PPD. The agreement between the Q-IFN test and the PPD-ST was good (Cohen's kappa = 0.73). The Q-IFN assay can be a useful tool in further studies of immune responses to M. tuberculosis antigens.

Serodiagnostic potential of culture filtrate antigens of Mycobacterium tuberculosis.

Our studies of the humoral responses of tuberculosis (TB) patients have defined the repertoire of culture filtrate antigens of Mycobacterium tuberculosis that are recognized by antibodies from cavitary and noncavitary TB patients and demonstrated that the profile of antigens recognized changes with disease progression (K. Samanich et al., J. Infect. Dis. 178:1534-1538, 1998). We have identified several antigens with strong serodiagnostic potential. In the present study we have evaluated the reactivity of cohorts of human immunodeficiency virus (HIV)-negative, smear-positive; HIV-negative, smear-negative; and HIV-infected TB patients, with three of the candidate antigens, an 88-kDa protein, antigen (Ag) 85C, and MPT32, and compared the reactivity of the same patient cohort with the 38-kDa antigen and Ag 85A. We have also compared the reactivity of native Ag 85C and MPT32 with their recombinant counterparts. The evaluation of the reactivity was done by a modified enzyme-linked immunosorbent assay described earlier (S. Laal et al., Clin. Diag. Lab. Immunol. 4:49-56, 1997), in which all sera are preadsorbed against Escherichia coli lysates to reduce the levels of cross-reactive antibodies. Our results demonstrate that (i) antigens identified on the basis of their reactivity with TB patients' sera provide high sensitivities for serodiagnosis, (ii) recombinant Ag 85C and MPT32, expressed in E. coli, show reduced reactivity with human TB sera, and (iii) of the panel of antigens tested, the 88-kDa protein is the most promising candidate for serodiagnosis of TB in HIV-infected individuals. Moreover, these results reaffirm that both the extent of the disease and the bacterial load may play a role in determining the antigen profile recognized by antibodies.

MTSA-10, the product of the Rv3874 gene of Mycobacterium tuberculosis, elicits tuberculosis-specific, delayed-type hypersensitivity in guinea pigs.

In a search for new skin test reagents specific for tuberculosis, we found that the antigen encoded by gene Rv3874 of Mycobacterium tuberculosis elicited delayed-type hypersensitivity in M. tuberculosis-infected guinea pigs but not in control animals immunized with Mycobacterium bovis bacillus Calmette-Guérin (BCG) or Mycobacterium avium. The antigen, which was named MTSA-10 (for M. tuberculosis-specific antigen 10), is a prime candidate for a component of a new tuberculin that will allow discrimination by a skin test of latent M. tuberculosis infection from vaccination with BCG or from sensitization with environmental, nontuberculous mycobacteria.

A successful and necessary evolution to shared leadership: a hospital's story.

The success and adaptability of healthcare organizations will depend more and more on their ability to draw on the capabilities of their people. Tillsonburg District Memorial Hospital, a rural Ontario hospital, has evolved an organization and culture based on shared leadership and decision-making responsibility. Today this extends to front-line teams. This did not come about, however, without continuous effort. Successful transition takes preparation, guidance, much thought, commitment and patience.