Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix completion”

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 829 records · Page 46

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Lightweight Chambers for Thrust Cell Applications

Traditional metals like steel and copper alloys have been used for many years to fabricate injector and chamber components of thruster assemblies. While the materials perform well, reducing engine weights would help existing and future vehicles gain performance and payload capability. It may now be possible to reduce current thruster weights up to 50% by applying composite materials. In this task, these materials are being applied to an existing thrust cell design to demonstrate new fabrication processes and potential weight savings. Two ceramic matrix composite (CMC) designs, three polymer matrix composite (PMC) designs, and two metal matrix composite (MMC) designs are being fabricated as small chamber demonstration units. In addition, a new alloy of copper, chrome, and niobium (Cu-8Cr-4Nb) is being investigated for thrust chamber liners since it offers higher strength and increased cycle life over traditional alloys. This new alloy is being used for the liner in each MMC and PMC demonstration unit. During June-August of 2000, hot-fire testing of each unit is planned to validate designs in an oxygen/hydrogen environment at chamber pressures around 850 psi. Although the weight savings using CMC materials is expected to be high, they have proven to be much harder to incorporate into chamber designs based on current fabrication efforts. However, the PMC & MMC concepts using the Cu-8Cr-4Nb liner are nearly complete and ready for testing. Additional efforts intend to use the PMC & MMC materials to fabricate a full size thrust chamber (60K lb(sub f) thrust class). The fabrication of this full size unit is expected to be complete by October 2000, followed by hot-fire testing in November-December 2000.

Elam, S.↗

Osteoclast cytosolic calcium, regulated by voltage-gated calcium channels and extracellular calcium, controls podosome assembly and bone resorption

The mechanisms of Ca2+ entry and their effects on cell function were investigated in cultured chicken osteoclasts and putative osteoclasts produced by fusion of mononuclear cell precursors. Voltage-gated Ca2+ channels (VGCC) were detected by the effects of membrane depolarization with K+, BAY K 8644, and dihydropyridine antagonists. K+ produced dose-dependent increases of cytosolic calcium ([Ca2+]i) in osteoclasts on glass coverslips. Half-maximal effects were achieved at 70 mM K+. The effects of K+ were completely inhibited by dihydropyridine derivative Ca2+ channel blocking agents. BAY K 8644 (5 X 10(-6) M), a VGCC agonist, stimulated Ca2+ entry which was inhibited by nicardipine. VGCCs were inactivated by the attachment of osteoclasts to bone, indicating a rapid phenotypic change in Ca2+ entry mechanisms associated with adhesion of osteoclasts to their resorption substrate. Increasing extracellular Ca2+ ([Ca2+]e) induced Ca2+ release from intracellular stores and Ca2+ influx. The Ca2+ release was blocked by dantrolene (10(-5) M), and the influx by La3+. The effects of [Ca2+]e on [Ca2+]i suggests the presence of a Ca2+ receptor on the osteoclast cell membrane that could be coupled to mechanisms regulating cell function. Expression of the [Ca2+]e effect on [Ca2+]i was similar in the presence or absence of bone matrix substrate. Each of the mechanisms producing increases in [Ca2+]i, (membrane depolarization, BAY K 8644, and [Ca2+]e) reduced expression of the osteoclast-specific adhesion structure, the podosome. The decrease in podosome expression was mirrored by a 50% decrease in bone resorptive activity. Thus, stimulated increases of osteoclast [Ca2+]i lead to cytoskeletal changes affecting cell adhesion and decreasing bone resorptive activity.

NASA Discipline Musculoskeletal↗

Modified Coulomb-Dipole Theory for 2e Photoionization

