Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Interactions”

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,207 records · Page 67

Basic research in solar physics

This grant, dating back more than 20 years has supported a variety of investigations of the chromospheres and coronae of the Sun and related cool stars by the Principal Investigator, his postdocs and graduate students, and colleagues at other institutions. This work involved studies of radiative transfer and spectral line formation theory, and the application of these techniques to the analysis of spectra obtained from space and ground-based observatories in the optical, ultraviolet, x-ray and radio portions of the spectrum. Space observations have included the analysis of spectra from OSO-7, Skylab, SMM, and the HRTS rocket experiments. Recent work has concentrated on the interaction of magnetic fields, plasma and radiation in the outer atmospheres of the Sun and other magnetically active stars with different fundamental parameters. Our study of phenomena common to the Sun and stars, the 'solar-stellar connection', can elucidate the fundamental physics, because spatially-resolved observations of the Sun provide us with the 'groundtruth,' while interpretation of stellar data permit us to isolate those parameters critical to stellar activity. Recently, we have studied the differences in physical properties between solar regions of high magnetic flux density and the surrounding plasma. High-resolution CN and CO spectroheliograms have been used to model the thermal inhomogeneities driven by unstable CO cooling, and we have analyzed spatially resolved UV spectra from HRTS to model the thermal structure and energy balance of small-scale structures. The study of nonlinear relations between atmospheric radiative losses and the photospheric magnetic flux density has been continued. We have also proposed a new model for the decay of plages by random walk diffusion of magnetic flux. Our analysis of phenomena common to the Sun and stars included the application of available spectroscopic diagnostics, establishing evidence that the atmospheres of the least active stars are heated at a 'basal' rate that is also found in the centers of solar supergranules, and using the Doppler-imaging technique to measure the position, size, and brightness of stellar active regions. We are computing multi-component models for solar and stellar atmospheres, and models for coronal loops and for the transition-region down flows. The study of solar and stellar flares permits us to assess the role of turbulent energy transport, to pinpoint the mechanism behind Type I radio bursts, to determine whether plasma radiation or cyclotron maser is responsible for microwave flares on M dwarfs, and to extend our knowledge of the basic physics pertinent to cyclotron-maser processes operating on the Sun.

Linsky, Jeffrey L.↗

Soft-sphere simulations of a planar shock interaction with a granular bed

Here we consider the problem of shock propagation through a layer of spherical particles. A point particle force model is used to capture the shock-induced aerodynamic force acting upon the particles. The discrete element method (DEM) code liggghts is used to implement the shock-induced force as well as to capture the collisional forces within the system. A volume-fraction-dependent drag correction is applied using Voronoi tessellation to calculate the volume of fluid around each individual particle. A statistically stationary frame is chosen so that spatial and temporal averaging can be performed to calculate ensemble-averaged macroscopic quantities, such as the granular temperature. A parametric study is carried out by varying the coefficient of restitution for three sets of multiphase shock conditions. A self-similar profile is obtained for the granular temperature that is dependent on the coefficient of restitution. A traveling wave structure is observed in the particle concentration downstream of the shock and this instability arises from the volume-fraction-dependent drag force. The intensity of the traveling wave increases significantly as inelastic collisions are introduced. Downstream of the shock, the variance in Voronoi volume fraction is shown to have a strong dependence upon the coefficient of restitution, indicating clustering of particles induced by collisional dissipation. Statistics of the Voronoi volume are computed upstream and downstream of the shock and compared to theoretical results for randomly distributed hard spheres.

42 ENGINEERING↗

Dynamical Decoupling for Measuring and Suppressing Crosstalk

