Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “adaptive sampling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75

Auto-Gopher-II – An Autonomous Wireline Rotary-Hammer Ultrasonic Drill

Developing technologies that would enable future NASA exploration missions to penetrate deeper into the subsurface of planetary bodies for sample collection is of great importance. Performing these tasks while using minimal mass/volume systems and with low energy consumption is another set of requirements imposed on such technologies. A deep drill, called Auto-Gopher II, is currently being developed as a joint effort between JPL’s NDEAA laboratory and Honeybee Robotics Corp. The Auto-Gopher II is a wireline rotary-hammer drill that combines formation breaking by hammering using an ultrasonic actuator and cuttings removal by rotating a fluted auger bit. The hammering mechanism is based on the Ultrasonic/Sonic Drill/Corer (USDC) mechanism that has been developed as an adaptable tool for many drilling and coring applications. The USDC uses an intermediate free-flying mass to transform high frequency vibrations of a piezoelectric transducer horn tip into sonic hammering of the drill bit. The USDC concept was used in a previous task to develop an Ultrasonic/Sonic Ice Gopher and then integrated into a rotary hammer device to develop the Auto-Gopher-I. The lessons learned from these developments are being integrated into the development of the Auto-Gopher-II, an autonomous deep wireline drill with integrated cuttings and sample management and drive electronics. Subsystems of the wireline drill are being developed in parallel at JPL and Honeybee Robotics Ltd. This paper presents the development efforts of the piezoelectric actuator, cuttings removal and retention flutes and drive electronics.

Badescu, Mircea↗

Characterizing particle-based thermal storage performance using optical methods for use in next generation concentrating solar power plants

Concentrating Solar Power (CSP) generation is an attractive option for low-emission power generation; however, the high costs of thermal storage associated with concentrating solar create a large barrier for their use and adaptation into modern life. Lowering their operation costs, while maintaining high thermal storage and transfer performance is essential. Solid particle-based heat exchange systems can reduce CSP cost but are often less efficient. Efforts to increase their performance have led to use of binary size particle mixes. Presented is an optical-based thermal analysis technique used to measure near-wall thermal conductivity of particle beds essential in determining their heat exchanger efficiency. Modulated Photothermal Radiometry is used to make dynamic temperature measurements, allowing for the extraction of the most relevant thermal properties like thermal conductivity, specific heat, and effusivity. The system uses a modulated laser source causing a damped periodic heat flux, resulting in a frequency and thermal property dependent surface temperature, of which is measured using radiometry. Lock-In techniques are used to extrapolate the amplitude of the signal. Plotting the amplitude against the root angular frequency allows for effusivity measurement by ratio to a known sample. Using specific heat measurements from literature and density measurements, the thermal conductivity of the particle mixes can be calculated. The simplicity of MPTR to probe through the depth of the bed is ideal for use in CSP for dynamic thermal performance monitoring.

Corona, Javier↗

Lewis Acid Site Assisted Bifunctional Activity of Tin Doped Gallium Oxide and Its Application in Rechargeable Zn‐Air Batteries

The enhanced safety, superior energy, and power density of rechargeable metal‐air batteries make them ideal energy storage systems for application in energy grids and electric vehicles. However, the absence of a cost‐effective and stable bifunctional catalyst that can replace expensive platinum (Pt)‐based catalyst to promote oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) at the air cathode hinders their broader adaptation. Here, in this work, it is demonstrated that Tin (Sn) doped β‐gallium oxide (β‐Ga 2 O 3 ) in the bulk form can efficiently catalyze ORR and OER and, hence, be applied as the cathode in Zn‐air batteries. The Sn‐doped β‐Ga 2 O 3 sample with 15% Sn (Sn x =0.15 ‐Ga 2 O 3 ) displayed exceptional catalytic activity for a bulk, non‐noble metal‐based catalyst. When used as a cathode, the excellent electrocatalytic bifunctional activity of Sn x =0.15 ‐Ga 2 O 3 leads to a prototype Zn‐air battery with a high‐power density of 138 mW cm −2 and improved cycling stability compared to devices with benchmark Pt‐based cathode. The combined experimental and theoretical exploration revealed that the Lewis acid sites in β‐Ga 2 O 3 aid in regulating the electron density distribution on the Sn‐doped sites, optimize the adsorption energies of reaction intermediates, and facilitate the formation of critical reaction intermediate (O*), leading to enhanced electrocatalytic activity.

