Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “level set methods”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Dynamic Server-Based KML Code Generator Method for Level-of-Detail Traversal of Geospatial Data

Web-based geospatial client applications such as Google Earth and NASA World Wind must listen to data requests, access appropriate stored data, and compile a data response to the requesting client application. This process occurs repeatedly to support multiple client requests and application instances. Newer Web-based geospatial clients also provide user-interactive functionality that is dependent on fast and efficient server responses. With massively large datasets, server-client interaction can become severely impeded because the server must determine the best way to assemble data to meet the client applications request. In client applications such as Google Earth, the user interactively wanders through the data using visually guided panning and zooming actions. With these actions, the client application is continually issuing data requests to the server without knowledge of the server s data structure or extraction/assembly paradigm. A method for efficiently controlling the networked access of a Web-based geospatial browser to server-based datasets in particular, massively sized datasets has been developed. The method specifically uses the Keyhole Markup Language (KML), an Open Geospatial Consortium (OGS) standard used by Google Earth and other KML-compliant geospatial client applications. The innovation is based on establishing a dynamic cascading KML strategy that is initiated by a KML launch file provided by a data server host to a Google Earth or similar KMLcompliant geospatial client application user. Upon execution, the launch KML code issues a request for image data covering an initial geographic region. The server responds with the requested data along with subsequent dynamically generated KML code that directs the client application to make follow-on requests for higher level of detail (LOD) imagery to replace the initial imagery as the user navigates into the dataset. The approach provides an efficient data traversal path and mechanism that can be flexibly established for any dataset regardless of size or other characteristics. The method yields significant improvements in userinteractive geospatial client and data server interaction and associated network bandwidth requirements. The innovation uses a C- or PHP-code-like grammar that provides a high degree of processing flexibility. A set of language lexer and parser elements is provided that offers a complete language grammar for writing and executing language directives. A script is wrapped and passed to the geospatial data server by a client application as a component of a standard KML-compliant statement. The approach provides an efficient means for a geospatial client application to request server preprocessing of data prior to client delivery. Data is structured in a quadtree format. As the user zooms into the dataset, geographic regions are subdivided into four child regions. Conversely, as the user zooms out, four child regions collapse into a single, lower-LOD region. The approach provides an efficient data traversal path and mechanism that can be flexibly established for any dataset regardless of size or other characteristics.

Baxes, Gregory↗
DOE OSTI2013

Solution of resource allocation problem for identification of cost-effective measures to reduce nuclear proliferation risks

This report presents a methodology of selection of cost-effective measures to reduce nuclear proliferation risks. The methodology relies on a graded security model used in practice in different applications. The method is based on the controlled finite Markov chain approach set in combination with discrete dynamic programming and MCDM (Multi Criteria Decision Making) techniques that enables the expert to select the cost-effective measures to reduce nuclear proliferation risks depending on availability of resources. The analysis performed with different number of possible measures confirms the conclusions that the implementation of extra-large costs may not produce the required effect, and the increase in resources above a certain level does not appear sensitive. Diversification in improving the effectiveness of other measures seems more rational and efficient for the whole system than the unlimited improvement of the effectiveness of only one measure.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2019

Toward Fast and Reliable Potential Energy Surfaces for Metallic Pt Clusters by Hierarchical Delta Neural Networks

