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 a system of carbon nanotubes

Orientation cues for high-flying nocturnal insect migrants: do turbulence-induced temperature and velocity fluctuations indicate the mean wind flow?

Migratory insects flying at high altitude at night often show a degree of common alignment, sometimes with quite small angular dispersions around the mean. The observed orientation directions are often close to the downwind direction and this would seemingly be adaptive in that large insects could add their self-propelled speed to the wind speed, thus maximising their displacement in a given time. There are increasing indications that high-altitude orientation may be maintained by some intrinsic property of the wind rather than by visual perception of relative ground movement. Therefore, we first examined whether migrating insects could deduce the mean wind direction from the turbulent fluctuations in temperature. Within the atmospheric boundary-layer, temperature records show characteristic ramp-cliff structures, and insects flying downwind would move through these ramps whilst those flying crosswind would not. However, analysis of vertical-looking radar data on the common orientations of nocturnally migrating insects in the UK produced no evidence that the migrants actually use temperature ramps as orientation cues. This suggests that insects rely on turbulent velocity and acceleration cues, and refocuses attention on how these can be detected, especially as small-scale turbulence is usually held to be directionally invariant (isotropic). In the second part of the paper we present a theoretical analysis and simulations showing that velocity fluctuations and accelerations felt by an insect are predicted to be anisotropic even when the small-scale turbulence (measured at a fixed point or along the trajectory of a fluid-particle) is isotropic. Our results thus provide further evidence that insects do indeed use turbulent velocity and acceleration cues as indicators of the mean wind direction

Intriguing Electrostatic Potential of CO: Negative Bond-ends and Positive Bond-cylindrical-surface

The strong electronegativity of O dictates that the ground state of singlet CO has positively charged C and negatively charged O, in agreement with ab initio charge analysis, but in disagreement with the dipole direction. Though this unusual phenomenon has been fairly studied, the study of electrostatic potential (EP) for noncovalent interactions of CO is essential for better understanding. Here we illustrate that both C and O atom-ends show negative EP (where the C end gives more negative EP), favoring positively charged species, whereas the cylindrical surface of the CO bond shows positive EP, favoring negatively charged ones. This is demonstrated from the interactions of CO with Na+, Cl-, H2O, CO and benzene. It can be explained by the quadrupole driven electrostatic nature of CO (like N2) with very weak dipole moment. The EP is properly described by the tripole model taking into account the electrostatic multipole moments, which has a large negative charge at a certain distance protruded from C, a large positive charge on C, and a small negative charge on O. We also discuss the EP of the first excited triplet COopen0

Controlling infectious disease outbreaks: Lessons from mathematical modelling.

Accepted versio

The role of rapid diagnostics in managing Ebola epidemics

Ebola emerged in West Africa around December 2013 and swept through Guinea, Sierra Leone and Liberia, giving rise to 27,748 confirmed, probable and suspected cases reported by 29 July 2015. Case diagnoses during the epidemic have relied on polymerase chain reaction-based tests. Owing to limited laboratory capacity and local transport infrastructure, the delays from sample collection to test results being available have often been 2 days or more. Point-of-care rapid diagnostic tests offer the potential to substantially reduce these delays. We review Ebola rapid diagnostic tests approved by the World Health Organization and those currently in development. Such rapid diagnostic tests could allow early triaging of patients, thereby reducing the potential for nosocomial transmission. In addition, despite the lower test accuracy, rapid diagnostic test-based diagnosis may be beneficial in some contexts because of the reduced time spent by uninfected individuals in health-care settings where they may be at increased risk of infection; this also frees up hospital beds. We use mathematical modelling to explore the potential benefits of diagnostic testing strategies involving rapid diagnostic tests alone and in combination with polymerase chain reaction testing. Our analysis indicates that the use of rapid diagnostic tests with sensitivity and specificity comparable with those currently under development always enhances control, whether evaluated at a health-care-unit or population level. If such tests had been available throughout the recent epidemic, we estimate, for Sierra Leone, that their use in combination with confirmatory polymerase chain-reaction testing might have reduced the scale of the epidemic by over a third

External validation, update and development of prediction models for pre-eclampsia using an Individual Participant Data (IPD) meta-analysis: the International Prediction of Pregnancy Complication Network (IPPIC pre-eclampsia) protocol.

