Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “condensed phase systems”

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 865 records · Page 48

Simulation of a PCM Integrated Heat Pump Using Time-of-Use Utility Structure based Control Strategy for Demand Response

To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries and constructing new peaking power plants. The aggressive carbon reduction goals of the United States are accelerating the electrification of building equipment. If space heating and cooling in the US were electrified using conventional heat pump technologies, the winter electric grid peak would approximately double. Integrating TES with heat pump (HP) enables electrification of space cooling and heating devices without overtaxing the grid.This paper evaluates the benefits of a novel phase change material (PCM) integrated heat pump configuration via simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation.To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.

Li, Zhenning↗

Multidimensional coherent spectroscopy of light-driven states and their collective modes in multiband superconductors

We present a comprehensive theory of light-controlled multiband superconductivity, and apply it to predict distinctive signatures of light-driven superconducting (SC) states in terahertz multidimensional coherent spectroscopy (THz-MDCS) experiments. We first derive gauge-invariant Maxwell-Bloch equations for multiband BCS superconductors with spatial fluctuations. We consider driving electromagnetic fields determined self-consistently by Maxwell's equations. By calculating the THz-MDCS spectra measured experimentally in the clean SC limit, we identify unique signatures of finite-momentum Cooper-pairing states that live longer than the laser pulse. They are controlled by a pair of THz laser pulses with well-defined relative phase (pulse pair). The pseudospin oscillators that describe the properties of these SC states are parametrically driven by both finite-momentum Cooper pairing and by time oscillations of the order-parameter relative phase. We show that such strong parametric driving leads to drastic changes in the THz-MDCS spectral shape from the predictions of third-order nonlinear susceptibility calculations. These spectral changes strongly depend on the interband-to-intraband interaction ratio and on the collective modes of the light-driven state. For negligible interband interaction, the spectra show a transition with increasing field, from traditional pump-probe, four-wave-mixing, and third-harmonic generation peaks determined by the laser frequency to sidebands determined by the excitations of the driven system. These sidebands emerge from difference-frequency Raman processes in the nonequilibrium SC state. For interband couplings weaker than the intraband pairing, we show that the Leggett phase collective mode leads to harmonic sidebands around the traditional pump-probe peaks. Additional Higgs collective mode peaks result from light-induced inversion-symmetry breaking in a thin-film geometry. For strong interband coupling, we find a transition from a nonequilibrium finite Cooper-pair momentum state characterized by hybrid-Higgs amplitude mode peaks in THz-MDCS spectra to a driven state identified experimentally by the emergence of Floquet-type sidebands at bi-Higgs frequencies. Those dominant bi-Higgs-frequency satellites are manifestations of a new order parameter relative phase collective mode that characterizes the nonequilibrium SC state. The predicted interaction- and field-dependent transitions in the spectral profile allow us to propose THz-MDCS experiments for quantum tomography of light-driven superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Driven superconducting vortex dynamics in systems with twofold anisotropy in the presence of pinning

Abstract We examine the dynamics of superconducting vortices with twofold anisotropic interaction potentials driven over random pinning, and compare the behavior under drives applied along the hard and the soft anisotropy directions. As the driving force increases, the number of topological defects reaches a maximum near the depinning threshold, and then decreases as the vortices form one-dimensional (1D) chains. This coincides with a transition from a pinned nematic to a moving smectic aligned with the soft anisotropy direction. The system is generally more ordered when the drive is applied along the soft direction of the anisotropy. For driving along the hard direction, there is a critical value of the twofold anisotropy above which the system remains aligned with the soft direction. Hysteretic behavior appears upon cycling the driving force, with 1D vortex chains persisting during the decreasing leg below the threshold for chain formation for increasing drive. More anisotropic systems have a greater amount of structural disorder in the moving state. For lower anisotropy, the system forms a moving smectic-A state, while at higher anisotropy, a moving nematic state appears instead.

36 MATERIALS SCIENCE↗