Data-driven machine learning force fields (MLF) are more and more popular in atomistic simulations, and exploit machine learning methods to predict energies and forces for unknown structures based on the knowledge learned from an existing reference database. The latter usually comes from density functional theory calculations. One main drawback of MLFs is that physical laws are not incorporated in the machine learning models and instead, MLFs are designed to be very flexible to simulate complex quantum chemistry potential energy surface (PES). In general, MLFs have poor transferability, and hence a very large trainset is required to span all the target feature space to get a reliable MLF. This procedure becomes more troublesome when the PES is complicated, with a large number of degrees of freedom, in which building a large database is inevitable and very expensive, especially when accurate but costly exchange-correlation functionals have to be used. In this manuscript, we exploit a high dimensional neural network potential (HDNNP) on Pt clusters of size 6 to 20 as one example. Furthermore, our standard level of energy calculation is DFT GGA (PBE) using a plane wave basis set. We introduce an approximate but fast level with the PBE functional and a minimal atomic orbital basis set, then a more accurate but expensive level, using a hybrid functional or non-local vdw functional and a plane wave basis set, is reliably predicted by learning the difference with HDNNP. Findings show that such a differential approach (named ΔHDNNP) can deliver very accurate predictions (error < 10 meV/atom) in reference to converged basis set energies as well as more accurate but expensive xc functional. The overall speedup can be as large as 900 for 20 atom Pt cluster. More importantly, ΔHDNNP shows much better transferability due to the intrinsic smoothness of delta potential energy surface, and accordingly one can use much smaller trainset data to obtain better accuracy than the conventional HDNNP. A multi-layer ΔHDNNP is thus proposed to obtain very accurate predictions versus expensive non-local vdW functional calculations in which the required trainset is further reduced. The approach can be easily generalized to any other machine learning methods and opens a path to study the structure and dynamics of Pt clusters and nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A novel method for quantifying irradiation damage in nuclear graphite using Raman spectroscopy

Raman spectroscopy has long been used in studying irradiation damage in carbon/graphite materials. It is however unclear if the measurements from different types of materials are directly comparable. Further, decoupling the contribution from irradiation temperature and fluence to the total damage in nuclear graphite possessing irradiation damage gradient is currently not readily feasible as it requires a complete set of test reactor irradiation experiments over the relevant temperature and fluence range. A novel methodology has therefore been proposed and developed to quantify total irradiation damage evolution at crystal level based on graphite Raman G-band position shift. Specifically, G-band positions derived from a large number of spectra collected from the fractured surface of microfine-grained POCO ZXF-5Q graphite possessing proton irradiation damage gradient were plotted with open literature Raman data on HOPG, BEPO (AXGP graphite), PCEA and IG-110 graphite as a function of dpa. The total damage level within 2σ beam radius in this POCO graphite was estimated to be equivalent to ~2–5 dpa at ~350–370 °C. Derived G-band positions were then mapped to the three-stage amorphization trajectory model indicating beam centre area has entered the second stage, i.e., transitioning from nanocrystalline graphite into amorphous carbon. G-band position relative shift (ΔG) curve with a ‘turn-around’ peak as a function of total damage can be used for in-service lifetime prediction. Importantly, the developed methodology has the potential to ‘unify’ total damage levels across different grades of nuclear graphite caused by ion, neutron and proton irradiation at different temperatures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Couch height–based patient setup for abdominal radiation therapy

There are 2 methods commonly used for patient positioning in the anterior-posterior (A-P) direction: one is the skin mark patient setup method (SMPS) and the other is the couch height–based patient setup method (CHPS). This study compared the setup accuracy of these 2 methods for abdominal radiation therapy. The enrollment for this study comprised 23 patients with pancreatic cancer. For treatments (539 sessions), patients were set up by using isocenter skin marks and thereafter treatment couch was shifted so that the distance between the isocenter and the upper side of the treatment couch was equal to that indicated on the computed tomographic (CT) image. Setup deviation in the A-P direction for CHPS was measured by matching the spine of the digitally reconstructed radiograph (DRR) of a lateral beam at simulation with that of the corresponding time-integrated electronic portal image. For SMPS with no correction (SMPS/NC), setup deviation was calculated based on the couch-level difference between SMPS and CHPS. SMPS/NC was corrected using 2 off-line correction protocols: no action level (SMPS/NAL) and extended NAL (SMPS/eNAL) protocols. Margins to compensate for deviations were calculated using the Stroom formula. A-P deviation > 5 mm was observed in 17% of SMPS/NC, 4% of SMPS/NAL, and 4% of SMPS/eNAL sessions but only in one CHPS session. For SMPS/NC, 7 patients (30%) showed deviations at an increasing rate of > 0.1 mm/fraction, but for CHPS, no such trend was observed. The standard deviations (SDs) of systematic error (Σ) were 2.6, 1.4, 0.6, and 0.8 mm and the root mean squares of random error (σ) were 2.1, 2.6, 2.7, and 0.9 mm for SMPS/NC, SMPS/NAL, SMPS/eNAL, and CHPS, respectively. Margins to compensate for the deviations were wide for SMPS/NC (6.7 mm), smaller for SMPS/NAL (4.6 mm) and SMPS/eNAL (3.1 mm), and smallest for CHPS (2.2 mm). Achieving better setup with smaller margins, CHPS appears to be a reproducible method for abdominal patient setup.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2016

