Engineering topics
Judson, Richard S.
Publications and source records attributed to Judson, Richard S..
Evaluation of a rapid, generic human gestational dose model
Chemical risk assessment considers potentially susceptible populations including pregnant women and developing fetuses. Humans encounter thousands of chemicals in their environments, few of which have been fully characterized. Toxicokinetic (TK) information is needed to relate chemical exposure to potentially bioactive tissue concentrations. Observational data describing human gestational exposures are unavailable for most chemicals, but physiologically based TK (PBTK) models estimate such exposures. Development of chemical-specific PBTK models requires considerable time and resources. As an alternative, generic PBTK approaches describe a standardized physiology and characterize chemicals with a set of standard physical and TK descriptors – primarily plasma protein binding and hepatic clearance. Here we report and evaluate a generic PBTK model of a human mother and developing fetus. We used a published set of formulas describing the major anatomical and physiological changes that occur during pregnancy to augment the High-Throughput Toxicokinetics (httk) software package. We simulated the ratio of concentrations in maternal and fetal plasma and compared to literature in vivo measurements. We evaluated the model with literature in vivo time-course measurements of maternal plasma concentrations in pregnant and non-pregnant women. Finally, we prioritized chemicals measured in maternal serum based on predicted fetal brain concentrations. This new model can be used for TK simulations of 859 chemicals with existing human-specific in vitro TK data as well as any new chemicals for which such data become available. Furthermore, this gestational model may allow for in vitro to in vivo extrapolation of point of departure doses relevant to reproductive and developmental toxicity.
State-to-state rates for the D + H2(v = 1, j = 1) yield HD(v-prime, j-prime) + H reaction - Predictions and measurements
A fully quantal wavepacket approach to reactive scattering in which the best available H3 potential energy surface was used enabled a comparison with experimentally determined rates for the D + H2(v = 1, j = 1) yield HD(v-prime = 0, 1, 2; j-prime) + H reaction at significantly higher total energies (1.4 to 2.25 electron volts) than previously possible. The theoretical results are obtained over a sufficient range of conditions that a detailed simulation of the experiment was possible, thus making this a definitive comparison of experiment and theory. Good to excellent agreement is found for the vibrational branching ratios and for the rotational distributions within each product vibrational level. However, the calculated rotational distributions are slightly hotter than the experimentally measured ones.
A comparison of three time-dependent wave packet methods for calculating electron-atom elastic scattering cross sections
Three time-dependent wave packet methods for performing elastic scattering calculations from screened Coulomb potentials are compared. The three methods are the time-dependent amplitude density method (TDADM), a Cayley-transform method (CTM), and the Chebyshev propagation method of Tal-Ezer and Kosloff. Both the TDADM and the CTM are based on a time-dependent integral equation for the wave function. In the first, the time-dependent amplitude density is propagated, while in the other two, the wave function is propagated. As a numerical example, phase shifts and cross sections are calculated using a screened Coulomb, Yukawa type potential over the range 200-1000 eV. It is found that, in most cases, all three methods yield comparable accuracy and are about equally efficient computationally. However for l = 0, where the Coulomb well is not screened by the centrifugal potential, the TDADM requires smaller grid spacings to maintain accuracy.
Total integral reactive cross sections for F + H2 yielding HF + H - Comparison of converged quantum, quasiclassical trajectory and experimental results
The paper reports converged quantum total integral reactive cross sections for the reaction F + H2 yielding HF + H, for initial rotational states j sub i = 0 and 1, using a time-dependent method. The results are compared to classical results and to the experimental results of Neumark et al. (1985). Strong quantum effects are found in the threshold region for both initial states (i.e., in the dependence of the reaction on initial state for low energies). The classical results agree better with experiment than do the quantum results; this appears to be due to errors in the potential used.
The application of time-dependent wavepacket methods to reactive scattering
Several methods for performing numerically exact reactive scattering calculations using time-dependent wavepackets are reviewed. The basic idea is to take the multiarrangement reactive problem and reformulate it as one or more inelastic ones. In the simplest method, total reaction probabilities are extracted by calculating the flux of the wavepacket as it leaves the interaction region in the direction of the reactive arrangement. To make this practical, complex potentials that absorb the wavepacket before it reaches the numerical grid boundary are used. Methods that generate observables ranging from total, energy-averaged reaction probabilities up to energy- and state-resolved S-matrix elements are used. Techniques for efficiently performing the necessary inelastic wavepacket propagation are also reviewed.
Time-dependent treatment of scattering. II - Novel integral equation approach to quantum wave packets
The novel wave-packet propagation scheme presented is based on the time-dependent form of the Lippman-Schwinger integral equation and does not require extensive matrix inversions, thereby facilitating application to systems in which some degrees of freedom express the potential in a basis expansion. The matrix to be inverted is a function of the kinetic energy operator, and is accordingly diagonal in a Bessel function basis set. Transition amplitudes for various orbital angular momentum quantum numbers are obtainable via either Fourier transform of the amplitude density from the time to the energy domain, or the direct analysis of the scattered wave packet.
A time-dependent wave packet approach to atom-diatom reactive collision probabilities - Theory and application to the H + H2(J = 0) system
This paper describes a new approach to the study of atom-diatom reactive collisions in three dimensions employing wave packets and the time-dependent Schroedinger equation. The method uses a projection operator approach to couple the inelastic and reactive portions of the total wave function and optical potentials to circumvent the necessity of using product arrangement coordinates. Reactive transition probabilities are calculated from the state resolved flux of the wave packet as it leaves the interaction region in the direction of the reactive arrangement channel. The present approach is used to obtain such vibrationally resolved probabilities for the three-dimensional H + H2 (J = 0) hydrogen exchange reaction, using a body-fixed system of coordinates.
Time-dependent wave-packet method for the complete determination of S-matrix elements for reactive molecular collisions in three dimensions
An alternative time-dependent wave-packet method for treating three-dimensional gas phase reactive atom-diatom collisions is presented. The method employs a nonreactive body-frame wave packet propagation procedure, made possible by judicious use of absorbing optical potentials, a novel scheme for interpolating the wave function from coordinates in one arrangement to those in another and the fact that the time-dependent Schroedinger equation is an initial-value problem. The last feature makes possible a computationally viable and accurate procedure for changing from one arrangement's coordinates to another. In addition, the method allows the determination of S-matrix elements over a wide range of energies from a single wave-packet propagation. The method is illustrated by carrying out detailed calculations of inelastic and reactive scattering in the H + H2 system using the Liu-Siegbahn-Truhlar-Horowitz potential surface.
Time-dependent treatment of scattering - Integral equation approaches using the time-dependent amplitude density
The time-dependent form of the Lippmann-Schwinger integral equation is used as the basis of several new wave packet propagation schemes. These can be formulated in terms of either the time-dependent wave function or a time-dependent amplitude density. The latter is nonzero only in the region of configuratiaon space for which the potential is nonzero, thereby in principle obviating the necessity of large grids or the use of complex absorbing potentials when resonances cause long collision times (leading, consequently, to long propagation times). Transition amplitudes are obtained in terms of Fourier transforms of the amplitude density from the time to the energy domain. The approach is illustrated by an application to a standard potential scattering model problem where, as in previous studies, the action of the kinetic energy operator is evaluated by fast Fourier transform (FFT) techniques.
Time dependent three-dimensional body frame quantal wave packet treatment of the H + H2 exchange reaction on the Liu-Siegbahn-Truhlar-Horowitz (LSTH) surface
The first successful application of the three-dimensional quantum body frame wave packet approach to reactive scattering is reported for the H + H2 exchange reaction on the LSTH potential surface. The method used is based on a procedure for calculating total reaction probabilities from wave packets. It is found that converged, vibrationally resolved reactive probabilities can be calculated with a grid that is not much larger than required for the pure inelastic calculation. Tabular results are presented for several energies.