Dynamical decoupling (DD) is a noise-mitigating strategy in which sequences of pulses are applied to single qubits to average out their interaction with the environment. DD has been extensively studied and demonstrated for suppressing single-qubit decoherence and can be tailored for different noise spectrum. We report another important adaptation of DD where crosstalk between qubits are suppressed. We demonstrate the efficiency of this procedure on quantum circuits on superconducting transmon-based quantum devices. We designed a family of syncopated DD sequences that effectively suppress ZZ coupling between qubit pairs, which is the dominating crosstalk form on the device. We insert DD to a quantum circuit whenever single qubits are idle (often during two-qubits gates on other qubits). While standard periodic DD suppress crosstalk between these qubits and their neighbors, the syncopated DD further decouples crosstalk between these qubits. We further designed short sequences that maximize the application of DD without adding time to the quantum circuit execution. Such DD sequences yield significant improvement of the performance of the algorithm on the hardware. The performance is further boosted by combining DD with another mitigation strategy, randomized compilation. Our work demonstrated that syncopated DD is effective and practical way to suppress crosstalk in quantum circuits and serves as a great probe to characterize the crosstalk and inform hardware design.

Quantum Computing↗

Evaluating ecosystem water use efficiency under drought stress: a case study of the Helan Mountain region, northwest China

Context Water use efficiency (WUE) is a fundamental ecological indicator links carbon assimilation and water loss in terrestrial ecosystems. Understanding its responses to drought stress is essential for adaptive ecosystem management, particularly in climate-sensitive mountain landscapes. Objectives This study aimed to investigate drought-driven variations in WUE across major vegetation types in the Helan Mountain region of Northwest China. Specifically, we sought to identify dominant ecological drivers of WUE variability and to disentangle their relative importance and causal pathways. Methods We quantified WUE using the Moderate Resolution Imaging Spectroradiometer (MODIS) products and the Drought Severity Index (DSI) data from 2001 to 2020. To examine WUE – drought relationships across contrasting vegetation types, we employed a spatially explicit analytical framework integrating Random Forest (RF) modeling, partial correlation analysis, and structural equation modeling (SEM). Results Regional WUE exhibited relatively stable interannual dynamics, yet pronounced spatial heterogeneity that was strongly modulated by drought conditions. Vegetation properties, particularly Leaf Area Index (LAI) and Normalized Difference Vegetation Index (NDVI), emerged as the dominant determinants of WUE, with NDVI alone explaining over 20% of its spatial variance in forest and grassland during non-drought periods. SEM analyses revealed that climate forcing influenced WUE mainly through indirect pathways mediated by soil moisture availability and vegetation structural dynamics, rather than through direct climatic controls. Among all regulating factors, LAI acted as the central control node governing ecosystem carbon–water coupling. In contrast, short-term climatic stress, especially atmospheric demand and drought duration, exerted weak or negative direct effects on WUE. Ecosystem-specific responses were observed, with croplands mainly regulated by soil water availability, whereas forests and grasslands showed more sensitive to atmospheric drought stress. Together, these results reveal a hierarchical control framework where soil–vegetation interactions mediate climate impacts on WUE, driving strong spatial heterogeneity in drought responses across mountain landscapes. Conclusions Our findings highlight the pivotal role of indirect drought effects mediated by vegetation and soil processes in shaping ecosystem WUE. The identified soil–vegetation–climate regulatory hierarchy provides mechanistic insight into landscape–scale drought sensitivity and supports integrated modeling approaches for evaluating ecosystem resilience and sustainable management in arid mountain regions.

China↗

Néel domain walls with bistable chirality in a perpendicularly magnetized ferrimagnetic insulator

Field-free spin-orbit torque-driven domain wall motion in magnetic thin films with perpendicular magnetic anisotropy (PMA) requires the domain walls to have Néel character. Conventionally, Néel domain walls are stabilized by the Dzyaloshinskii-Moriya interaction (DMI) in ultrathin films. Here, in a europium iron garnet thin film with PMA and an additional uniaxial in-plane anisotropy, we demonstrate two bistable Néel domain wall states in the absence of DMI, and the capability to toggle the wall states with an in-plane field pulse and consequently their directions of motion under a current pulse. We present a phase diagram for the bistable Néel domain wall states as a function of in-plane field pulse width and amplitude. By fitting the experimental data to an analytical model of Néel wall reversal through the nucleation and propagation of Bloch lines, we extract the length of the initial reversed domain wall segment and Bloch line nucleation energy barrier. Current-driven motion of in-plane anisotropy stabilized Néel walls is qualitatively different from that of DMI-stabilized ones owing to the different symmetry of the effective fields that stabilize the Néel configuration. Furthermore, we present a proof of principle demonstration for 2-bit random number generation based on the stochastic reversal of domain wall chirality. These results provide critical insight into the topological energy barrier of Bloch lines and identify paths towards domain wall-based memory and computing devices.