Additive Manufacturing, Design, Testing, and Fabrication: A Full Engineering Experience at JSC

I worked on several projects this term. While most projects involved additive manufacturing, I was also involved with two design projects, two testing projects, and a fabrication project. The primary mentor for these was Richard Hagen. Secondary mentors were Hai Nguyen, Khadijah Shariff, and fabrication training from James Brown. Overall, my experience at JSC has been successful and what I have learned will continue to help me in my engineering education and profession long after I leave. My 3D printing projects ranged from less than a 1 cubic centimeter to about 1 cubic foot and involved several printers using different printing technologies. It was exciting to become familiar with printing technologies such as industrial grade FDM (Fused Deposition Modeling), the relatively new SLA (Stereolithography), and PolyJet. My primary duty with the FDM printers was to model parts that came in from various sources to print effectively and efficiently. Using methods my mentor taught me and the Stratasys Insight software, I was able to minimize imperfections, hasten build time, improve strength for specific forces (tensile, shear, etc...), and reduce likelihood of a print-failure. Also using FDM, I learned how to repair a part after it was printed. This is done by using a special kind of glue that chemically melts the two faces of plastic parts together to form a fused interface. My first goal with SLA technology was to bring the printer back to operational readiness. In becoming familiar with the Pegasus SLA printer, I researched the leveling, laser settings, and different vats to hold liquid material. With this research, I was successfully able to bring the Pegasus back online and have successfully printed multiple sample parts as well as functional parts. My experience with PolyJet technology has been focused on an understanding of the abilities/limits, costs, and the maintenance for daily use. Still upcoming will be experience with using a composite printer that uses FDM technology to print plastic while laying an internal filament of Kevlar or carbon-fiber inside the printed material. It has been incredible being exposed to this range of technologies and I feel very fortunate to be ready for virtually any kind of printing technology I come across in the future. Design work played a part in my internship this term as well. Working with Hai Nguyen, I was able to design a set of testing tips and a test frame for use with an Arcjet. The testing tips will be made of several different materials that will possibly be used in the heat shield of the Orion space craft. These designs included technical drawings that were presented to the fabrication shop. The frame design was created from 80/20 (a popular brand of frame construction equipment) and included an order form with pricing for fabrication. An independent design was also done for the virtual reality lab. This design was to create a hand-controller based on a previous design. This final design was sent directly to a 3D printer without technical drawings. Overall, my design work has given me experience with using 80/20, helped improve my CAD (Computer Aided Design) proficiency, and increased my knowledge of how to set up technical drawings for fabrication. The final role I have played in this internship has been to assist with testing of the inflatable technology materials working Khadijah Shariff. I began the internship assisting with permeability testing with the initial plan to continue the testing independently after training. Unfortunately, the testing apparatus suffered a technical failure and had funding pulled which cancelled that portion of the project. Further testing with inflatable technology continued with tensile testing of various stitching methods. This testing took place over a two-days and concluded successfully. Final testing was to be more tensile testing but of straps used to connect various inflatable sections. Unfortunately, the needed grips for the tensile tests could not be located and put the testing on hold. It is possible this round of testing will still take place by the end of the internship if the grips can be found. Overall, this portion of the internship has helped me become familiar with one kind of permeability test as well as a popular tensile/compression testing machine. Finally, I also had the chance to be trained using a metal lathe for making very small tips for a soldering iron. These tips will be used to melt threaded brass inserts into 3D printed plastic. Working with James Brown, I was able to successfully machine a brass rod down to as little as 0.064 inches plus or minus 0.001 inches. It was very rewarding to learn how to best use the machine and become familiar with a skill that will undoubtedly be used again in the future. I have been told by several professional engineers that learning to use a lathe and mill will be invaluable skills in the field. This knocks 50 percent of that goal off and I look forward to learning the mill at some point in the future. As is apparent with this list of projects, my internship was not focused on a single over-arching goal. Instead, I was able to gain experience in a myriad of very different areas. I feel like my time here was spent bouncing from one project to the next. Though sometimes difficult to switch gears, it was very rewarding to be a part of so much in so little time. My career and education will be positively impacted by what I have learned at JSC. My experience with 3D printing has improved my ability to handle many issues that may come up in the future with multiple different technologies. The design work I took part in, especially creating technical drawings, will help me better present designs to any engineer or shop I will encounter. My testing experience has helped me become familiar with a popular kind of tensile test machine that will likely be similar to the kinds I will encounter in the future. Finally, my experience with fabrication has given me a rare opportunity, as an engineer, to take part in the fabrication of a part. This experience will help me better tailor my future designs for the manufacturing process and has given me an appreciation for detailed/delicate machining work. My experience at JSC has been successful and will continue to assist me for a long time within the engineering field.