Lewis acid sites↗

Development of a portable high-T c spherical neutron polarimetry device at the Oak Ridge National Laboratory

Spherical neutron polarimetry is a powerful polarized neutron scattering technique used to determine complex magnetic structures which are only partly accessible by other methods. This technique measures the full neutron polarization change upon scattering from a sample by fully decoupling the incoming and outgoing neutron polarization with a zero-field chamber placed at the sample position. Recent advancements and testing are presented for a new spherical neutron polarimetry device utilizing high-T c superconducting YBCO films, PHiTPAD, at the High Flux Isotope Reactor at Oak Ridge National Laboratory. Furthermore, we introduce a conceptual design that utilizes wavelength-independent adiabatic transitions to adapt spherical neutron polarimetry for use with pulsed neutron sources, thereby expanding its potential applications in neutron scattering research.

47 OTHER INSTRUMENTATION↗

Low-rank matrix decomposition and spatio-temporal sparse recovery for STAP radar

We develop space-time adaptive processing (STAP) methods by leveraging the advantages of sparse signal processing techniques in order to detect a slowly-moving target. We observe that the inherent sparse characteristics of a STAP problem can be formulated as the low-rankness of clutter covariance matrix when compared to the total adaptive degrees-of-freedom, and also as the sparse interference spectrum on the spatio-temporal domain. By exploiting these sparse properties, we propose two approaches for estimating the interference covariance matrix. In the first approach, we consider a constrained matrix rank minimization problem (RMP) to decompose the sample covariance matrix into a low-rank positive semidefinite and a diagonal matrix. The solution of RMP is obtained by applying the trace minimization technique and the singular value decomposition with matrix shrinkage operator. Our second approach deals with the atomic norm minimization problem to recover the clutter response-vector that has a sparse support on the spatio-temporal plane. We use convex relaxation based standard sparse-recovery techniques to find the solutions. With extensive numerical examples, we demonstrate the performances of proposed STAP approaches with respect to both the ideal and practical scenarios, involving Doppler-ambiguous clutter ridges, spatial and temporal decorrelation effects. As a result, the low-rank matrix decomposition based solution requires secondary measurements as many as twice the clutter rank to attain a near-ideal STAP performance; whereas the spatio-temporal sparsity based approach needs a considerably small number of secondary data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Shortgrass prairie spectral measurements

The spectral methods of vegetation analysis not only measure herbage biomass on a nondestructive basis but also can be adapted to aircraft and satellite devices to map the spatial distribution over an area in an efficient and economical fashion. This study reviews the ground-based in situ field spectrometry in the 0.350-0.800 micron region of the spectrum. A statistical analysis of in situ spectroreflectance data from sample plots of the shortgrass prairie shows that green biomass, chlorophyll concentration, and leaf water content are directly interrelated to that composite property of the plot which is called functioning green biomass. Spectrocorrelation data indicate the spectral regions of optimum sensitivity for a remote estimation of the green biomass, chlorophyll, and leaf water content. The near-infrared region of the spectrum shows a high positive spectrocorrelation to these three sample parameters, regardless of the amount of standing dead vegetation.

Tucker, C. J.↗

Adaptive identification for the dynamics of large space structures

A new approach for identifying the dynamics of large space structures using noisy measurements is presented. In this approach, a modal system description is avoided and, hence, the a priori specification of model order. Instead, the system's order is recursively determined on-line, along with its parameters, using recursive lattice filters which are widely used in adaptive signal processing. The order determination is carried out based on statistical hypothesis testing theory. The approach is illustrated for the problem of identifying the structural dynamic characteristics of a free-free beam. Results on the effectiveness and accuracy of the identification scheme with respect to noise levels, sampling rates, and data window sizes are presented

Sundararajan, N.↗

Highly ordered vertical structure of Synechococcus populations within the one-millimeter-thick photic zone of a hot spring cyanobacterial mat