36 MATERIALS SCIENCE↗

Unit-Sphere Multiaxial Stochastic-Strength Model Applied to Anisotropic and Composite Materials

Models that predict the failure probability of brittle materials under multiaxial loading have been developed by authors such as Batdorf, Evans, and Matsuo. These "unit-sphere" models assume that the strength-controlling flaws are randomly oriented, noninteracting planar microcracks of specified geometry but of variable size. This methodology has been extended to predict the multiaxial strength response of transversely isotropic brittle materials, including polymer matrix composites (PMCs), by considering (1) flaw-orientation anisotropy, whereby a preexisting microcrack has a higher likelihood of being oriented in one direction over another direction, and (2) critical strength, or K (sub Ic) orientation anisotropy, whereby the level of critical strength or fracture toughness for mode I crack propagation, K (sub Ic), changes with regard to the orientation of the microstructure. In this report, results from finite element analysis of a fiber-reinforced-matrix unit cell were used with the unit-sphere model to predict the biaxial strength response of a unidirectional PMC previously reported from the World-Wide Failure Exercise. Results for nuclear-grade graphite materials under biaxial loading are also shown for comparison. This effort was successful in predicting the multiaxial strength response for the chosen problems. Findings regarding stress-state interactions and failure modes also are provided.

Nemeth, Noel, N.↗

Glove-Enabled Computer Operations (GECO): Design and Testing of an Extravehicular Activity Glove Adapted for Human-Computer Interface

The Glove-Enabled Computer Operations (GECO) system enables an extravehicular activity (EVA) glove to be dual-purposed as a human-computer interface device. This paper describes the design and human participant testing of a right-handed GECO glove in a pressurized glove box. As part of an investigation into the usability of the GECO system for EVA data entry, twenty participants were asked to complete activities including (1) a Simon Says Games in which they attempted to duplicate random sequences of targeted finger strikes and (2) a Text Entry activity in which they used the GECO glove to enter target phrases in two different virtual keyboard modes. In a within-subjects design, both activities were performed both with and without vibrotactile feedback. Participants mean accuracies in correctly generating finger strikes with the pressurized glove were surprisingly high, both with and without the benefit of tactile feedback. Five of the subjects achieved mean accuracies exceeding 99 in both conditions. In Text Entry, tactile feedback provided a statistically significant performance benefit, quantified by characters entered per minute, as well as reduction in error rate. Secondary analyses of responses to a NASA Task Loader Index (TLX) subjective workload assessments reveal a benefit for tactile feedback in GECO glove use for data entry. This first-ever investigation of employment of a pressurized EVA glove for human-computer interface opens up a wide range of future applications, including text chat communications, manipulation of procedureschecklists, cataloguingannotating images, scientific note taking, human-robot interaction, and control of suit andor other EVA systems.

Adams, Richard J.↗

Glove-Enabled Computer Operations (GECO): Design and Testing of an Extravehicular Activity Glove Adapted for Human-Computer Interface

The Glove-Enabled Computer Operations (GECO) system enables an extravehicular activity (EVA) glove to be dual-purposed as a human-computer interface device. This paper describes the design and human participant testing of a right-handed GECO glove in a pressurized glove box. As part of an investigation into the usability of the GECO system for EVA data entry, twenty participants were asked to complete activities including (1) a Simon Says Games in which they attempted to duplicate random sequences of targeted finger strikes and (2) a Text Entry activity in which they used the GECO glove to enter target phrases in two different virtual keyboard modes. In a within-subjects design, both activities were performed both with and without vibrotactile feedback. Participants' mean accuracies in correctly generating finger strikes with the pressurized glove were surprisingly high, both with and without the benefit of tactile feedback. Five of the subjects achieved mean accuracies exceeding 99% in both conditions. In Text Entry, tactile feedback provided a statistically significant performance benefit, quantified by characters entered per minute, as well as reduction in error rate. Secondary analyses of responses to a NASA Task Loader Index (TLX) subjective workload assessments reveal a benefit for tactile feedback in GECO glove use for data entry. This first-ever investigation of employment of a pressurized EVA glove for human-computer interface opens up a wide range of future applications, including text "chat" communications, manipulation of procedures/checklists, cataloguing/annotating images, scientific note taking, human-robot interaction, and control of suit and/or other EVA systems.