Zusack, Steven↗
DOE OSTI2015

CONSTRAINTS FROM ASYMMETRIC HEATING: INVESTIGATING THE EPSILON AURIGAE DISK

Epsilon Aurigae is a long-period eclipsing binary that likely contains an F0Ia star and a circumstellar disk enshrouding a hidden companion, assumed to be a main-sequence B star. High uncertainty in its parallax has kept the evolutionary status of the system in question and, hence, the true nature of each component. This unknown, as well as the absence of solid state spectral features in the infrared, requires an investigation of a wide parameter space by means of both analytic and Monte Carlo radiative transfer (MCRT) methods. The first MCRT models of epsilon Aurigae that include all three system components are presented here. We seek additional system parameter constraints by melding analytic approximations with MCRT outputs (e.g., dust temperatures) on a first-order level. The MCRT models investigate the effects of various parameters on the disk-edge temperatures; these include two distances, three particle size distributions, three compositions, and two disk masses, resulting in 36 independent models. Specifically, the MCRT temperatures permit analytic calculations of effective heating and cooling curves along the disk edge. These are used to calculate representative observed fluxes and corresponding temperatures. This novel application of thermal properties provides the basis for utilization of other binary systems containing disks. We find degeneracies in the model fits for the various parameter sets. However, the results show a preference for a carbon disk with particle size distributions ≥10 μm. Additionally, a linear correlation between the MCRT noon and basal temperatures serves as a tool for effectively eliminating portions of the parameter space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗
DOE OSTI2012

Impact of Schedule Duration on Head and Neck Radiotherapy: Accelerated Tumor Repopulation Versus Compensatory Mucosal Proliferation