A variety of contemporary techniques were used to investigate the vertical distribution of thermophilic unicellular cyanobacteria, Synechococcus spp., and their activity within the upper 1-mm-thick photic zone of the mat community found in an alkaline siliceous hot spring in Yellowstone National Park in Wyoming. Detailed measurements were made over a diel cycle at a 61 degrees C site. Net oxygenic photosynthesis measured with oxygen microelectrodes was highest within the uppermost 100- to 200-microm-thick layer until midmorning, but as the day progressed, the peak of net activity shifted to deeper layers, stabilizing at a depth of 300 microm from midday throughout the afternoon. Examination of vertical thin sections by bright-field and autofluorescence microscopy revealed the existence of different populations of Synechococcus which form discrete bands at different vertical positions. Denaturing gradient gel electrophoresis analysis of PCR-amplified 16S rRNA gene segments from horizontal cryosections obtained at 100-microm-thick vertical intervals also suggested vertical stratification of cyanobacterial, green sulfur bacterium-like, and green nonsulfur bacterium-like populations. There was no evidence of diel migration. However, image analysis of vertical thin sections revealed the presence of a narrow band of rod-shaped Synechococcus cells in which the cells assumed an upright position. These upright cells, located 400 to 800 microm below the surface, were observed only in mat samples obtained around noon. In mat samples obtained at other time points, the cells were randomly oriented throughout the mat. These combined observations reveal the existence of a highly ordered structure within the very thin photic zone of this hot spring microbial mat, consisting of morphologically similar Synechococcus populations that are likely to be differentially adapted, some co-occurring with green sulfur bacterium-like populations, and all overlying green nonsulfur bacterium-like populations.

Non-NASA Center↗

Surveys of ISS Returned Hardware for MMOD Impacts

Since February 2001, the Hypervelocity Impact Technology (HVIT) group at the Johnson Space Center in Houston has performed 26 post-flight inspections on space exposed hardware that have been returned from the International Space Station. Data on 1,024 observations of MMOD damage have been collected from these inspections. Survey documentation typically includes impact feature location and size measurements as well as microscopic photography (25-200x). Sampling of impacts sites for projectile residue was performed for the largest features. Results of Scanning Electron Microscopy (SEM) analysis to discern impactor source is included in the database. This paper will summarize the post-flight MMOD inspections, and focus on two inspections in particular: (1) Pressurized Mating Adapter-2 (PMA-2) cover returned in 2015 after 1.6 years exposure with 26 observed damages, and (2) Airlock shield panels returned in 2010 after 8.7 years exposure with 58 MMOD damages. Feature sizes from the observed data are compared to predictions using the Bumper risk assessment code.

Hyde, James↗

Variable Pressure - Scanning Electron Microscopy (VP-SEM)

Variable Pressure (or Environmental) Scanning Electron Microscopy (VP-SEM) combined with Energy Dispersive X-ray Spectroscopy (EDS) is one of the most powerful methods for characterizing the sub-micron topography and chemical composition of uncoated samples. Terrestrially, VP-SEM is extensively used to non-destructively study geologic and manufactured materials with high spatial resolution (tens of nanometers) and large depth-of-field. An SEM offers a geologist a first survey of microscopic mineral phases via secondary electron imaging (SEI), which provides a topographic look at a sample, as well as backscattered electron (BEI) imaging, which contrasts the phases present based on their geochemistry (atomic number). A VP-SEM utilizes a gas in the sample chamber as a charge dissipation and signal amplification method, allowing analysis of a sample without preparation in the form of a conductive coating. A miniaturized VP-SEM operating in-situ on a lander or rover would be able to use the CO2-rich Martian atmosphere (e.g., Nier et al., 1976; Williams 2016) as an imaging medium for this purpose. Adaptation of a VP-SEM for in-situ Mars surface studies will provide a new imaging capability (via SEI) that is at least an order of magnitude higher resolution than the Mars Hand Lens Imager (MAHLI) on the Mars Science Laboratory (Williams et al., 2015), and equal to or better than the achieved resolution of the Atomic Force Microscope on the Phoenix Mars lander (Pike et al., 2011). In addition, the MVP-SEM is capable of BSI and simultaneous chemical analysis of the imaged region. In this sense, the MVP-SEM can be regarded as an instrument suite that will provide a new set of information not achievable by any other instrument. The Miniaturized Variable Pressure Scanning Electron Microscope (MVP-SEM) was designed, built, and benchtop tested by a team at NASA’s Marshall Space Flight Center and Jet Propulsion Laboratory (JPL), Jacobs Space Exploration Group, Applied Physics Technologies, Inc. (AP-Tech), and Creare LLC, working with a team of technical and science collaborators. Benchtop testing was successful, proving concept feasibility. To date, the MVP-SEM has achieved an imaging resolution of <100 nm in the lab, and with continued optimization, even better performance (~50 nm resolution) is possible. Use in-situ on the lunar or other planetary surfaces would require some redesign to optimize instrument performance, but such an instrument would be equally useful.