Background: Pre-eclampsia, a condition with raised blood pressure and proteinuria is associated with an increased risk of maternal and offspring mortality and morbidity. Early identification of mothers at risk is needed to target management. Methods/design: We aim to systematically review the existing literature to identify prediction models for pre-eclampsia. We have established the International Prediction of Pregnancy Complication Network (IPPIC), made up of 72 researchers from 21 countries who have carried out relevant primary studies or have access to existing registry databases, and collectively possess data from more than two million patients. We will use the individual participant data (IPD) from these studies to externally validate these existing prediction models and summarise model performance across studies using random-effects meta-analysis for any, late (after 34 weeks) and early (before 34 weeks) onset pre-eclampsia. If none of the models perform well, we will recalibrate (update), or develop and validate new prediction models using the IPD. We will assess the differential accuracy of the models in various settings and subgroups according to the risk status. We will also validate or develop prediction models based on clinical characteristics only; clinical and biochemical markers; clinical and ultrasound parameters; and clinical, biochemical and ultrasound tests. Discussion: Numerous systematic reviews with aggregate data meta-analysis have evaluated various risk factors separately or in combination for predicting pre-eclampsia, but these are affected by many limitations. Our large-scale collaborative IPD approach encourages consensus towards well developed, and validated prognostic models, rather than a number of competing non-validated ones. The large sample size from our IPD will also allow development and validation of multivariable prediction model for the relatively rare outcome of early onset pre-eclampsia. Trial registration: The project was registered on Prospero on the 27 November 2015 with ID: CRD42015029349

Improved Measurement of the Pseudoscalar Decay Constant fDsf_{D_{s}}

We present a new determination of the Ds decay constant, f_{Ds} using 5 million continuum charm events obtained with the CLEO II detector. Our value is derived from our new measured ratio of widths for Ds -> mu nu/Ds -> phi pi of 0.173+/- 0.021 +/- 0.031. Taking the branching ratio for Ds -> phi pi as (3.6 +/- 0.9)% from the PDG, we extract f_{Ds} = (280 +/- 17 +/- 25 +/- 34){MeV}. We compare this result with various model calculations.Comment: 23 page postscript file, postscript file also available through http://w4.lns.cornell.edu/public/CLN

Search for the Decays B^0 -> D^{(*)+} D^{(*)-}

Using the CLEO-II data set we have searched for the Cabibbo-suppressed decays B^0 -> D^{(*)+} D^{(*)-}. For the decay B^0 -> D^{*+} D^{*-}, we observe one candidate signal event, with an expected background of 0.022 +/- 0.011 events. This yield corresponds to a branching fraction of Br(B^0 -> D^{*+} D^{*-}) = (5.3^{+7.1}_{-3.7}(stat) +/- 1.0(syst)) x 10^{-4} and an upper limit of Br(B^0 -> D^{*+} D^{*-}) D^{*\pm} D^\mp and B^0 -> D^+ D^-, no significant excess of signal above the expected background level is seen, and we calculate the 90% CL upper limits on the branching fractions to be Br(B^0 -> D^{*\pm} D^\mp) D^+ D^-) < 1.2 x 10^{-3}.Comment: 12 page postscript file also available through http://w4.lns.cornell.edu/public/CLNS, submitted to Physical Review Letter

Anisotropic Charge Distribution and Anisotropic van der Waals Radius Leading to Intriguing Anisotropic Noncovalent Interactions

Although group (IV-VII) nonmetallic elements do not favor interacting with anionic species, there are counterexamples including the halogen bond. Such binding is known to be related to the charge deficiency because of the adjacent atom&apos;s electron withdrawing effect, which creates s/p-holes at the bond-ends. However, a completely opposite behavior is exhibited by N-2 and O-2, which have electrostatically positive/negative character around cylindrical-bond-surface/bond-ends. Inspired by this, here we elucidate the unusual features and origin of the anisotropic noncovalent interactions in the ground and excited states of the 2nd and 3rd row elements belonging to groups IV-VII. The anisotropy in charge distributions and van der Waals radii of atoms in such molecular systems are scrutinized. This provides an understanding of their unusual molecular configuration, binding and recognition modes involved in new types of molecular assembling and engineering. This work would lead to the design of intriguing molecular systems exploiting anisotropic noncovalent interactions.open