WavePacket: A Matlab package for numerical quantum dynamics. III: Quantum-classical simulations and surface hopping trajectories

WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I [Comp. Phys. Comm. 213, 223-234 (2017)] and Part II [Comp. Phys. Comm. 228, 229-244 (2018)] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagations to WavePacket. In those simulations classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces trajectories may switch between surfaces. To model these transitions, two classes of stochastic algorithms have been implemented: (1) J. C. Tully's fewest switches surface hopping and (2) Landau-Zener based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring non-adiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.0.2 which is essentially an object-oriented rewrite of previous versions, allowing to perform fully classical, quantum-classical and quantum-mechanical simulations on an equal footing, i.e., for the same physical system described by the same WavePacket input. The software package is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics are available

Leveraging legacy codes to distributed problem solving environments: A web service approach

This paper describes techniques used to leverage high performance legacy codes as CORBA components to a distributed problem solving environment. It first briefly introduces the software architecture adopted by the environment. Then it presents a CORBA oriented wrapper generator (COWG) which can be used to automatically wrap high performance legacy codes as CORBA components. Two legacy codes have been wrapped with COWG. One is an MPI-based molecular dynamic simulation (MDS) code, the other is a finite element based computational fluid dynamics (CFD) code for simulating incompressible Navier-Stokes flows. Performance comparisons between runs of the MDS CORBA component and the original MDS legacy code on a cluster of workstations and on a parallel computer are also presented. Wrapped as CORBA components, these legacy codes can be reused in a distributed computing environment. The first case shows that high performance can be maintained with the wrapped MDS component. The second case shows that a Web user can submit a task to the wrapped CFD component through a Web page without knowing the exact implementation of the component. In this way, a user’s desktop computing environment can be extended to a high performance computing environment using a cluster of workstations or a parallel computer

Imaging Molecules from Within: Ultra-fast, {\AA}ngstr\"om Scale Structure Determination of Molecules via Photoelectron Holography using Free Electron Lasers

A new scheme based on (i) upcoming brilliant X-ray Free Electron Laser (FEL) sources, (ii) novel energy and angular dispersive, large-area electron imagers and (iii) the well-known photoelectron holography is elaborated that provides time-dependent three-dimensional structure determination of small to medium sized molecules with {\AA}ngstr\"om spatial and femtosecond time resolution. Inducing molecular dynamics, wave-packet motion, dissociation, passage through conical intersections or isomerization by a pump pulse this motion is visualized by the X-ray FEL probe pulse launching keV photoelectrons within few femtoseconds from specific and well-defined sites, deep core levels of individual atoms, inside the molecule. On their way out the photoelectrons are diffracted generating a hologram on the detector that encodes the molecular structure at the instant of photoionization, thus providing 'femtosecond snapshot images of the molecule from within'. Detailed calculations in various approximations of increasing sophistication are presented and three-dimensional retrieval of the spatial structure of the molecule with {\AA}ngstr\"om spatial resolution is demonstrated. Due to the large photo-absorption cross sections the method extends X-ray diffraction based, time-dependent structure investigations envisioned at FELs to new classes of samples that are not accessible by any other method. Among them are dilute samples in the gas phase such as aligned, oriented or conformer selected molecules, ultra-cold ensembles and/or molecular or cluster objects containing mainly light atoms that do not scatter X-rays efficiently.Comment: 18 pages, 11 figure

Reactive Programming of Simulations in Physics

We consider the Reactive Programming (RP) approach to simulate physical systems. The choice of RP is motivated by the fact that RP genuinely offers logical parallelism, instantaneously broadcast events, and dynamic creation/destruction of parallel components and events. To illustrate our approach, we consider the implementation of a system of Molecular Dynamics, in the context of Java with the Java3D library for 3D visualisation

An Introduction to Programming for Bioscientists: A Python-based Primer

Computing has revolutionized the biological sciences over the past several decades, such that virtually all contemporary research in the biosciences utilizes computer programs. The computational advances have come on many fronts, spurred by fundamental developments in hardware, software, and algorithms. These advances have influenced, and even engendered, a phenomenal array of bioscience fields, including molecular evolution and bioinformatics; genome-, proteome-, transcriptome- and metabolome-wide experimental studies; structural genomics; and atomistic simulations of cellular-scale molecular assemblies as large as ribosomes and intact viruses. In short, much of post-genomic biology is increasingly becoming a form of computational biology. The ability to design and write computer programs is among the most indispensable skills that a modern researcher can cultivate. Python has become a popular programming language in the biosciences, largely because (i) its straightforward semantics and clean syntax make it a readily accessible first language; (ii) it is expressive and well-suited to object-oriented programming, as well as other modern paradigms; and (iii) the many available libraries and third-party toolkits extend the functionality of the core language into virtually every biological domain (sequence and structure analyses, phylogenomics, workflow management systems, etc.). This primer offers a basic introduction to coding, via Python, and it includes concrete examples and exercises to illustrate the language's usage and capabilities; the main text culminates with a final project in structural bioinformatics. A suite of Supplemental Chapters is also provided. Starting with basic concepts, such as that of a 'variable', the Chapters methodically advance the reader to the point of writing a graphical user interface to compute the Hamming distance between two DNA sequences.Comment: 65 pages total, including 45 pages text, 3 figures, 4 tables, numerous exercises, and 19 pages of Supporting Information; currently in press at PLOS Computational Biolog

Improvements to the APBS biomolecular solvation software suite

The Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve the equations of continuum electrostatics for large biomolecular assemblages that has provided impact in the study of a broad range of chemical, biological, and biomedical applications. APBS addresses three key technology challenges for understanding solvation and electrostatics in biomedical applications: accurate and efficient models for biomolecular solvation and electrostatics, robust and scalable software for applying those theories to biomolecular systems, and mechanisms for sharing and analyzing biomolecular electrostatics data in the scientific community. To address new research applications and advancing computational capabilities, we have continually updated APBS and its suite of accompanying software since its release in 2001. In this manuscript, we discuss the models and capabilities that have recently been implemented within the APBS software package including: a Poisson-Boltzmann analytical and a semi-analytical solver, an optimized boundary element solver, a geometry-based geometric flow solvation model, a graph theory based algorithm for determining p$K_a$ values, and an improved web-based visualization tool for viewing electrostatics