Light-weight wolf trap study Final report
Light scattering cells to measure particles in liquid solutions
SEARCH · Engineering Papers
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.
Light scattering cells to measure particles in liquid solutions
As the performance of passive atomic frequency standards improves, a new limitation is encountered due to frequency fluctuations in an ancillary local oscillator (L.O.). The effect is due to time variation in the gain of the feedback which compensates L.O. frequency fluctuations. The high performance promised by new microwave and optical trapped ion standards may be severely compromised by this effect. Researchers present an analysis of this performance limitation for the case of sequentially interrogated standards. The time dependence of the sensitivity of the interrogation process to L.O. frequency fluctuations is evaluated for single-pulse and double-pulse Ramsey RF interrogation and for amplitude modulated pulses. The effect of these various time dependencies on performance of the standard is calculated for an L.O. with frequency fluctuations showing a typical 1/f spectral density. A limiting 1/sq. root gamma dependent deviation of frequency fluctuations is calculated as a function of pulse lengths, dead time, and pulse overlap. Researchers also present conceptual and hardware-oriented solutions to this problem which achieve a much more nearly constant sensitivity to L.O. fluctuations. Solutions involve use of double-pulse interrogation; alternate interrogation of multiple traps so that the dead time of one trap can be covered by operation of the other; and the use of double-pulse interrogation for two traps, so that during the time of the RF pulses, the increasing sensitivity of one trap tends to compensate for the decreasing sensitivity of the other. A solution making use of amplified-modulated pulses is also presented which shows nominally zero time variation.
An experimental study was devised to determine the profile drag on a subsonic wind tunnel model partially covered with various compliant surface materials. The model consisted of a large section of constant thickness bounded fore and aft by symmetric airfoil fairings. A flat rigid plate, the control surface, could be exchanged for four contiguous complaint panels. The flexible media generally consisted of thin polyvinylchloride membranes stretched to various tensions over trapped fluid cavities of air, water, or polyethelene oxide solution, or over dry or flooded open-celled polyurethane foams. On the basis of very accurate direct-force cell measurements, all configurations were found to yield profile drags equal to or slightly higher than that obtained for the conventional rigid surface case. The results indicate that the potential for practical drag-reducing applications will be limited to either high-density or high-speed flows where more sensible pressure perturbations occur.
Over the last several decades the city of Los Angeles, California, has been experiencing increased temperatures resulting from the urban heat island effect. This is largely due to the expansion of developed areas which allow for the trapping of heat, posing dangerous health risks. As a solution, many organizations have turned to urban greening and tree planting initiatives to help cool vulnerable communities. NASA DEVELOP has partnered with City Plants and the City of Los Angeles, Office of Forest Management to study the role of trees in urban environments and their relation to the mitigation of local urban heat islands. This team used NASA Earth observation data spanning from 2016 to 2022, including land surface temperature and Normalized Difference Vegetation Index (NDVI) data collected from Landsat 8 Thermal Infrared Sensor (TIRS) and International Space Station (ISS) ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS), respectively. Data from the National Agriculture Imagery Program (NAIP) were also used to obtain a supervised classification of tree canopy cover. Our analysis reveals a spatial and temporal connection between temperature and vegetation, suggesting that areas with more vegetation are less likely to suffer high summertime temperatures. Results also highlight the impacts of tree planting programs, such as the Vermont Corridor planting project, which increased tree canopy cover by up to 5% in the community between 2016 and 2022. These findings support the implementation of urban greening practices and inform residents and officials about how investing in trees will help mitigate increasing heat within Los Angeles.
Ccurrent natural gas condensing furnaces generate substantial acidic condensate and SOx/NOx/CO/HC/methane emissions, causing long-term environment issues related to air, water and soil. Acidic condensate occurs when flue gases within the furnaces cool down below the dew points of acidic gases, which upon condensation, combine with water vapor to produce acidic solutions. The acidic condensate causes serious corrosion and fouling problems in natural gas condensing furnaces. This paper presents a novel solution to use an advanced adsorption technology based on monolithic acidic gas trap (AGT) adsorbers for SOx trapping, NOx redox, and formic acid/CO/HC/methane oxidation, enabling a new natural gas furnace with utra-clean flue gas and neutral condensate. When the AGT adsorber is appropriately integrated with natural gas furnaces, the acidic gases of SOx, NOx and formic gas/CO/HC will be continuously adsorbed, redoxed or oxidized, respectively, from the flue gas during a regular heating season, except for occasional regeneration activities for SOx trapping. In the paper, we present results from a prototype furnace employing the novel AGT component, which enables a pH of the condensate =7, NOx emissions of 1-2 nanograms/joule and an annual fuel utilization efficiency (AFUE) of 96%, providing a cost-effective technology for condensing furnaces that is eco-friendly to the environment.
Overview This dataset compiles dam sediment management and monitoring information from surveys, case studies, and journal articles. Additionally, features described by the National Inventory of Dams (i.e., presence of sluice gates) are included to indicate known infrastructure features that may address sediment releases. The location and description of records from downstream monitoring gages are catalogued in order to help with tracking conditions over time (e.g., before and after management actions, as operations change, etc.). The data help address national scale understanding of challenges and solutions related to the accumulation of sediment behind a dam as well as downstream passage. Sediment trapping causes problems as it reduces storage capacity, disrupts dam and reservoir function, impedes access for recreation, alters water quality/habitat conditions, and contributes to riverbank and coastal erosion within the reservoir. Data compilation from a variety of sources is a first step towards assessing system-wide efficacy of management solutions. This dataset was developed under the Water Power Technologies Office funded effort which began as a Seedling on Reservoir Sedimentation Data, and was supported by the Reservoir Sedimentation Modeling Framework and Data Analysis project. These projects have addressed challenges in describing sediment transport, trapping, and management at dams throughout the US. Methodology An outer join on dams/reservoirs with surveys and survey reports (documented in the RESSED database, USBR or USACE databases, project websites, etc.) with the National Inventory of Dams, based on the NIDID to determine dams with documented management and/or sluice gates. Additional dams with documented management activity were identified through review of technical articles from the past 25 years in Journal of Hydrology, Journal of Water Resources Planning and Management, Geomorphology, Journal of Hydraulic Engineering, Water, Journal of Cleaner Production, International Journal of Sediment Research, Nature Scientific Reports, Earth Surface Processes and Landforms, and Environmental Science and Pollution Research. Individual records were created for each survey or management activity documented. To evaluate downstream sediment monitoring records, the nhdPlusTools and dataRetrieval packages in R were used to find gages within 10km of each dam in the management database. Length of record and location of matched gages were retrieved for those parameters relevant to sediment concentration or total sediment discharge.
There is a growing demand to develop high-throughput and high-sensitivity mass spectrometry methods for single-cell proteomics. The commonly used isobaric labeling-based multiplexed single-cell proteomics approach suffers from distorted protein quantification due to co-isolated interfering ions during MS/MS fragmentation, also known as ratio compression. We reasoned that the use of MS3-based quantification could mitigate ratio compression and provide better quantification. However, previous studies indicated reduced proteome coverages in the MS3 method, likely due to long duty cycle time and ion losses during multilevel ion selection and fragmentation. Here, in this paper, we described an improved MS acquisition method for MS3-based single-cell proteomics by employing a linear ion trap to measure reporter ions. We demonstrated that linear ion trap can increase the proteome coverages for single-cell-level peptides with even higher gain obtained via the MS3 method. The optimized real-time search MS3 method was further applied to study the immune activation of single macrophages. Among a total of 126 single cells studied, over 1200 and 1000 proteins were quantifiable when at least 50 and 75% nonmissing data were required, respectively. Our evaluation also revealed several limitations of the low-resolution ion trap detector for multiplexed single-cell proteomics and suggested experimental solutions to minimize their impacts on single-cell analysis.
Organic semiconductors (OSCs) are a promising candidate to produce low-cost, large area, and flexible electronics. There has previously been successful development of organic field effect transistors (OFETs), photovoltaics (OPV), and bioelectronics using OSCs. Solution processing of these materials allows for the fabrication of large-area devices in the kinetic crystallization regime. The charge transport capabilities have been shown to depend on the morphology and molecular packing of the OSC within thin films. Our hypothesis is that solution processing conditions may significantly impact OSC morphology and as a result the charge transport. Therefore, this work has focused on obtaining a fundamental understanding of the morphology and molecular packing of OSCs for optimal device performance. Through our work, we have gained a better understanding of how these structural conditions were influenced by the solution processing conditions. Our studies have led to a more systematic understanding of the various parameters that impact OSC morphology. Solution processing leads to non-equilibrium films, thus allowing for the formation of diverse morphologies that are inaccessible by other fabrication methods. Previously, our group has focused on tuning the morphology of small-molecular OSCs [53, 55-57]. However, little work had been devoted to polymer OSCs. There was a lack of detailed studies characterizing the solution-state of polymer OSCs in terms of their conformation, degree of entanglement, polymer aggregation, and chain relaxation dynamics. The solution state properties are also likely to be influenced by the rigidity and molecular weight of the polymer OSCs investigated. Thus, we have focused our attention on gaining a better insight of polymer OSC films and solution processing methods. Through this proposal, our approach is to investigate the correlation between solution-state properties and the morphology of the resulting polymer OSC films. We worked three specific aims: investigate the effects of 1) polymer OSC solution-state properties on final film morphology; 2) molecular additives on the solution-state and final film properties; 3) controlled pre-aggregation in the solution-state in the final film morphology. More rigid and planar polymer backbones should promote interchain charge transport and more efficient interchain hopping between polymer chains. Through tailoring the polymer backbone, polymer sidechains, and molecular weight, we expected the altered solution-state properties to affect the final film morphology. In addition, reducing the entanglements and promoting chain alignment will likely prevent charge carrier trapping through conformation disorders. We thus studied different mechanisms, such as molecular additives, to reduce entanglements in solution. Devices fabricated from these solutions were expected to have improved charge transport abilities. In addition to tailoring the solution-state characteristics of the polymer OSCs, we investigated the effect of solution processing methods on the molecular packing and morphology of polymer OSC films. Two main solution processing techniques, e.g., spin-coating and solution shearing, were employed to fabricate OFETs. Spin coating was employed to prepare OSC films, which produce isotropic films. This technique creates several parameters to tune such as spin coating speed, acceleration, and time. Additionally, our research group developed the solution shearing method, which consists of the solution initially sandwiched between two plates. By sliding the top plate, the solution front is exposed and drying begins. This technique allows for the creation of aligned large crystalline domains. Solution shearing consists of different processing parameters which can affect the final morphology, such as shearing speed, substrate temperature and temperature gradient, distance between both plates, and the tilt angle of the top plate. Due to the complexity of the polymer systems investigated, numerous techniques were employed to characterize the polymer OSC solution state and films. All the materials were characterized using the DOE supported synchrotron X-ray scattering facilities at the Stanford Synchrotron Radiation Lightsource (SSRL). Using grazing incidence X-ray scattering (GIXS) and Near Edge X-ray Absorption Fine Structure (NEXAFS) techniques, the crystalline structure and molecular orientations of the thin films were measured. Optical absorption (UV-Vis) spectroscopy was employed to determine the aggregation state in solution and films of the polymer systems. Polarized UV-Vis also allowed for the determination of the relative degree of polymer chain alignment for solution sheared films. Additionally, various other techniques were used to investigate other properties within the film, such as atomic force microscopy (AFM) and solution rheology. Finally, the device performance is quantified through the fabrication and characterization of OFETs, which will highlight the effects of morphology on charge transport.
The general concepts and theories of crystal growth are proven to work for biomolecular crystallization. This allowed us to extract basic parameters controlling growth kinetics - free surface energy, alpha, and kinetic coefficient, beta, for steps. Surface energy per molecular site in thermal units, alpha(omega)(sup 2/3)/kT approx. = 1, is close to the one for inorganic crystals in solution (omega is the specific molecular volume, T is the temperature). Entropic restrictions on incorporation of biomolecules into the lattice reduce the incorporation rate, beta, by a factor of 10(exp 2) - 10(exp 3) relative to inorganic crystals. A dehydration barrier of approx. 18kcal/mol may explain approx. 10(exp -6) times difference between frequencies of adding a molecule to the lattice and Brownian attempts to do so. The latter was obtained from AFM measurements of step and kink growth rates on orthorhombic lysozyme. Protein and many inorganic crystals typically do not belong to the Kossel type, thus requiring a theory to account for inequivalent molecular positions within its unit cell. Orthorhombic lysozyme will serve as an example of how to develop such a theory. Factors deteriorating crystal quality - stress and strain, mosaicity, molecular disorder - will be reviewed with emphasis on impurities. Dimers in ferritin and lysozyme and acetylated lysozyme, are microheterogeneous i.e. nearly isomorphic impurities that are shown to be preferentially trapped by tetragonal lysozyme and ferritin crystals, respectively. The distribution coefficient, K defined as a ratio of the (impurity/protein) ratios in crystal and in solution is a measure of trapping. For acetylated lysoyzme, K = 2.15 or, 3.42 for differently acetylated forms, is independent of both the impurity and the crystallizing protein concentration. The reason is that impurity flux to the surface is constant while the growth rate rises with supersaturation. About 3 times lower dimer concentration in space grown ferritin and lysozyme crystals might be examples explaining higher quality of the space grown protein crystal. Depletion of solution with respect to isomorphic impurities around a growing crystal may be K times deeper than with respect to the crystallizing protein.
Gibbsite [α-Al(OH) 3 ] is the solubility limiting phase for aluminum across a wide pH range, and it is a common mineral phase with many industrial applications. The growth mechanism of this layered-structure material, however, remains incompletely understood. Synthesis of gibbsite at low to circumneutral pH yields nanoplates with substantial interlayer disorder. Here we examine defects in this material in detail, and the effects of recrystallization in highly alkaline sodium hydroxide solution at 80 °C. We employed a multimodal approach, including scanning electron microscopy, magic-angle spinning nuclear magnetic resonance (MAS-NMR), Raman and infrared spectroscopies, X-ray diffraction (XRD), and X-ray total scattering pair distribution function (XPDF) analysis to characterize the ageing of the nanoplates over several days. Additionally, XRD and XPDF indicate that gibbsite nanoplates precipitated at circumneutral pH contain dense, truncated sheets imparting a local difference in interlayer distance. These interlayer defects appear well described by flat Al 13 aluminum hydroxide nanoclusters nearly isostructural with gibbsite sheets present under synthesis conditions and trapped as interlayer inclusions during growth. Ageing at elevated temperature in alkaline solutions gradually improves crystallinity, showing a gradual increase in H-bonding between interlayer OH groups. Between 7 to 8 vol% of the initial gibbsite nanoparticles exhibit this defect, with the majority of differences disappearing after 2–4 hours of recrystallization in alkaline solution. The results not only identify the source of disorder in gibbsite formed under acidic/neutral conditions but also point to a possible cluster-mediated growth mechanism evident through inclusion of relict oligomers with gibbsite-like topology trapped in the interlayer spaces.
In bacterial defense and genome editing applications, the CRISPR-associated protein Cas9 searches millions of DNA base pairs to locate a 20-nucleotide, guide RNA-complementary target sequence that abuts a protospacer-adjacent motif (PAM). Target capture requires Cas9 to unwind DNA at candidate sequences using an unknown ATP-independent mechanism. Here we show that Cas9 sharply bends and undertwists DNA on PAM binding, thereby flipping DNA nucleotides out of the duplex and toward the guide RNA for sequence interrogation. Cryogenic-electron microscopy (cryo-EM) structures of Cas9–RNA–DNA complexes trapped at different states of the interrogation pathway, together with solution conformational probing, reveal that global protein rearrangement accompanies formation of an unstacked DNA hinge. Bend-induced base flipping explains how Cas9 ‘reads’ snippets of DNA to locate target sites within a vast excess of nontarget DNA, a process crucial to both bacterial antiviral immunity and genome editing. This mechanism establishes a physical solution to the problem of complementarity-guided DNA search and shows how interrogation speed and local DNA geometry may influence genome editing efficiency.
Imaging using X-ray absorption contrast by synchrotron nano-computed tomography (47.4 nm/voxel, 7keV) of a Ryugu pyrrhotite crystal revealed probable fluid inclusions. The inclusions were completely encapsuled in pyrrhotite, suggesting the fluids were trapped in the early stages of crystal growth, and filled with low-Z material. We performed Time-of Flight-Secondary Ion Mass Spectrometry (TOF-SIMS) depth profiling and high-resolution mapping at -120°C to exposure and measure the composition of the trapped fluids in a frozen state. These measurements revealed that the ancient, trapped fluids on Ryugu’s parent asteroid were saline aqueous solutions containing H2O, CO2, sulfur species, and nitrogen- and chlorine-bearing organic compounds identified by representative secondary ion species including O(-), OH(-), CO(-), S(-), Cl(-), C2(-), C2H(-), and CN(-). The inferred presence of CO2 indicates formation of the sulfides, and by implication the Ryugu parent body, beyond the H2O and CO2 snow lines of the early solar system, i.e. > 3-4 au from the Sun.
Savannah River Remediation (SRR) manages and operates the liquid waste facilities at Savannah River Site (SRS) for the Department of Energy (DOE). Stored liquid waste is a complex mixture of insoluble solids (sludge) and soluble salts in an alkaline solution. SRR has deployed the Tank Closure Cesium Removal (TCCR) system, a tank-side ion exchange process, to remove radioactive cesium from salt waste and enable onsite disposal of the resulting decontaminated salt solution as low-level waste at the Saltstone facilities. The TCCR system consists of two prefilters, four ion exchange (IX) columns, one resin trap, and a ventilation system. The IX process uses a form of inorganic crystalline silicotitanate (CST), which has a high affinity for cesium and other alkali metals, strontium, and actinides. This process is currently deployed utilizing salt feed from Tank 10, with future plans to dissolve solid salt in Tank 9 and transfer the salt solution to Tank 10 for processing through TCCR. The feed for TCCR must be created from salt-cake in Tank 10 through a dissolution process. Once enough salt has been dissolved, a qualification process is entered. This process characterizes the feed and ensures the cesium loading on the columns will not cause boiling of waste within the columns during or after processing. Once the batch has been qualified, salt waste is fed to the TCCR system through a transfer pump in the center of the tank. The waste is filtered through a set of two shielded, dead-end prefilters that prevent solids buildup in the columns. The filtered salt solution then travels to the shielded IX columns, which can be operated individually or in series, where the cesium is sorbed on the CST media. The decontaminated salt solution (DSS) then travels through a resin trap and out of the module to Tank 11. TCCR has successfully processed approximately 795,000 L of Tank 10H radioactive salt waste over two batches to date. There has not yet been a system induced shutdown. The prefilters performed as expected with only minor degradation in recovery of differential pressure after a backflush sequence. The time between backflushes decreased as each batch reached the end of processing. The hydraulics in the IXCs mostly performed as expected at all flow rates, except for one IXC that will be further investigated during Batch 3 processing. The TCCR system has shown some opportunities for more efficient processing during the length of the demonstration so far. For future processing of material from Tank 9H through Tank 10H and the TCCR unit, TCCR 1A will implement changes to the prefilters and the IXCs. The prefilters will have an increased surface area and a new filter media in an effort to increase time between filter swaps and improve backwashing cleaning capability. The IXCs will have a reduced diameter to allow for increased heat transfer out of the column and increased loading of Cs-137. Additionally, a new form of CST with an increased kinetic performance is being investigated for use during TCCR 1A operation. (authors)
A method has been developed for unambiguously measuring the exact magnetic field experienced by trapped mercury ions contained within an atomic clock intended for space applications. In general, atomic clocks are insensitive to external perturbations that would change the frequency at which the clocks operate. On a space platform, these perturbative effects can be much larger than they would be on the ground, especially in dealing with the magnetic field environment. The solution is to use a different isotope of mercury held within the same trap as the clock isotope. The magnetic field can be very accurately measured with a magnetic-field-sensitive atomic transition in the added isotope. Further, this measurement can be made simultaneously with normal clock operation, thereby not degrading clock performance. Instead of using a conventional magnetometer to measure ambient fields, which would necessarily be placed some distance away from the clock atoms, first order field-sensitive atomic transition frequency changes in the atoms themselves determine the variations in the magnetic field. As a result, all ambiguity over the exact field value experienced by the atoms is removed. Atoms used in atomic clocks always have an atomic transition (often referred to as the clock transition) that is sensitive to magnetic fields only in second order, and usually have one or more transitions that are first-order field sensitive. For operating parameters used in the (199)Hg(+) clock, the latter can be five orders of magnitude or more sensitive to field fluctuations than the clock transition, thereby providing an unambiguous probe of the magnetic field strength.
Perovskite‐based direct radiation detectors offer a compelling platform for room‐temperature gamma spectroscopy due to their high sensitivity, tunable optoelectronic properties, and compatibility with solution processing. While prior studies on formamidinium lead tribromide (FAPbBr 3 ) have shown that surface passivation mitigates deep trap state impact on the charge collection efficiency (CCE), the influence of surface‐based shallow traps on time‐resolved detector response remains largely unexplored. In this article, a digital signal processing (DSP) method is applied to analyze waveform rising edge dynamics in FAPbBr 3 single‐crystal detectors, identifying prompt and delayed charge collection components. Using the second derivative of the collected charge signal, the onset of shallow trap reemission is isolated, occurring 1.1 µs after carrier generation. Surface passivation suppressed deep trap states, shifting carrier occupancy toward shallower traps and extending the temporal window of delayed charge collection. This delayed contribution correlates with measurable changes in the temporal and amplitude distributions of detector events. By leveraging time‐domain filtering informed by these dynamics, a selective enhancement in signal quality and gamma peak‐to‐Compton ratio (PCR) is achieved, while preserving over 90% of photoelectric events. An improvement in the PCR to 1.80, from the original 1.09, is enabled by focusing on high collection efficiency events corresponding to reemitted carriers in the 1.1 to 3 µs range after the generation event. These results demonstrate a pathway for leveraging shallow trap dynamics to enhance spectral fidelity in time‐resolved perovskite detectors.
Next generation ion traps will likely need to support tens if not hundreds of ions in order to achieve several logical qubits. As we scale to those sizes, the same problems we face now – rf dissipation, control I/O, and optical access – will only grow and become more complicated. While many of these challenges can potentially be solved with technology integration, independently researching the feasibility of that integration and other solutions may help reduce the time and risk in scaling up to larger traps, by testing on smaller less complex devices. We should also consider other fabrication techniques that may help scale to larger devices, such as: through-substrate-vias (TSVs), different metal coatings, exotic rf routing, on chip laser sources, or even a secondary macroscopic trap to reload ions from. To have these technologies ready for full scale integration when we need them, ion traps with some of these capabilities need to be produced now. Developing the rigorous fabrication methods for producing reliable traps takes time and experimentation. We propose developing larger ion traps and reliable integrated technology in conjunction to make both available faster.
The acceleration of solar flare protons is considered by cyclotron damping of intense Alfven wave turbulence in a magnetic trap. The energy diffusion coefficient is computed for a near-isotropic distribution of super-Alfvenic protons and a steady-state solution for the particle spectrum is found for both transit-time and diffusive losses out of the ends of the trap. The acceleration time to a characteristic energy approximately 20 Mev/nucl can be as short as 10 sec. On the basis of phenomenological arguments an omega/2 frequency dependence for the Alfven wave spectrum is inferred. The correlation time of the turbulence lies in the range .0005 less than tau/corr less than .05s.
At least some protein crystals were found to preferentially trap microheterogeneous impurities. The latter are, for example, dimmer molecules of the crystallizing proteines (e.g. ferritin, lysozyme), or the regular molecules on which surfaces small molecules or ions are adsorbed (e.g. acetilated lysozyme) and modi@ molecular charge. Impurities may induce lattice defects and deteriorate structural resolution. Distribution of impurities between mother solution and gorwing crystal is defined by two interrelated distribution coefficients: kappa = rho(sup c, sub 2) and K = (rho(sup c, sub 2)/rho(sup c, sub 1)/rho(sub 2)/rho(sub 1). Here, rho(sub 2), rho(sub 1) and rho(sup c, sub 2) are densities of impurity (2) and regular protein (1) in solution at the growing interface and within the crystal ("c"). For the microheterogeneous impurities studied, K approx. = 2 - 4, so that kappa approx. - 10(exp 2) - 10(exp 3), since K = kappa (rho(sub 1)/rho(sup c, sub 1) and protein solubility ratio rho(sub 1)/rho(sub=p c, sub 2) much less than 1. Therefore, a crystal growing in absence of convection purifies mother solution around itself, grows cleaner and, probably, more perfect. If convection is present, the solution flow permanently brings new impurities to the crystal. This work theoretically addressed two subjects: 1) onset of convection, 2) distribution of impurities.