Adams, Richard J.↗

The importance of independent replication of treatments in plant science

Among other foci, the Journal of Experimental Botany aims to advance understanding of plant environment interactions, including abiotic stress, mineral nutrition, and the response of plants to global change. Advancing understanding in these areas often requires manipulation of the growth environment. Experiments range in scale and include: growth chambers (Yiotis et al., 2020), glasshouses (Rho et al., 2019), whole plant chambers (Sharwood et al., 2017), and open-air field experiments (Ruiz-Vera et al., 2020). The application of environmental manipulations can be constrained by available resources, operational costs, and time. In designing experiments, researchers attempt to balance logistical challenges with the desire to address scientific questions and maximize statistical power.

54 ENVIRONMENTAL SCIENCES↗

The Joint Assimilation of Remotely Sensed Leaf Area Index and Surface Soil Moisture into a Land Surface Model

This work tests the hypothesis that jointly assimilating satellite observations of leaf area index and surface soil moisture into a land surface model improves the estimation of land vegetation and water variables. An Ensemble Kalman Filter is used to test such hypothesis across the Contiguous United States during April 2015 – December 2018. The performance of the proposed methodology is assessed for several modeled vegetation and water variables (evapotranspiration, net ecosystem exchange, and soil moisture) in terms of random errors and anomaly correlation coefficients against a set of independent validation datasets (i.e., Global Land Evaporation Amsterdam Model, FLUXCOM, and International Soil Moisture Network). Results show that the assimilation of leaf area index mostly improves the estimation of evapotranspiration and net ecosystem exchange, whereas the assimilation of surface soil moisture alone improves surface soil moisture content, especially in the western US, in terms of both root mean squared error and anomaly correlation coefficient. The joint assimilation of vegetation and soil moisture information combines the results of individual vegetation and soil moisture assimilations and reduces errors (and increases correlations with the reference datasets) in evapotranspiration, net ecosystem exchange, and surface soil moisture simulated by the land surface model. However, because soil moisture satellite observations only provide information on the water content in the top 5 cm of the soil column, the impact of the proposed data assimilation technique on root zone soil moisture is limited. This work moves one step forward in the direction of improving our estimation and understanding of land surface interactions using a multi-variate data assimilation approach, which can be particularly useful in regions of the world where ground observations are sparse or missing altogether.

Data assimilation↗

Liquid-to-gas transfer of sodium in a liquid anode glow discharge

Solute transfer from the liquid to the gas phase in plasmas has been widely used in analytical chemistry often focused on liquid cathode discharges. In this study, we report on time-resolved absolute density measurements of sodium in a glow discharge plasma with a NaCl solution anode using two-photon absorption laser induced fluorescence. The averaged maximum sodium density is found to be similar to the electron density in the discharge consistent with plasma-induced vapor generation. Random sodium flashes were also observed extending into the plasma afterglow and are attributed to sodium particles falling back on the liquid surface. Cataphoresis can explain the location of the maximum sodium density near the cathode side of the electrode gap rather than near the liquid anode region.

chemical vapor generation↗

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗

Direct numerical simulation of open-channel flow over a heterogeneous particle bed at low relative submergence

In this article, this study investigates turbulent open-channel flows over beds of irregularly arranged particles, using direct numerical simulations at a friction Reynolds number of Re T = 300⁠. Two distinct cases are examined: a polydisperse bed (P800) composed of multiple layers of randomly distributed spheres of varying sizes, and a monodisperse bed (M1015) formed by a random distribution of uniform sized spheres, with a bottommost single layer of varied-sized particles to introduce realistic randomness. Our investigation unveils a rich network of low- and high-speed streaks within the flow field, exhibiting distinctive behaviors in different bed configurations. The P800 case presents a poorly organized flow pattern induced by the varied particle sizes and arrangements, while the M1015 case shows a more regular flow pattern, marked by larger streaks. We also observe that total wall shear stress is substantially influenced by surface roughness-induced drag, extending beyond the effects documented in existing studies of open-channel flows. The present study reveals intricate secondary flow patterns over irregular particle beds. Large-scale circulations are discerned around particle crests in the P800 case and localized circulations with increased turbulence in the M1015 case. Furthermore, analysis of Reynolds stress tensor components indicates that roughness disrupts coherent turbulent eddies, consequently mitigating peak stress. We quantify correlations between drag force and local fluid velocity fluctuations. Notably, a larger deviation in drag is observed in the P800 case compared to M1015, accentuating the influence of particle size and distribution on fluid–particle interactions.