In the light of recent experiment on 2e photoionization of Li near threshold, we have considered a modification of the Coulomb-dipole theory, retaining the basic assumption that the threshold is dominated by asymmetric events in phase space [implies r(sub 1), k(sub 1)) greater than or equal to 2(r(sub 2), k(sub )]. In this region [in a collinear model, 2/r(sub 12) approached + 2/(r(sub 1)+r(sub 2)] the interaction reduces to V(rIsub 1) is greater than or equal to 2r(sub 2) is identically equal to [(-Z/r(sub 2)-(A-1)/r(sub 1)] + [(-2r(sub 2)/r(sub 1 exp 2)] is identically equal to V(sub c)+[V(sub d)]. For two electron emission Z = 2, thus both electrons see a Coulomb potential (V(sub c)) asymptotically, albeit each seeing a different charge. The residual potential (V(sub d)) is dipole in character. Writing the total psi = psi (sub c) + psi(sub d) = delta psi, and noting that. (T+V(sub c)-E)psy(sub c) = 0 and (T+V(sub c))psi(sub d) = 0 can be solved exactly, we find, substituting psi into the complete Schrod. Eq., that delta psi = -(H-E)(exp -1)(V(sub d) psi(sub 0)+V(sub c psi (sub 1). Using the fact that the absolute value of V(sub c) is much more than the absolute value of V(sub d) in almost all of configuration space, we can replace H by H(sub 0) in 9H-E)(exp -1) to obtain an improved approximation psi (improved) = psi(sub c) + psi(sub d) -(H(sub 0)-E)(exp -1) (V(sub c) psi (sub 0) + V(sub c) psi(sub 1). Here's the Green's function (H(sub 0)-E)(exp -1), can be exhibited explicitly, but the last term in psi (improved) is small, compared to the first two terms. Inserting them into the transition matrix element, which one handles in the usual way, we obtain in the limit E approaches 0, the threshold law: Q(E) alpha E + M(E)E(exp 5/4) + higher order (Eq. 1a). The modulation function, M(E), is a well-defined (but very non-trivial integral, but it is expected to be well approximated by a sinusoidal function containing a dipole phase term (M(E) = c sin[alpha log (E) + micron] (Eq. 1b). Experimental results show definite modulations, and are well fitted by Eqs (1).

Source record↗

Time-temperature-stress capabilities of composite materials for advanced supersonic technology application, phase 1

Implementation of metal and resin matrix composites into supersonic vehicle usage is contingent upon accelerating the demonstration of service capacity and design technology. Because of the added material complexity and lack of extensive service data, laboratory replication of the flight service will provide the most rapid method of documenting the airworthiness of advanced composite systems. A program in progress to determine the time temperature stress capabilities of several high temperature composite materials includes thermal aging, environmental aging, fatigue, creep, fracture, and tensile tests as well as real time flight simulation exposure. The program has two parts. The first includes all the material property determinations and aging and simulation exposures up through 10,000 hours. The second continues these tests up to 50,000 cumulative hours. Results are presented of the 10,000 hour phase, which has now been completed.

Kerr, J. R.↗

Meteor ablation spheres from deep-sea sediments

Spheres from mid-Pacific abyssal clays (0 to 500,000 yrs old), formed from particles that completely melted and subsequently recrystallized as they separated from their meteoroid bodies, or containing relict grains of parent meteoroids that did not experience any melting were analyzed. The spheres were readily divided into three groups using their dominant mineralogy. The Fe-rich spheres were produced during ablation of Fe and metal-rich silicate meteoroids. The glassy spheres are considerably more Fe-rich than the silicate spheres. They consist of magnetite and an Fe glass which is relatively low in Si. Bulk compositions and relict grains are useful for determining the parent meteoroid types for the silicate spheres. Bulk analyses of recrystallized spheres show that nonvolatile elemental abundances are similar to chondrite abundances. Analysis of relict grains identified high temperature minerals associated with a fine-grained, low temperature, volatile-rich matrix. The obvious candidates for parent meteoroids of this type of silicate sphere is a carbonaceous chondrite.

Blanchard, M. B.↗

Investigation for improving Global Positioning System (GPS) orbits using a discrete sequential estimator and stochastic models of selected physical processes

GEODYNII is a conventional batch least-squares differential corrector computer program with deterministic models of the physical environment. Conventional algorithms were used to process differenced phase and pseudorange data to determine eight-day Global Positioning system (GPS) orbits with several meter accuracy. However, random physical processes drive the errors whose magnitudes prevent improving the GPS orbit accuracy. To improve the orbit accuracy, these random processes should be modeled stochastically. The conventional batch least-squares algorithm cannot accommodate stochastic models, only a stochastic estimation algorithm is suitable, such as a sequential filter/smoother. Also, GEODYNII cannot currently model the correlation among data values. Differenced pseudorange, and especially differenced phase, are precise data types that can be used to improve the GPS orbit precision. To overcome these limitations and improve the accuracy of GPS orbits computed using GEODYNII, we proposed to develop a sequential stochastic filter/smoother processor by using GEODYNII as a type of trajectory preprocessor. Our proposed processor is now completed. It contains a correlated double difference range processing capability, first order Gauss Markov models for the solar radiation pressure scale coefficient and y-bias acceleration, and a random walk model for the tropospheric refraction correction. The development approach was to interface the standard GEODYNII output files (measurement partials and variationals) with software modules containing the stochastic estimator, the stochastic models, and a double differenced phase range processing routine. Thus, no modifications to the original GEODYNII software were required. A schematic of the development is shown. The observational data are edited in the preprocessor and the data are passed to GEODYNII as one of its standard data types. A reference orbit is determined using GEODYNII as a batch least-squares processor and the GEODYNII measurement partial (FTN90) and variational (FTN80, V-matrix) files are generated. These two files along with a control statement file and a satellite identification and mass file are passed to the filter/smoother to estimate time-varying parameter states at each epoch, improved satellite initial elements, and improved estimates of constant parameters.

Goad, Clyde C.↗

InSitu SEM Investigation of Microstructural Damage Evolution and Strain Relaxation in a Melt Infiltrated SiC/SiC Composite

With CMC components poised to complete flight certification in turbine engines on commercial aircraft within the near future, there are many efforts within the aerospace community to model the mechanical and environmental degradation of CMCs. Direct observations of damage evolution are needed to support these modeling efforts and provide quantitative measures of damage parameters used in the various models. This study was performed to characterize the damage evolution during tensile loading of a melt infiltrated (MI) silicon carbide reinforced silicon carbide (SiC/SiC) composite. A SiC/SiC tensile coupon was loaded to a maximum global stress of 30 ksi in a tensile fixture within an SEM while observations were made at 5 ksi increments. Both traditional image analysis and DIC (digital image correlation) were used to quantify damage evolution. With the DIC analysis, microscale damage was observed at the fiber-matrix interfaces at stresses as low as 5 ksi. First matrix cracking took place between 20 and 25 ksi, accompanied by an observable relaxation in strain near matrix cracks. Matrix crack opening measurements at the maximum load ranged from 200 nm to 1.5 m. Crack opening along the fiber-matrix interface was also characterized as a function of load and angular position relative to the loading axis. This characterization was funded by NASA GRC and was performed to support NASA GRC modeling of SiC/SiC environmental degradation

testing↗

Helium causing disappearance of a/2<111> dislocation loops in binary Fe-Cr ferritic alloys

In this work, single and dual-beam self-ion irradiations were performed at 500°C on ultra-high purity Fe14%Cr alloy to ~0.33 displacements-per-atom (dpa), and 0 or 3030 atomic-parts-per-million (appm) helium/dpa, respectively. Using transmission electron microscopy, we reveal that helium can drastically modify the dislocation loop Burgers vector in Fe-Cr alloys. Helium co-implantation caused complete disappearance of a/2<111> type dislocation loops, and the microstructure consisted of only a<100> loops. Conversely, a/2<111> type loops were predominant without He co-implantation. The total loop density remained largely unaffected. The results strikingly contrast literature asserting that helium stabilizes a/2<111> type loops in bcc Fe alloys, based on low temperature irradiations. Collectively analyzing the results with literature suggest that the small positive interaction between helium and self-interstitial atoms (SIA) in Fe predicted by atomistic simulations maybe insufficient to holistically explain the dislocation loop microstructure development in presence of helium. Helium-SIA positive binding inadvertently implies elevated a/2<111> loop fraction and higher loop densities that the present results contradict. Helium induced high cavity density causing a preferential loss of highly glissile <111> clusters, leaving the matrix saturated with <100> type clusters is proposed as a potential mechanism. Further, the in-situ irradiations combined with Burgers vector analysis strengthened the evidence of Cr-induced dislocation loop mobility reduction that appears to stabilize the a/2<111> type loops and causes higher loop densities in Fe-Cr alloys.

36 MATERIALS SCIENCE↗

Seismic Study of the Solar Interior: Inferences from SOI/MDI Observations During Solar Activity

We have continued in collaboration with Dr. Eff-Darwich (University of La Laguna, Tenerife, Spain) the study of the structure, asphericity and dynamics of the solar interior from p-mode frequencies and frequency splittings. In March 2001, Dr. Eff-Darwich came for 3 weeks visit to CfA. During this visit we completed our work on the inversion of the internal solar rotation rate, and submitted a paper describing this work to the Astrophysical Journal. This paper has been recently revised in response to the referee comments and I expect that it will be accepted for publication very soon. We also have analyzed helioseismic data looking for temporal variations of the solar stratification near the base of the convection zone. We have expanded on the initial work that was presented at the SOHO-10/GONG-2000 meeting (October 2000, Tenerife), and are in the process of writing this up. Substantial progress towards the characterization of high-degree p-modes has been achieved. Indeed, in collaboration Dr. Rabello-Soares (Stanford University), we have gained a clear conceptual understanding of the various elements that affect the leakage matrix of the SOI/MDI instrument. This was presented in an invited talk at the SOHO-10/GONG-2000 meeting (October 2000, Tenerife). Once we will have successfully migrated from a qualitative to a quantitative assessment of these effects, we should be able to generate high-degree p-modes frequencies so crucial in the diagnostic of the layers just below solar surface.

Korzennik, Sylvain G.↗

Friction and wear of nuclear graphite exposed to molten FLiNaK salt environment

This report formally documents the completion of the Advanced Reactor Technologies Level 3 Milestone (M3TG-25OR0501104), “Initiate graphite wear studies of graphite samples exposed to molten salt environment,” due August 1, 2025. The report summarizes the ongoing activities aimed at characterizing the friction and wear behavior of graphite in molten LiF–NaF–KF (FLiNaK) salt environment. The wear and friction behavior of self-mated ET-10 nuclear graphite were tested at varying temperatures (550°C and 650°C), contact loads (40 and 80 N) and sliding speeds (1 and 10 mm/s) in a controlled argon environment. The results were compared to the wear and friction behavior of the same graphite material in a dry argon environment. This report presents initial studies on the friction and wear behavior of graphite in sliding contact mode, with future work focused on investigating the graphite matrix materials, as well as rolling and impact contact modes. The outcomes of this report could provide insights into assessing the integrity of the graphitic components of pebble-bed fluoride salt–cooled high-temperature reactors (PB-FHRs) to enhance their safe operation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Designing the Ares I Crew Launch Vehicle Upper Stage Element and Integrating the Stack at NASA's Marshall Space Flight Center

Fielding an integrated launch vehicle system entails many challenges, not the least of which is the fact that it has been over 30 years since the United States has developed a human-rated vehicle - the venerable Space Shuttle. Over time, whole generations of rocket scientists have passed through the aerospace community without the opportunity to perform such exacting, demanding, and rewarding work. However, with almost 50 years of experience leading the design, development, and end-to-end systems engineering and integration of complex launch vehicles, NASA's Marshall Space Flight Center offers the in-house talent - both junior- and senior-level personnel - to shape a new national asset to meet the requirements for safe, reliable, and affordable space exploration solutions.' These personnel are housed primarily in Marshall's Engineering Directorate and are matrixed into the programs and projects that reside at the rocket center. Fortunately, many Apollo era and Shuttle engineers, as well as those who gained valuable hands-on experience in the 1990s by conducting technology demonstrator projects such as the Delta-Clipper Experimental Advanced, X-33, X-34, and X-37, as well as the short-lived Orbital Space Plane, work closely with industry partners to advance the nation's strategic capability for human access to space. Currently, only three spacefaring nations have this distinction, including the United States, Russia, and, more recently, China. The U.S. National Space Policy of2006 directs that NASA provide the means to travel to space, and the NASA Appropriations Act of2005 provided the initial funding to begin in earnest to replace the Shuttle after the International Space Station construction is complete in 20 IO? These and other strategic goals and objectives are documented in NASA's 2006 Strategic Plan.3 In 2005, a team of NASA aerospace experts conducted the Exploration Systems Architecture Study, which recommended a two-vehicle approach to America's next space transportation system for missions to the International Space Station in the next decade and to explore the Moon and establish an outpost around the 2020 timeframe.4 Based on this extensive study, NASA selected the Ares I crew launch vehicle configuration and the heavy-lift Ares V cargo launch vehicle (fig 1). This paper will give an overview of NASA's approach to integrating the Ares I vehicle stack using capabilities and assets that are resident in Marshall's Engineering Directorate, working in partnership with other NASA Centers and the U.S. aerospace industry. It also will provide top-level details on the progress of the in-house design of the Ares I vehicle's upper stage element.

Lyles, Garry↗

Summary of Model-based Attitude Control of LUVOIR in Modular Dynamic Analysis (MDA) Simulink Environment

This document overviews the work I completed in the second half of my summer 2018 internship experience in Code 591 at NASA Goddard Space Flight Center. Please see the former memorandum for additional context on the project. The take away point is that LUVOIR A has been modeled as three rigid bodies linked by 1-DOF rotary joints. An LQR attitude controller was designed for precision pointing of the spacecraft with the design requirement that steady state oscillations have amplitude less than a miliarcsecond. The performance of the algorithm was tested in the Modular Dynamic Analysis Simulink library developed by J. Roger Chen at NASA Goddard. While the initial simulations were of rigid bodies, subsequent simulations included flexible body modes on all three bodies of the model. The simulated structural deformations initially destabilized the LQR controller thus requiring a redesign of the K gain matrix. The final simulation demonstrated a stabilizing feedback gain with miliarcsecond precision within 400 minutes. This run used 20 modes on the spacecraft, 9 on the bus, and 100 on the payload. Future recommended work includes the development of a Vibration Isolation and Precision Pointing System (VIPPS) Simulink model. Armed with this model, the system should be modeled as four bodies whose associated flexible files must be generated from the spacecraft through Gimbal 1, Gimbal 1 through Gimbal 2, Gimbal 2 through the VIPPS interface, and the VIPPS interface through the payload.

William Bentz↗

Recrystallization, cracking, and erosion of dispersoid-strengthened tungsten materials during exposure to divertor plasmas

In this study, we investigated the effects of combined intense particle and heat flux exposure on advanced tungsten plasma-facing materials within the DIII-D fusion facility. Our test matrix included two types of dispersoid-strengthened tungsten (containing either 100 nm diameter TiO 2 or Ni particles), along with high-purity polycrystalline tungsten as a reference. This experiment relied on a sample geometry angled at 15° relative to the divertor surface, thereby allowing the surfaces to intercept steady-state perpendicular heat fluxes (q ⟂ ) ranging from 10.1 to 19.6 MW/m 2 . During each shot, the samples were exposed to 42 Hz edge-localized modes (ELMs), allowing us to test the material response to transient heating. We correlated the exposure conditions with extensive post-test surface composition analysis and microscopy to determine how the plasma modified each surface. The angled specimens closest to the strike point received the highest combined heat and particle flux and melted midway through the experiment. EBSD analysis revealed they were completely recrystallized throughout, with an average grain size >100 µm. On the other hand, the specimens that received a lower steady state heat flux survived with more superficial surface damage. Whereas the high-purity polycrystalline tungsten exhibited a higher surface roughness, the dispersoid-strengthened material exhibited more extensive shallow inter-granular cracking. In addition, the surface was depleted of dispersoids following plasma exposure, possibly because of evaporation and/or sputtering. The results described here provide insights into the performance of these materials in a fusion environment which can guide further optimization for use in long-pulse devices.

Kolasinski, Robert D. [Sandia National Laboratorie↗

Flexible Composite-Material Pressure Vessel

A proposed lightweight pressure vessel would be made of a composite of high-tenacity continuous fibers and a flexible matrix material. The flexibility of this pressure vessel would render it (1) compactly stowable for transport and (2) more able to withstand impacts, relative to lightweight pressure vessels made of rigid composite materials. The vessel would be designed as a structural shell wherein the fibers would be predominantly bias-oriented, the orientations being optimized to make the fibers bear the tensile loads in the structure. Such efficient use of tension-bearing fibers would minimize or eliminate the need for stitching and fill (weft) fibers for strength. The vessel could be fabricated by techniques adapted from filament winding of prior composite-material vessels, perhaps in conjunction with the use of dry film adhesives. In addition to the high-bias main-body substructure described above, the vessel would include a low-bias end substructure to complete coverage and react peak loads. Axial elements would be overlaid to contain damage and to control fiber orientation around side openings. Fiber ring structures would be used as interfaces for connection to ancillary hardware.

Brown, Glen↗

Tritium generation, release, and retention from in-core fluoride salt irradiations

Further understanding of tritium transport mechanisms in the combined molten fluoride salt and graphite environment is necessary for the design and licensing of a Fluoride-Salt-Cooled High-Temperature Reactor (FHR). The three in-core fluoride salt irradiations completed at the Massachusetts Institute of Technology Reactor (MITR) are a useful parallel for studying transport phenomena expected in a FHR environment. During the irradiations, evolution of tritium from the flibe salt was monitored and compared to the calculated total generation rate. A difference of 22 ± 10% between the integrated calculated tritium generation rate and the total release was measured for the third MITR irradiation (FS-3). The fraction of tritium which was not released from the salt could be explained by tritium retention in graphite. Additionally, for post irradiation examination, a thermal desorption furnace was used to heat nuclear graphite samples in order to release and measure retained tritium. The desorption analysis in this work utilized seven subsections of graphite from the second salt irradiation (FS-2); three from a disc of IG-110U and four from ARB matrix graphite. Observed desorption versus temperature as well as total tritium content in the samples after irradiation indicate that the graphites were not volumetrically saturated with tritium, but rather tritium retention was likely limited to the near-surface region. Measurements of the samples resulted in 2.90 ± 0.29 μCi/mm 2 of tritium retained by IG-110U and 1.83 ± 0.31 μCi/mm 2 for ARB during the 300 h FS-2 in-core irradiation. Based on the desorption measurements, the estimated total tritium retention in graphite from the FS-2 samples is consistent with the tritium release measurements from the FS-3 experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The accuracies of effective interactions in downfolding coupled-cluster approaches for small-dimensionality active spaces

Here, this paper evaluates the accuracy of the Hermitian form of the downfolding procedure using the double unitary coupled cluster (DUCC) ansatz on the benchmark systems of linear chains of hydrogen atoms, H6 and H8. The computational infrastructure employs the occupation-number-representation codes to construct the matrix representation of arbitrary second-quantized operators, allowing for the exact representation of exponentials of various operators. The tests demonstrate that external amplitudes from standard single-reference coupled cluster methods that sufficiently describe external (out-of-active-space) correlations reliably parameterize the Hermitian downfolded effective Hamiltonians in the DUCC formalism. The results show that this approach can overcome the problems associated with losing the variational character of corresponding energies in the corresponding SR-CC theories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Echo mapping of active galactic nuclei broad-line regions: Fundamental algorithms

We formulate and test a series of algorithms for echo mapping the emission-line regions near active galactic nuclei from measurements of correlated variability in their line and continuum light curves. The linear regularization method (LRM) employs a direct inversion of evenly spaced light-curve data, with a regularization parameter that can be used to control the trade-off between noise and resolution. Matrix formulas express the formal solution as well as its variance and covariance in terms of uncertainties in the measurements. Unlike the maximum-entropy method (MEM), LRM applies to kernels with both positive and negative values, but the results are somewhat limited by ringing effects. A positivity constraint proves effective in controlling the ringing. MEM combines regularization and positivity in a natural way, but similar results are also found using positivity constraints with nonentropic regularization functions. Direct inversions of unevenly sampled light curves require interpolating the noisy data. In this case better results are found by solving for both the continuum light curve and kernel function in a simultaneous fit to the data. Our conclusion is that while echo mapping currently gives ambiguous results, the algorithms are not the limiting factor. Progress depends on efforts to increase the accuracy and completeness of sampling of the observed light curves.

Vio, Roberto↗