Variable Pressure Scanning Electron Microscopy↗

Miniature Variable Pressure - Scanning Electron Microscopy (MVP-SEM)

Variable Pressure (or Environmental) Scanning Electron Microscopy (VP-SEM) combined with Energy Dispersive X-ray Spectroscopy (EDS) is one of the most powerful methods for characterizing the sub-micron topography and chemical composition of uncoated samples. Terrestrially, VP-SEM is extensively used to non-destructively study geologic and manufactured materials with high spatial resolution (tens of nanometers) and large depth-of-field. An SEM offers a geologist a first survey of microscopic mineral phases via secondary electron imaging (SEI), which provides a topographic look at a sample, as well as backscattered electron (BEI) imaging, which contrasts the phases present based on their geochemistry (atomic number). A VP-SEM utilizes a gas in the sample chamber as a charge dissipation and signal amplification method, allowing analysis of a sample without preparation in the form of a conductive coating. A miniaturized VP-SEM operating in-situ on a lander or rover would be able to use the CO2-rich Martian atmosphere (e.g., Nier et al., 1976; Williams 2016) as an imaging medium for this purpose. Adaptation of a VP-SEM for in-situ Mars surface studies will provide a new imaging capability (via SEI) that is at least an order of magnitude higher resolution than the Mars Hand Lens Imager (MAHLI) on the Mars Science Laboratory (Williams et al., 2015), and equal to or better than the achieved resolution of the Atomic Force Microscope on the Phoenix Mars lander (Pike et al., 2011). In addition, the MVP-SEM is capable of BSI and simultaneous chemical analysis of the imaged region. In this sense, the MVP-SEM can be regarded as an instrument suite that will provide a new set of information not achievable by any other instrument. The Miniaturized Variable Pressure Scanning Electron Microscope (MVP-SEM) was designed, built, and benchtop tested by a team at NASA’s Marshall Space Flight Center and Jet Propulsion Laboratory (JPL), Jacobs Space Exploration Group, Applied Physics Technologies, Inc. (AP-Tech), and Creare LLC, working with a team of technical and science collaborators. Benchtop testing was successful, proving concept feasibility. To date, the MVP-SEM has achieved an imaging resolution of <100 nm in the lab, and with continued optimization, even better performance (~50 nm resolution) is possible. Use in-situ on the lunar or other planetary surfaces would require some redesign to optimize instrument performance, but such an instrument would be equally useful.

Variable Pressure Scanning Electron Microscopy↗

Well Design, Operation and Maintenance and Proactive Rehabilitation to Optimize Remediation Well Performance - 18605