Direct observation of twisted stacking domains in the van der Waals magnet CrI 3

Van der Waals (vdW) stacking is a powerful technique to achieve desired properties in condensed matter systems through layer-by-layer crystal engineering. A remarkable example is the control over the twist angle between artificially-stacked vdW crystals, enabling the realization of unconventional phenomena in moiré structures ranging from superconductivity to strongly correlated magnetism. Here, we report the appearance of unusual 120° twisted faults in vdW magnet CrI 3 crystals. In exfoliated samples, we observe vertical twisted domains with a thickness below 10 nm. The size and distribution of twisted domains strongly depend on the sample preparation methods, with as-synthesized unexfoliated samples showing tenfold thicker domains than exfoliated samples. Cooling induces changes in the relative populations among different twisting domains, rather than the previously assumed structural phase transition to the rhombohedral stacking. The stacking disorder induced by sample fabrication processes may explain the unresolved thickness-dependent magnetic coupling observed in CrI 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Interfacially Enhanced Superconductivity in Fe(Te,Se)/Bi 4 Te 3 Heterostructures

Realizing topological superconductivity by integrating high-transition-temperature (T C ) superconductors with topological insulators can open new paths for quantum computing applications. Here, a new approach is reported for increasing the superconducting transition temperature $T^{onset}_C$ by interfacing the unconventional superconductor Fe(Te,Se) with the topological insulator Bi–Te system in the low-Se doping regime, near where superconductivity vanishes in the bulk. The critical finding is that the $T^{onset}_C$ of Fe(Te,Se) increases from nominally non-superconducting to as high as 12.5 K when Bi 2 Te3 is replaced with the topological phase Bi 4 Te 3 . Interfacing Fe(Te,Se) with Bi4Te3 is also found to be critical for stabilizing superconductivity in monolayer films where $T^{onset}_C$ can be as high as 6 K. Measurements of the electronic and crystalline structure of the Bi 4 Te 3 layer reveal that a large electron transfer, epitaxial strain, and novel chemical reduction processes are critical factors for the enhancement of superconductivity. This novel route for enhancing T C in an important epitaxial system provides new insight on the nature of interfacial superconductivity and a platform to identify and utilize new electronic phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic properties of epitaxial La 1- x Sr x RhO 3 thin films

Here, we report on the synthesis and electronic properties of epitaxial perovskite La 1- x Sr x RhO 3 thin films. Thin films with a Sr content ranging from $\textit{x}$ = 0 to 0.5 have been grown using molecular beam epitaxy. Transport and x-ray photoemission spectroscopy data reveal an insulator-metal-insulator transition, accompanied by a $\textit{p-}$ to $\textit{n}-$ type carrier change observed in Hall measurements. Combined with theoretical calculations, we find that the addition of Sr does not directly dope carriers into the conduction band, but rather induces localized Rh 4$\textit{d}$ states within the LaRhO 3 band gap. The bandwidth of the impurity band increases with Sr content, eventually causing the valence band (VB) and the localized Rh 4$\textit{d}$ band to overlap, which explains the first insulator-to-metal transition occurring at $\textit{x}$ = 0.35. For Sr content $\textit{x}$ > 0.4, possible cation ordering results in an increase of the gap between the VB and the Rh 4$\textit{d}$ band, leading to the second metal-to-insulator transition. We map out the electronic phase diagram of the Sr-doped LaRhO 3 system and suggest strategies to engineer the electronic states in rhodate systems via delocalizing the Rh 3+ states.

36 MATERIALS SCIENCE↗

Solubility-limited depolymerization kinetics in the glycolysis of carbonyl-containing polymers