42 ENGINEERING↗

Numerical Test of Analytical Theories for Perpendicular Diffusion in Small Kubo Number Turbulence

In the literature, one can find various analytical theories for perpendicular diffusion of energetic particles interacting with magnetic turbulence. Besides quasi-linear theory, there are different versions of the nonlinear guiding center (NLGC) theory and the unified nonlinear transport (UNLT) theory. For turbulence with high Kubo numbers, such as two-dimensional turbulence or noisy reduced magnetohydrodynamic turbulence, the aforementioned nonlinear theories provide similar results. For slab and small Kubo number turbulence, however, this is not the case. In the current paper, we compare different linear and nonlinear theories with each other and test-particle simulations for a noisy slab model corresponding to small Kubo number turbulence. We show that UNLT theory agrees very well with all performed test-particle simulations. In the limit of long parallel mean free paths, the perpendicular mean free path approaches asymptotically the quasi-linear limit as predicted by the UNLT theory. For short parallel mean free paths we find a Rechester and Rosenbluth type of scaling as predicted by UNLT theory as well. The original NLGC theory disagrees with all performed simulations regardless what the parallel mean free path is. The random ballistic interpretation of the NLGC theory agrees much better with the simulations, but compared to UNLT theory the agreement is inferior. We conclude that for this type of small Kubo number turbulence, only the latter theory allows for an accurate description of perpendicular diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗

Remote Linkages to Anomalous Winter Atmospheric Ridging over the Northeastern Pacific

Severe drought in California between 2013 and 2016 has been linked to the multiyear persistence of anomalously high atmospheric pressure over the northeastern Pacific Ocean, which deflected the Pacific storm track northward and suppressed regional precipitation during California's winter 'rainy season.' Multiple hypotheses have emerged regarding why this high pressure ridge near the west coast of North America was so resilient-including unusual sea surface temperature patterns in the Pacific Ocean, reductions in Arctic sea ice, random atmospheric variability, or some combination thereof. Here we explore relationships between previously documented atmospheric conditions over the North Pacific and several potential remote oceanic and cryospheric influences using both observational data and a large ensemble of climate model simulations. Our results suggest that persistent wintertime atmospheric ridging similar to that implicated in California's 2013-2016 drought can at least partially be linked to unusual Pacific sea surface temperatures, and that Pacific Ocean conditions may offer some degree of cool-season foresight in this region despite the presence of substantial internal variability.

California drought;Atmospheric ridging;Ocean-atmos↗

Contributions of vegetation heterogeneity within tower footprint to CO 2 flux estimations through graph neural network modeling

Net ecosystem exchange of CO 2 (Fc) measured directly by eddy covariance towers is based on various assumptions, including large, flat and homogenous land cover type. In reality, often a tower site is not large enough for flux measurements, and landscapes consist of patches of different land cover types within the flux footprint. In addition, some portions of fluxes are contributed by different cover types when a footprint exceeds the size of the target ecosystem. The contributions of non-dominant patches to Fc are often ignored. Here, in this study, we propose a novel integrated modeling framework that combines random forest (RF) and XGBoost with a residual correction module based on a deep graph convolutional network (DeeperGCN) to simulate Fc for seven flux measurement sites in southwest Michigan. High-resolution remote sensing vegetation indices, soil properties, meteorological variables, and footprint-weighted spatial features were used as model inputs at three spatial resolutions (10 m, 20 m, 30 m), and their importance in predicting Fc with DeeperGCN was assessed. We found that residual correction using DeeperGCN significantly improved prediction accuracy, with the R 2 increasing from 0.9098 to 0.9479 for RF and from 0.9235 to 0.9433 for XGBoost. At site level, the maximum improvement in R 2 reached 0.1617. Paired t-tests confirmed that these improvements were statistically significant (p < 0.05). Among all predictors, leaf area index and incoming shortwave radiation emerged as the dominant drivers of spatial residual variation, followed by precipitation, relative humidity, and selected vegetation indices. The 20 m resolution yielded the best balance between model performance and computational efficiency. In conclusion, our modeling framework effectively captures both spatial heterogeneity and nonlinear interactions, offering a robust solution for spatially explicit flux modeling in structurally diverse ecosystems beyond the study sites.