Optimizing performance of extraction, injection and monitoring wells, which are integral components of groundwater remediation systems, result in lower life cycle costs and a higher probability of meeting project objectives. These wells provide the means to extract contaminated groundwater, inject reagents, treated water and clean water to create groundwater reaction/flushing zones, and to provide groundwater level, contaminant concentration and geochemical data. Initial costs associated with well installation and associated pipelines and treatment systems are significant project expenditures related to site remediation. However, the life cycle costs for wells are dominated by the long-term operation and maintenance (O and M) expenditures. Remediation wells are particularly susceptible to various types of fouling that impact well performance and, therefore, drive up operating costs. Simply defined, fouling is the result of any biological, geochemical or physical process that negatively impacts the intended performance of a remediation well. Fouling decreases well performance by effectively reducing hydraulic communication between the well and the surrounding formation. Oftentimes this occurs via clogging screens and/or gravel pack/formation pore space by bacterial or chemical byproducts and gases, which can lead to corrosion of metal components and spreading fouling deposits to pipelines and surface treatment systems. Ultimately, this results in reduced hydraulic performance in extraction and injection wells and result in erroneous water level and chemical parameter concentration measurements in monitoring wells. Remediation well fouling cannot be eliminated, and natural fouling processes are exacerbated by the installation and operation of remedial systems. Steps required to extend the life of remediation wells and to lower life cycle costs are discussed in this paper in order to establish the required approach to effective system management. These optimization approaches have been developed through years of operation and maintenance of complex large-scale plume remediation systems utilizing a variety of subsurface in-situ and ex-situ technologies. Remediation well performance can be improved by: - careful site characterization, - appropriate well design, - proper well installation and development, - thoughtful and adaptable operation and maintenance, - well performance maintenance tracking and, - proactive well rehabilitation. Clear understanding of the site hydro-stratigraphy and geochemical settings places constraints on the compatibility of well materials and the configurations of well screens. Proper well design, including collection and analysis (sieve/grain-size) of screen zone formation samples, assures that wells can be optimally installed, developed and redeveloped during well rehabilitation. Design considerations should include proper sizing of filter packs and well screens, casing/screen diameter, well material selection/compatibility, expected life/use, and well head configuration to allow for future access. Applying appropriate and effective well installation and development methods are of paramount importance and rely on a good understanding of site characteristics (e.g., geologic setting) as well as local drilling company expertise and experience. Different approaches to operation and maintenance can greatly affect the degree of fouling, these include limiting reagent concentrations, maintaining submerged well screens and filtering injection water to prevent the introduction of fouling enhancing constituents. Tracking well performance by monitoring changes in system operation and metrics such as specific capacity or injectivity for extraction and injection wells, respectively, or turbidity and chemical analyses for monitoring wells is key to identifying well maintenance issues early when well rehabilitation can be most effective. Proactive well rehabilitation can include mechanical brushing, surging and jetting, chemical treatment or other specialized approaches; however, the approach needs to be based on a site-specific understanding of the predominant fouling mechanisms. Collection and analyses of fouling deposits can provide guidance for selection of rehabilitation methods and chemicals that will be most effective. Application of these approaches have led to extended operational life for reagent injection wells, which can foul quickly, from months to years. Additionally, proper maintenance will often lead to increased and extended productivity of the well. For extraction wells, this translates to a higher sustainable yield, which means shorter cleanup times. Similarly, for injection wells, proper maintenance allows for the introduction of a greater volume of reagents over a shorter period, thus, shorter cleanup times that reduce overall project costs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Locoregional Control of Non-Small Cell Lung Cancer in Relation to Automated Early Assessment of Tumor Regression on Cone Beam Computed Tomography

Purpose: Large interindividual variations in volume regression of non-small cell lung cancer (NSCLC) are observable on standard cone beam computed tomography (CBCT) during fractionated radiation therapy. Here, a method for automated assessment of tumor volume regression is presented and its potential use in response adapted personalized radiation therapy is evaluated empirically. Methods and Materials: Automated deformable registration with calculation of the Jacobian determinant was applied to serial CBCT scans in a series of 99 patients with NSCLC. Tumor volume at the end of treatment was estimated on the basis of the first one third and two thirds of the scans. The concordance between estimated and actual relative volume at the end of radiation therapy was quantified by Pearson's correlation coefficient. On the basis of the estimated relative volume, the patients were stratified into 2 groups having volume regressions below or above the population median value. Kaplan-Meier plots of locoregional disease-free rate and overall survival in the 2 groups were used to evaluate the predictive value of tumor regression during treatment. Cox proportional hazards model was used to adjust for other clinical characteristics. Results: Automatic measurement of the tumor regression from standard CBCT images was feasible. Pearson's correlation coefficient between manual and automatic measurement was 0.86 in a sample of 9 patients. Most patients experienced tumor volume regression, and this could be quantified early into the treatment course. Interestingly, patients with pronounced volume regression had worse locoregional tumor control and overall survival. This was significant on patient with non-adenocarcinoma histology. Conclusions: Evaluation of routinely acquired CBCT images during radiation therapy provides biological information on the specific tumor. This could potentially form the basis for personalized response adaptive therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Optimal Nested Test Plan for Combinatorial Quantitative Group Testing

We consider the quantitative group testing problem where the objective is to identify defective items in a given population based on results of tests performed on subsets of the population. Under the quantitative group testing model, the result of each test reveals the number of defective items in the tested group. The minimum number of tests achievable by nested test plans was established by Aigner and Schughart in 1985 within a minimax framework. Furthermore, the optimal nested test plan offering this performance, however, was not obtained. In this paper, we establish the optimal nested test plan in a closed form. This optimal nested test plan is also order optimal among all test plans as the population size approaches infinity. Using heavy-hitter detection as a case study, we show via simulation examples orders of magnitude improvement of the group testing approach over two prevailing sampling-based approaches in detection accuracy and counter consumption. Other applications include anomaly detection and wideband spectrum sensing in cognitive radio systems.