Purpose: To determine how modelled maximum tumor control rates, achievable without exceeding mucositis tolerance (tcp{sub max-early}) vary with schedule duration for head and neck squamous cell carcinoma (HNSCC). Methods and materials: Using maximum-likelihood techniques, we have fitted a range of tcp models to two HNSCC datasets (Withers' and British Institute of Radiology [BIR]), characterizing the dependence of tcp on duration and equivalent dose in 2 Gy fractions (EQD{sub 2}). Models likely to best describe future data have been selected using the Akaike information criterion (AIC) and its quasi-AIC extension to overdispersed data. Setting EQD{sub 2}s in the selected tcp models to levels just tolerable for mucositis, we have plotted tcp{sub max-early} against schedule duration. Results: While BIR dataset tcp fits describe dose levels isoeffective for tumor control as rising significantly with schedule protraction, indicative of accelerated tumor repopulation, repopulation terms in fits to Withers' dataset do not reach significance after accounting for overdispersion of the data. The tcp{sub max-early} curves calculated from tcp fits to the overall Withers' and BIR datasets rise by 8% and 0-4%, respectively, between 20 and 50 days duration; likewise, tcp{sub max-early} curves calculated for stage-specific cohorts also generally rise slowly with increasing duration. However none of the increases in tcp{sub max-early} calculated from the overall or stage-specific fits reach significance. Conclusions: Local control rates modeled for treatments which lie just within mucosal tolerance rise slowly but insignificantly with increasing schedule length. This finding suggests that whereas useful gains may be made by accelerating unnecessarily slow schedules until they approach early reaction tolerance, little is achieved by shortening schedules further while reducing doses to remain within mucosal tolerance, an approach that may slightly worsen outcomes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Feasibility study of prompt gamma neutron activation for NDT measurement of moisture in stone and brick

The conservation of stone and brick architecture or sculpture often involves damage caused by moisture. The feasibility of a NDT method based on prompt gamma neutron activation (PGNA) for measuring the element hydrogen as an indication of water is being evaluated. This includes systematic characterization of the lithology and physical properties of seven building stones and one brick type used in the buildings of the Smithsonian Institution in Washington, D.C. To determine the required dynamic range of the NDT method, moisture-related properties were measured by standard methods. Cold neutron PGNA was also used to determine chemically bound water (CBW) content. The CBW does not damage porous masonry, but creates an H background that defines the minimum level of detection of damaging moisture. The CBW was on the order of 0.5% for all the stones. This rules out the measurement of hygric processes in all of the stones and hydric processed for the stones with fine scale pore-size distributions The upper bound of moisture content, set by porosity through water immersion, was on the order of 5%. The dynamic range is about 10–20. The H count rates were roughly 1–3 cps. Taking into account differences in neutron energies and fluxes and sample volume between cold PGNA and a portable PGNA instrument, it appears that it is feasible to apply PGNA in the field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

A support vector machine classifier reduces interscanner variation in the HRCT classification of regional disease pattern in diffuse lung disease: Comparison to a Bayesian classifier

Purpose: To investigate the effect of using different computed tomography (CT) scanners on the accuracy of high-resolution CT (HRCT) images in classifying regional disease patterns in patients with diffuse lung disease, support vector machine (SVM) and Bayesian classifiers were applied to multicenter data. Methods: Two experienced radiologists marked sets of 600 rectangular 20 Multiplication-Sign 20 pixel regions of interest (ROIs) on HRCT images obtained from two scanners (GE and Siemens), including 100 ROIs for each of local patterns of lungs-normal lung and five of regional pulmonary disease patterns (ground-glass opacity, reticular opacity, honeycombing, emphysema, and consolidation). Each ROI was assessed using 22 quantitative features belonging to one of the following descriptors: histogram, gradient, run-length, gray level co-occurrence matrix, low-attenuation area cluster, and top-hat transform. For automatic classification, a Bayesian classifier and a SVM classifier were compared under three different conditions. First, classification accuracies were estimated using data from each scanner. Next, data from the GE and Siemens scanners were used for training and testing, respectively, and vice versa. Finally, all ROI data were integrated regardless of the scanner type and were then trained and tested together. All experiments were performed based on forward feature selection and fivefold cross-validation with 20 repetitions. Results: For each scanner, better classification accuracies were achieved with the SVM classifier than the Bayesian classifier (92% and 82%, respectively, for the GE scanner; and 92% and 86%, respectively, for the Siemens scanner). The classification accuracies were 82%/72% for training with GE data and testing with Siemens data, and 79%/72% for the reverse. The use of training and test data obtained from the HRCT images of different scanners lowered the classification accuracy compared to the use of HRCT images from the same scanner. For integrated ROI data obtained from both scanners, the classification accuracies with the SVM and Bayesian classifiers were 92% and 77%, respectively. The selected features resulting from the classification process differed by scanner, with more features included for the classification of the integrated HRCT data than for the classification of the HRCT data from each scanner. For the integrated data, consisting of HRCT images of both scanners, the classification accuracy based on the SVM was statistically similar to the accuracy of the data obtained from each scanner. However, the classification accuracy of the integrated data using the Bayesian classifier was significantly lower than the classification accuracy of the ROI data of each scanner. Conclusions: The use of an integrated dataset along with a SVM classifier rather than a Bayesian classifier has benefits in terms of the classification accuracy of HRCT images acquired with more than one scanner. This finding is of relevance in studies involving large number of images, as is the case in a multicenter trial with different scanners.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-206-08: How Should Focal Spot Be Chosen for Optimized CT Imaging with Dose Modulation?

Purpose: To choose the preferred focal spot for achieving optimized CT image quality with balanced tube heating considerations. Methods: An anthropomorphic pelvic phantom was scanned using a GE Discovery CT750 HD at 120 and 140kVp, 0.8s rotation time, and pitch of 0.984. “Smart mA” was enabled to simulate a routine abdomen–pelvis CT scan. Permissible mA values at 120 and 140 kVp were obtained from the Serial Load Rating table (for mimicking a busy CT clinical operation) in the scanner Technical Reference Manual. At each kVp station and two Noise Index levels, the mA Upper Limit was set above/below the permissible mA values. Scanned mA values and focal spot (FS) used were obtained from the DICOM header of each image, and the FS-mA relationship was analyzed. For visual confirmation beyond recorded FS information, a CatPhan with a fat-ring attached for mimicking patient size/shape was scanned at 120kVp. A group of radiologists/physicists compared a pair of CatPhan images qualitatively. Lastly, a number of patient cases were evaluated to confirm the FS-mA relationship. Results: When preset Upper Limit values were above the permissible mA values, the Large FS (labeled 1.2) was used in scans, even if the maximum scanned mA values were much lower than the permissible values at both 120 and 140 kVp. Otherwise the Small FS (labeled 0.7) was used. Visual evaluation of the high contrast module of CatPhan and additional analysis of patient cases further confirmed that the preset Upper Limit determines which focal spot is to be used, not the actual maximum mA value to be scanned. Conclusion: Specific FS can be selected by setting up appropriate mA Upper Limit in a protocol. CT protocols could be optimized by selecting appropriate FS for improving patient image quality, especially benefiting the small size and pediatric patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Search for the Standard Model Higgs boson produced in association with a W Boson in the isolated-track charged-lepton channel using the Collider Detector at Fermilab

The Higgs boson is the only elementary particle predicted by the Standard Model (SM) that has not yet been observed experimentally. If it exists, it explains the spontaneous electroweak symmetry breaking and the origin of mass for gauge bosons and fermions. We test the validity of the SM by performing a search for the associated production of a Higgs boson and a W boson in the channel where the Higgs boson decays to a bottom-antibottom quark pair and the W boson decays to a charged lepton and a neutrino (the WH channel). We study a dataset of proton-antiproton collisions at a centre-of-mass energy √s = 1.96 TeV provided by the Tevatron accelerator, corresponding to an integrated luminosity of 5.7 fb -1 , and recorded using the Collider Detector at Fermilab (CDF).We select events consistent with the signature of exactly one charged lepton (electron or muon), missing transverse energy due to the undetected neutrino (MET) and two collimated streams of particles (jets), at least one of which is required to be identified as originating from a bottom quark. We improve the discrimination of Higgs signal from backgrounds through the use of an artificial neural network. Using a Bayesian statistical inference approach, we set for each hypothetical Higgs boson mass in the range 100-150 GeV/c 2 with 5 GeV/c 2 increments a 95% credibility level (CL) upper limit on the ratio between the Higgs production cross section times branching fraction and the SM prediction. Our main original contributions are the addition of a novel charged lepton reconstruction algorithm with looser requirements (ISOTRK) with respect the electron or muon tight criteria (TIGHT), as well as the introduction of a novel trigger-combination method that allows to maximize the event yield while avoiding trigger correlations and that is used for the ISOTRK category. The ISOTRK candidate is a high-transverse-momentum good-quality track isolated from other activity in the tracking system and not required to match a calorimeter cluster, as for a tight electron candidate, or an energy deposit in the muon detector, as for a tight muon candidate. The ISOTRK category recovers real charged leptons that otherwise would be lost in the non-instrumented regions of the detector. This allows the reconstruction of more W boson candidates, which in turn increases the number of reconstructed WH signal candidate events, and therefore improves the sensitivity of the WH search. For the TIGHT charged lepton categories, we employ charged-lepton-dedicated triggers to improve the rate of WH signal acceptance during data taking. Since there is no ISOTRK-dedicated trigger at CDF, for the ISOTRK charged lepton category we employ three MET-plus-jets-based triggers. For each trigger we first identify the jet selection where the trigger efficiency is flat with respect to jet information (transverse energy and direction of motion in the transverse plane for the two jets in the event) and then we parametrize the trigger efficiency as a function of trigger MET. On an event-by-event basis, for each trigger we compute a trigger efficiency as a function of trigger parametrization, trigger MET, jet information, trigger prescale and information about whether the trigger is defined or not. For the ISOTRK category we combine the three triggers using a novel method, which allows the combination of any number of triggers in order to maximize the event yield while avoiding trigger correlations. On an event-by-event basis, only the trigger with the largest efficiency is used. By avoiding a logical 'OR' between triggers, the loss in the yield of events accepted by the trigger combination is compensated by a smaller and easier-to-compute corresponding systematic uncertainty. The addition of the ISOTRK charged lepton category to the TIGHT category produces an increase of 33% in the WH signal yield and a decrease of 15.5% to 19.0% in the median expected 95% CL cross-section upper limits across the entire studied Higgs mass interval. The improvement in analysis sensitivity is smaller than the improvement in signal yield because the ISOTRK category has a smaller signal over background ratio than the TIGHT category, due to the looser ISOTRK reconstruction criteria. The observed (median expected) 95% CL SM Higgs upper limits on cross section times branching ratio vary between 2.39 x SM (2.73 x SM) for a Higgs mass of 100 GeV/c 2 to 31.1 x SM (31.2 x SM) for a Higgs mass of 150 GeV/c 2 , while the value for a 115 GeV/c 2 Higgs boson is that of 5.08 x SM (3.79 x SM). The novel trigger combination method is already in use by several CDF analyses. It is applicable to any analysis that uses triggers based on MET and jets, such as supersymmetry searches at the ATLAS and CMS experiments at the Large Hadron Collider. In its most general form, the method can be used by any analysis that combines any number of different triggers.

43 PARTICLE ACCELERATORS↗
NASA NTRS1989

Multi-level bandwidth efficient block modulation codes

The multilevel technique is investigated for combining block coding and modulation. There are four parts. In the first part, a formulation is presented for signal sets on which modulation codes are to be constructed. Distance measures on a signal set are defined and their properties are developed. In the second part, a general formulation is presented for multilevel modulation codes in terms of component codes with appropriate Euclidean distances. The distance properties, Euclidean weight distribution and linear structure of multilevel modulation codes are investigated. In the third part, several specific methods for constructing multilevel block modulation codes with interdependency among component codes are proposed. Given a multilevel block modulation code C with no interdependency among the binary component codes, the proposed methods give a multilevel block modulation code C which has the same rate as C, a minimum squared Euclidean distance not less than that of code C, a trellis diagram with the same number of states as that of C and a smaller number of nearest neighbor codewords than that of C. In the last part, error performance of block modulation codes is analyzed for an AWGN channel based on soft-decision maximum likelihood decoding. Error probabilities of some specific codes are evaluated based on their Euclidean weight distributions and simulation results.

Lin, Shu↗
DOE OSTI2026

Identifying Differential Equations in Fourier Domain (FourierIdent)

We investigate identifying differential equations in the frequency domain. Fourier analysis is an important tool in theoretical analysis and numerical solvers of differential equations, yet there is limited work in exploring this connection in the identification of differential equations. This paper aims to identify the underlying differential equation in the frequency domain, from a given single realization of the differential equation perturbed by noise. Such setting imposes difficulties which are different from other identification methods where computation is carried out in the physical domain. We propose several ways to mitigate the challenges arising from noise in data and large differences in the magnitudes of frequency responses. The main takeaways are that identifying differential equations solely in the frequency domain is challenging, the method we propose is based on a form of domain partitions in the frequency domain, and this method shows benefits for complex data even with high level of noise. We introduce a Fourier feature denoising, and define the meaningful data region and the core regions of features to reduce the effect of noise in the frequency domain and to enhance the accuracy in coefficient identification. The proposed method is tested on various differential equations with linear, nonlinear, and high-order derivative feature terms, and shows advantages on complex data with many frequency modes, even under high level of noise.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Using a Coarse-Grained Modeling Framework to Identify Oligomeric Motifs with Tunable Secondary Structure

Coarse-grained modeling can be used to explore general theories that are independent of specific chemical detail. In this paper, we present cg_openmm, a Python-based simulation framework for modeling coarse-grained hetero-oligomers and screening them for structural and thermodynamic characteristics of cooperative secondary structures. cg_openmm facilitates the building of coarse-grained topology and random starting configurations, setup of GPU-accelerated replica exchange molecular dynamics simulations with the OpenMM software package, and features a suite of postprocessing thermodynamic and structural analysis tools. In particular, native contact analysis, heat capacity calculations, and free energy of folding calculations are used to identify and characterize cooperative folding transitions and stable secondary structures. In this work, we demonstrate the capabilities of cg_openmm on a simple 1–1 Lennard-Jones coarse-grained model, in which each residue contains 1 backbone and 1 side-chain bead. By scanning both nonbonded and bonded force-field parameter spaces at the coarse-grained level, we identify and characterize sets of parameters which result in the formation of stable helices through cooperative folding transitions. Furthermore, we show that the geometries and stabilities of these helices can be tuned by manipulating the force-field parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Projected Commutator DIIS Method for Accelerating Hybrid Functional Electronic Structure Calculations

The commutator direct inversion of the iterative subspace (commutator DIIS or C-DIIS) method developed by Pulay is an efficient and the most widely used scheme in quantum chemistry to accelerate the convergence of self-consistent field (SCF) iterations in Hartree–Fock theory and Kohn–Sham density functional theory. The C-DIIS method requires the explicit storage of the density matrix, the Fock matrix, and the commutator matrix. Hence, the method can only be used for systems with a relatively small basis set, such as the Gaussian basis set. We develop a new method that enables the C-DIIS method to be efficiently employed in electronic structure calculations with a large basis set such as planewaves for the first time. The key ingredient is the projection of both the density matrix and the commutator matrix to an auxiliary matrix called the gauge-fixing matrix. The resulting projected commutator-DIIS method (PC-DIIS) only operates on matrices of the same dimension as that consists of Kohn–Sham orbitals. The cost of the method is comparable to that of standard charge mixing schemes used in large basis set calculations. The PC-DIIS method is gauge-invariant, which guarantees that its performance is invariant with respect to any unitary transformation of the Kohn–Sham orbitals. We demonstrate that the PC-DIIS method can be viewed as an extension of an iterative eigensolver for nonlinear problems. We use the PC-DIIS method for accelerating Kohn–Sham density functional theory calculations with hybrid exchange–correlation functionals, and demonstrate its superior performance compared to the commonly used nested two-level SCF iteration procedure. Furthermore, we demonstrate that in the context of ab initio molecular dynamics (MD) simulation with hybrid functionals one can extrapolate the gauge-fixing matrix to achieve the goal of extrapolating the entire density matrix implicitly along the MD trajectory. Here, numerical results indicate that the new method significantly reduces the number of SCF iterations per MD step, compared to the commonly used strategy of extrapolating the electron density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

First-Principles Integrated Adsorption Modeling for Selective Capture of Uranium from Seawater by Polyamidoxime Sorbent Materials

Nuclear power is a relatively carbon-free energy source that has the capacity to be utilized today in an effort to stem the tides of global warming. The growing demand for nuclear energy, however, could put significant strain on our uranium ore resources, and the mining activities utilized to extract that ore can leave behind long-term environmental damage. A potential solution to enhance the supply of uranium fuel is to recover uranium from seawater using amidoximated adsorbent fibers. This technology has been studied for decades but is currently plagued by the material’s relatively poor selectivity of uranium over its main competitor vanadium. In this work, we investigate the binding schemes between uranium, vanadium, and the amidoxime functional groups on the adsorbent surface. Using quantum chemical methods, binding strengths are approximated for a set of complexation reactions between uranium and vanadium with amidoxime functionalities. Those approximations are then coupled with a comprehensive aqueous adsorption model developed in this work to simulate the adsorption of uranium and vanadium under laboratory conditions. Experimental adsorption studies with uranium and vanadium over a wide pH range are performed, and the data collected are compared against simulation results to validate the model. It was found that coupling ab initio calculations with process level adsorption modeling provides accurate predictions of the adsorption capacity and selectivity of the sorbent materials. Furthermore, this work demonstrates that this multiscale modeling paradigm could be utilized to aid in the selection of superior ligands or ligand compositions for the selective capture of metal ions. Furthermore, this first-principles integrated modeling approach opens the door to the in silico design of next-generation adsorbents with potentially superior efficiency and selectivity for uranium over vanadium in seawater.

36 MATERIALS SCIENCE↗