Chemical recycling of condensation polymers is often rationalized on the basis of the intrinsic reactivity of ester and carbonate functional groups. However, under heterogeneous conditions relevant to plastic waste processing and environmental degradation, bulk depolymerization rates often diverge from trends predicted by homogeneous chemistry. Here, we investigate how polymer–solvent compatibility, catalyst strength, and phase behavior govern the heterogeneous glycolysis of carbonyl-containing polymers. Using poly(ethylene terephthalate) (PET), glycol-modified PET (PETG), and bisphenol-A polycarbonate (PC) as model systems, we examine depolymerization kinetics at 180 °C with ethylene glycol and bisphenol A as diols under both amphoteric organosalt (TBD : MSA) and strong base (TBD) catalysis. Despite substantial differences in crystallinity and glycol uptake, PET and PETG depolymerize at comparable rates under organosalt catalysis, while PC depolymerizes significantly more slowly under identical conditions. Time-resolved molecular weight analysis and thermal characterization demonstrate that these rate differences do not arise from crystallinity, swelling, or inherent carbonyl reactivity, but instead reflect solubility-limited kinetics that constrain the transition from heterogeneous to homogeneous reaction regimes. When polymer solubility is low, depolymerization remains heterogeneous and slow; when solubility is enhanced—either through increased polymer–diol compatibility or stronger base catalysis—rapid homogeneous depolymerization is observed, reversing apparent reactivity trends. These results establish solubility and phase behavior as primary determinants of depolymerization kinetics in heterogeneous polymer recycling systems. By demonstrating how catalyst selection and solvent compatibility can expose or overcome solubility limitations, this work provides mechanistic insight to design more energy-efficient and selective chemical recycling processes. More broadly, these findings suggest that polymers with limited solvent or water compatibility may resist chemical degradation in the environment, favoring fragmentation and persistence as micro- and nanoplastics. Understanding solubility-controlled depolymerization offers a pathway toward more sustainable polymer design and end-of-life chemical recovery.