97 MATHEMATICS AND COMPUTING↗

The structure of CO{sub 2} hydrate between 0.7 and 1.0 GPa

A deuterated sample of CO{sub 2} structure I (sI) clathrate hydrate (CO{sub 2}·8.3 D{sub 2}O) has been formed and neutron diffraction experiments up to 1.0 GPa at 240 K were performed. The sI CO{sub 2} hydrate transformed at 0.7 GPa into the high pressure phase that had been observed previously by Hirai et al. [J. Phys. Chem. 133, 124511 (2010)] and Bollengier et al. [Geochim. Cosmochim. Acta 119, 322 (2013)], but which had not been structurally identified. The current neutron diffraction data were successfully fitted to a filled ice structure with CO{sub 2} molecules filling the water channels. This CO{sub 2}+water system has also been investigated using classical molecular dynamics and density functional ab initio methods to provide additional characterization of the high pressure structure. Both models indicate the water network adapts a MH-III “like” filled ice structure with considerable disorder of the orientations of the CO{sub 2} molecule. Furthermore, the disorder appears to be a direct result of the level of proton disorder in the water network. In contrast to the conclusions of Bollengier et al., our neutron diffraction data show that the filled ice phase can be recovered to ambient pressure (0.1 MPa) at 96 K, and recrystallization to sI hydrate occurs upon subsequent heating to 150 K, possibly by first forming low density amorphous ice. Unlike other clathrate hydrate systems, which transform from the sI or sII structure to the hexagonal structure (sH) then to the filled ice structure, CO{sub 2} hydrate transforms directly from the sI form to the filled ice structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Consequences of sludge composition on combustion performance derived from thermogravimetry analysis

Highlights: • Volatiles, particularly proteins, play a key role in sludge combustion. • Sludge combustion performance varies with different sludge organic concentrations. • Carbohydrates significantly affect the combustion rate in the second stage. • Combustion performance of digested sludge is more negative compared with others. - Abstract: Wastewater treatment plants produce millions of tons of sewage sludge. Sewage sludge is recognized as a promising feedstock for power generation via combustion and can be used for energy crisis adaption. We aimed to investigate the quantitative effects of various sludge characteristics on the overall sludge combustion process performance. Different types of sewage sludge were derived from numerous wastewater treatment plants in Beijing for further thermogravimetric analysis. Thermogravimetric–differential thermogravimetric curves were used to compare the performance of the studied samples. Proximate analytical data, organic compositions, elementary composition, and calorific value of the samples were determined. The relationship between combustion performance and sludge composition was also investigated. Results showed that the performance of sludge combustion was significantly affected by the concentration of protein, which is the main component of volatiles. Carbohydrates and lipids were not correlated with combustion performance, unlike protein. Overall, combustion performance varied with different sludge organic composition. The combustion rate of carbohydrates was higher than those of protein and lipid, and carbohydrate weight loss mainly occurred during the second stage (175–300 °C). Carbohydrates have a substantial effect on the rate of system combustion during the second stage considering the specific combustion feature. Additionally, the combustion performance of digested sewage sludge is more negative than the others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computing thermal Wigner densities with the phase integration method

We discuss how the Phase Integration Method (PIM), recently developed to compute symmetrized time correlation functions [M. Monteferrante, S. Bonella, and G. Ciccotti, Mol. Phys. 109, 3015 (2011)], can be adapted to sampling/generating the thermal Wigner density, a key ingredient, for example, in many approximate schemes for simulating quantum time dependent properties. PIM combines a path integral representation of the density with a cumulant expansion to represent the Wigner function in a form calculable via existing Monte Carlo algorithms for sampling noisy probability densities. The method is able to capture highly non-classical effects such as correlation among the momenta and coordinates parts of the density, or correlations among the momenta themselves. By using alternatives to cumulants, it can also indicate the presence of negative parts of the Wigner density. Both properties are demonstrated by comparing PIM results to those of reference quantum calculations on a set of model problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