footprint model↗

Structure Tuning, Strong Second Harmonic Generation Response, and High Optical Stability of the Polar Semiconductors Na 1- x K x AsQ 2

We report the mixed cation compounds Na 1-x K x AsSe 2 (x = 0.8, 0.65, 0.5) and Na 0.1 K 0.9 AsS 2 crystallize in the polar noncentrosymmetric space group Cc. The AAsQ(2) (A = alkali metals, Q = S, Se) family features one-dimensional (1D) 1 / ∞ [AQ 2 - ] chains comprising corner-sharing pyramidal AQ 3 units in which the packing of these chains is dependent on the alkali metals. The parallel 1 / ∞ [AQ(2) - ] chains interact via short As ∙∙∙Se contacts, which increase in length when the fraction of K atoms is increased. The increase in the As ∙∙∙Se interchain distance increases the band gap from 1.75 eV in γ-NaAsSe 2 to 2.01 eV in Na 0.35 K 0.65 AsSe 2 , 2.07 eV in Na 0.2 K 0.8 AsSe 2 , and 2.18 eV in Na 0.1 K 0.9 AsS 2 . The Na 1-x K x AsSe 2 (x = 0.8, 0.65) compounds melt congruently at approximately 316 °C. Wavelength-dependent second harmonic generation (SHG) measurements on powder samples of Na 1-x K x AsSe 2 (x = 0.8, 0.65, 0.5) and Na 0.1 K 0.9 AsS 2 suggest that Na 0.2 K 0.8 AsSe 2 and Na 0.1 K 0.9 AsS 2 have the highest SHG response and exhibit significantly higher laser-induced damage thresholds (LIDTs). Theoretical SHG calculations on Na 0.5 K 0.5 AsSe 2 confirm its SHG response with the highest value of d 33 = 22.5 pm/V χ 333 (2) = 45.0 pm/V). The effective nonlinearity for a randomly oriented powder is calculated to be d eff = 18.9 pm/V χ eff (2) = 37.8 pm/V), which is consistent with the experimentally obtained value of d eff = 16.5 pm/V χ eff (2) = 33.0 pm/V). Three-photon absorption is the dominant mechanism for the optical breakdown of the compounds under intense excitation at 1580 nm, with Na 0.2 K 0.8 AsSe 2 exhibiting the highest stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of Monte Carlo techniques to optimization of high-energy beam transport in a stochastic environment

An algorithm employing a modified sequential random perturbation, or creeping random search, was applied to the problem of optimizing the parameters of a high-energy beam transport system. The stochastic solution of the mathematical model for first-order magnetic-field expansion allows the inclusion of state-variable constraints, and the inclusion of parameter constraints allowed by the method of algorithm application eliminates the possibility of infeasible solutions. The mathematical model and the algorithm were programmed for a real-time simulation facility; thus, two important features are provided to the beam designer: (1) a strong degree of man-machine communication (even to the extent of bypassing the algorithm and applying analog-matching techniques), and (2) extensive graphics for displaying information concerning both algorithm operation and transport-system behavior. Chromatic aberration was also included in the mathematical model and in the optimization process. Results presented show this method as yielding better solutions (in terms of resolutions) to the particular problem than those of a standard analog program as well as demonstrating flexibility, in terms of elements, constraints, and chromatic aberration, allowed by user interaction with both the algorithm and the stochastic model. Example of slit usage and a limited comparison of predicted results and actual results obtained with a 600 MeV cyclotron are given.

Parrish, R. V.↗