Watson-Sanders, Shelby [Department of Chemistry, U↗

Origin of Ca-Al-rich inclusions. II - Sputtering and collisions in the three-ph8se interstellar medium

The theory put forward by Clayton (1977) for the formation of the Ca-Al-rich inclusions within C3 meteorites is extended to an evolutionary history in a three-phase interstellar medium. Widespread supersonic turbulence in the hot interstellar medium is maintained by supernova shock waves, giving rise to heavy sputtering of the refractory dust. Subsequent reaccumulation with varying dust/gas ratios or varying particle sizes produces isotopically fractionated Ca-Al-rich accumulates. It is thought that the Ca-Al-rich inclusions themselves are formed by the following sequence in the solar system: (1) cold accumulation of larger-than-average Ca-Al-rich particles containing supernova condensate cores into macroscopic (approximately 1 cm) Ca-Al-rich agglomerates, probably by sedimentation; and (2) fusion of the supernova condensates into macroscopic minerals by exothermic chemical reactions that begin when the accumulate has been warmed, thereby releasing energy from the unequilibrated forms accumulated from the interstellar medium.

Clayton, D. D.↗

Morphology and composition of condensates on Apollo 17 orange and black glass

Lunar soil sample 74220 and core samples 74001/2 consist mainly of orange glass droplets, droplet fragments, and their crystallized equivalents. These samples are now generally accepted to be pyroclastic ejecta from early lunar volcanic eruptions. It has been known since early examination of these samples that they contain surface coatings and material rich in volatile condensible phases, including S, Zn, F, Cl, and many volatile metals. The volatiles associated with these orange and black glasses (and the Apollo 15 green glasses) may provide important clues in understanding the differentiation and volcanic history of the Moon. In addition, condensible volatiles can be mobilized and concentrated by volcanic processes. We have reviewed many of our existing photomicrographs and energy dispersive analysis (EDXA) of grain surfaces and have reexamined some of our older SEM mounts using an improved EDXA system capable of light-element detection and analysis (oxygen, nitrogen, and carbon). The results from these investigations are presented.

Mckay, David S.↗

Automation of closed environments in space for human comfort and safety

The Environmental Control and Life Support System (ECLSS) for the Space Station Freedom and future colonization of the Moon and Mars presents new challenges for present technologies. Current plans call for a crew of 8 to live in a safe, shirt-sleeve environment for 90 days without ground support. Because of these requirements, all life support systems must be self-sufficient and reliable. The ECLSS is composed of six subsystems. The temperature and humidity control (THC) subsystem maintains the cabin temperature and humidity at a comfortable level. The atmosphere control and supply (ACS) subsystem insures proper cabin pressure and partial pressures of oxygen and nitrogen. To protect the space station from fire damage, the fire detection and suppression (FDS) subsystem provides fire sensing alarms and extinguishers. The waste management (WM) subsystem compacts solid wastes for return to Earth, and collects urine for water recovery. Because it is impractical, if not impossible, to supply the station with enough fresh air and water for the duration of the space station's extended mission, these elements are recycled. The atmosphere revitalization (AR) subsystem removes CO2 and other dangerous contaminants from the air. The water recovery and management (WRM) subsystem collects and filters condensate from the cabin to replenish potable water supplies, and processes urine and other waste waters to replenish hygiene water supplies. These subsystems are not fully automated at this time. Furthermore, the control of these subsystems is not presently integrated; they are largely independent of one another. A fully integrated and automated ECLSS would increase astronauts' productivity and contribute to their safety and comfort. The Kansas State University Advanced Design Team is in the process of researching and designing controls for the automation of the ECLSS for Space Station Freedom and beyond. The approach chosen to solve this problem is to divide the design into three phases. The first phase is to research the ECLSS as a whole system and then concentrate efforts on the automation of a single subsystem. The AR subsystem was chosen for our focus. During the second phase, the system control process will then be applied to the AR subsystem.

Source record↗

Calibration of High Frequency MEMS Microphones

Understanding and controlling aircraft noise is one of the major research topics of the NASA Fundamental Aeronautics Program. One of the measurement technologies used to acquire noise data is the microphone directional array (DA). Traditional direction array hardware, consisting of commercially available condenser microphones and preamplifiers can be too expensive and their installation in hard-walled wind tunnel test sections too complicated. An emerging micro-machining technology coupled with the latest cutting edge technologies for smaller and faster systems have opened the way for development of MEMS microphones. The MEMS microphone devices are available in the market but suffer from certain important shortcomings. Based on early experiments with array prototypes, it has been found that both the bandwidth and the sound pressure level dynamic range of the microphones should be increased significantly to improve the performance and flexibility of the overall array. Thus, in collaboration with an outside MEMS design vendor, NASA Langley modified commercially available MEMS microphone as shown in Figure 1 to meet the new requirements. Coupled with the design of the enhanced MEMS microphones was the development of a new calibration method for simultaneously obtaining the sensitivity and phase response of the devices over their entire broadband frequency range. Over the years, several methods have been used for microphone calibration. Some of the common methods of microphone calibration are Coupler (Reciprocity, Substitution, and Simultaneous), Pistonphone, Electrostatic actuator, and Free-field calibration (Reciprocity, Substitution, and Simultaneous). Traditionally, electrostatic actuators (EA) have been used to characterize air-condenser microphones for wideband frequency ranges; however, MEMS microphones are not adaptable to the EA method due to their construction and very small diaphragm size. Hence a substitution-based, free-field method was developed to calibrate these microphones at frequencies up to 80 kHz. The technique relied on the use of a random, ultrasonic broadband centrifugal sound source located in a small anechoic chamber. Phase calibrations of the MEMS microphones were derived from cross spectral phase comparisons between the reference and test substitution microphones and an adjacent and invariant grazing-incidence 1/8-inch standard microphone.

Shams, Qamar A.↗

Compositional phase stability of correlated electron materials within DFT+DMFT

Predicting the compositional phase stability of strongly correlated electron materials is an outstanding challenge in condensed matter physics, requiring precise computations of total energies. In this work, we employ the density functional theory plus dynamical mean-field theory (DFT+DMFT) formalism to address local correlations due to transition metal $\textit{d}$ electrons on compositional phase stability in the prototype rechargeable battery cathode material Li x CoO 2 , and detailed comparisons are made with the simpler DFT+$\textit{U}$ approach (i.e., the Hartree-Fock solution of the DMFT impurity problem). Local interactions are found to strongly impact the energetics of the band insulator LiCoO 2 , most significantly via the $E_g$ orbitals, which are partially occupied via hybridization with $\textit{O p}$ states. We find CoO 2 and Li 1/2 CoO 2 to be moderately correlated Fermi liquids with quasiparticle weights of 0.6–0.8 for the $T_{2g}$ states, which are most impacted by the interactions. Further, as compared to DFT+$\textit{U}$, DFT+DMFT considerably dampens the increase in total energy as $\textit{U}$ is increased, which indicates that dynamical correlations are important to describe this class of materials despite the relatively modest quasiparticle weights. Unlike DFT+$\textit{U}$, which can incorrectly drive Li x CoO 2 toward spurious phase separating or charge-ordered states, DFT+DMFT correctly captures the system's phase stability and does not exhibit a strong charge-ordering tendency. Most importantly, the error within DFT+$\textit{U}$ varies strongly as the composition changes, challenging the common practice of artificially tuning $\textit{U}$ within DFT+$\textit{U}$ to compensate the errors of Hartree-Fock. DFT+DMFT predicts the average intercalation voltage decreases relative to DFT, opposite to the result of DFT+$\textit{U}$, which would yield favorable agreement with experiment in conjunction with the overprediction of the voltage by the strongly constrained and appropriately normed DFT functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hybrid Symmetry Epitaxy of the Superconducting Fe(Te,Se) Film on a Topological Insulator

It is challenging to grow an epitaxial 4-fold compound superconductor (SC) on a 6-fold topological insulator (TI) platform due to the stringent lattice-matching requirement. Here, we demonstrate that Fe(Te,Se) can grow epitaxially on a TI (Bi 2 Te 3 ) layer due to accidental, uniaxial lattice match, which is dubbed as “hybrid symmetry epitaxy”. This new growth mode is critical to stabilizing robust superconductivity with TC as high as 13 K. Furthermore, the superconductivity in this FeTe 1–x Se x /Bi 2 Te 3 system survives in the Te-rich phase with Se content as low as x = 0.03 but vanishes at Se content above x = 0.56, exhibiting a phase diagram that is quite different from that of the conventional Fe(Te,Se) systems. This unique heterostructure platform that can be formed in both TI-on-SC and SC-on-TI sequences opens a route to unprecedented topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Two-component nematic superconductivity in 4Hb-TaS 2

Most superconductors have an isotropic, single component order parameter and are well described by the standard (BCS) theory for superconductivity. Unconventional, multiple-component superconductors are exceptionally rare and are much less understood. Here, we combine scanning tunneling microscopy and angle-resolved macroscopic transport for studying the candidate chiral superconductor, 4Hb-TaS 2 . We reveal quasi-periodic one-dimensional modulations in the tunneling conductance accompanied by two-fold symmetric superconducting critical field. The strong modulation of the in-plane critical field, H c2 , points to a nematic, unconventional order parameter. However, the imaged vortex core is isotropic at low temperatures. We suggest a model that reconciles this apparent discrepancy and takes into account previously observed spontaneous time-reversal symmetry breaking at low temperatures. The model describes a competition between a dominating chiral superconducting order parameter and a nematic one. The latter emerges close to the normal phase. Our results strongly support the existence of two-component superconductivity in 4Hb-TaS 2 and can provide valuable insights into other systems with coexistent charge order and superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ultrafast nonlinear absorption of Haldane model quantum dots

We study theoretically the nonlinear absorbance of Haldane model quantum dots (QDs) placed in the field of an ultrashort and strong optical pulse. The absorbance strongly depends on the frequency of the pulse. When the frequency of the pulse is much less than the QD bandgap, the absorbance shows strong dependence on the pulse amplitude and, as a function of an internal phase of the Haldane model, the absorbance has maxima at intermediate values of the phase. When the frequency of the pulse becomes closer to, but still less than, the bandgap, the absorbance has a weak dependence on the pulse amplitude and, as a function of the internal phase, it has a maximum at the phase of 900 when the QD bandgap also has the smallest value. Furthermore, nonlinear electron dynamics in such QD systems changes from almost reversible one at small pulse frequencies to highly irreversible dynamics at large frequencies of the pulse.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Finite-size subthermal regime in disordered SU ( N ) -symmetric Heisenberg chains

SU⁡(N) symmetry is incompatible with the many-body localized (MBL) phase, even when strong disorder is present. However, recent studies have shown that finite-size SU⁡(2) systems exhibit nonergodic, subthermal behavior, characterized by the breakdown of the eigenstate thermalization hypothesis, and by the excited eigenstates entanglement entropy that is intermediate between area and volume law. In this paper, we extend previous studies of the SU⁡(2)-symmetric disordered Heisenberg model to larger systems, using the time-dependent density matrix renormalization group (tDMRG) method. We simulate quench dynamics from weakly entangled initial states up to long times, finding robust subthermal behavior at stronger disorder. Although we find an increased tendency towards thermalization at larger system sizes, the subthermal regime persists at intermediate time scales, nevertheless, and therefore should be accessible experimentally. At weaker disorder, we observe signatures of thermalization; however, entanglement entropy exhibits slow sublinear growth, in contrast to conventional thermalizing systems. Furthermore, we study dynamics of the SU⁡(3)-symmetric disordered Heisenberg model. Similarly, strong disorder drives the system into subthermal regime, albeit thermalizing phase is broader compared to the SU⁡(2) case. Finally, our findings demonstrate the robustness of the subthermal regime in spin chains with non-Abelian continuous symmetry, and are consistent with eventual thermalization at large system sizes and long time scales, suggested by previous studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Controlling topology through targeted composite symmetry manipulation in magnetic systems

The possibility of selecting magnetic space groups by orienting the magnetization direction or tuning magnetic orders offers a vast playground for engineering symmetry-protected topological phases in magnetic materials. In this work, we study how selective tuning of symmetry and magnetism can influence and control the resulting topology in a two-dimensional magnetic system, and we illustrate such a procedure in the ferromagnetic monolayer MnPSe 3 . Density functional theory calculations reveal a symmetry-protected accidental semimetallic (SM) phase for out-of-plane magnetization, which becomes an insulator when the magnetization is tilted in-plane, reaching band-gap values close to 100 meV. We identify an order-2 composite antiunitary symmetry and threefold rotational symmetry that induce the band crossing, and we classify the possible topological phases using symmetry analysis, which we support with tight-binding and k · p models. Breaking of inversion symmetry opens a gap in the SM phase, giving rise to a Chern insulator. We demonstrate this explicitly in the isostructural Janus compound Mn 2 ⁢P 2 ⁢S 3 ⁢Se 3 , which naturally exhibits Rashba spin-orbit coupling that breaks inversion symmetry. Our results map out the phase space of topological properties of ferromagnetic transition-metal phosphorus trichalcogenides, and they demonstrate the potential of the magnetization-dependent metal-to-insulator transition as a spin switch in integrated two-dimensional electronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Benchmarks and results of the two-band Hubbard model from the Gutzwiller conjugate gradient minimization theory

Ground-state properties, such as energies and double occupancies, of a one-dimensional two-band Hubbard model are calculated using a first-principles Gutzwiller conjugate gradient minimization theory. The favorable agreement with the results from the density matrix renormalization group theory demonstrates the accuracy of our method. A rotationally invariant approach is further incorporated into the method to greatly reduce the computational complexity with a speedup of approximately 50 times. Moreover, we investigate the Mott transition between a metal and a Mott insulator by evaluating the charge gap. In conclusion, with greatly reduced computational effort, our method reproduces the phase diagram in reasonable agreement with the density matrix renormalization group theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