The impact of ice crystals on radiative forcing and remote sensing of arctic boundary-layer mixed-phase clouds

This PhD thesis is embedded into the Arctic Study of Tropospheric Aerosol, Clouds and Radiation (ASTAR) and investigates the radiative transfer through Arctic boundary-layer mixed-phase (ABM) clouds. For this purpose airborne spectral solar radiation measurements and simulations of the solar and thermal infrared radiative transfer have been performed. This work reports on measurements with the Spectral Modular Airborne Radiation measurement sysTem (SMART-Albedometer) conducted in the framework of ASTAR in April 2007 close to Svalbard. For ASTAR the SMART-Albedometer was extended to measure spectral radiance. The development and calibration of the radiance measurements are described in this work. In combination with in situ measurements of cloud particle properties provided by the Laboratoire de M¶et¶eorologie Physique (LaMP) and simultaneous airborne lidar measurements by the Alfred Wegener Institute for Polar and Marine Research (AWI) ABM clouds were sampled. The SMART-Albedometer measurements were used to retrieve the cloud thermodynamic phase by three different approaches. A comparison of these results with the in situ and lidar measurements is presented in two case studies. Beside the dominating mixed-phase clouds pure ice clouds were found in cloud gaps and at the edge of a large cloud field. Furthermore the vertical distribution of ice crystals within ABM clouds was investigated. It was found that ice crystals at cloud top are necessary to describe the observed SMART-Albedometer measurements. The impact of ice crystals on the radiative forcing of ABM clouds is in vestigated by extensive radiative transfer simulations. The solar and net radiative forcing was found to depend on the ice crystal size, shape and the mixing ratio of ice crystals and liquid water droplets.

Patterns in and kinetics of phase separation in binary mixtures

This work focused mainly on two aspects of kinetics of phase separation in binary mixtures. In the first part, we studied the interplay of hydrodynamics and the phase separation of binary mixtures. A considerably flat container (a laterally extended geometry), at an aspect ratio of 14:1 (diameter: height) was chosen, so that any hydrodynamic instabilities, if they arise, could be tracked. Two binary mixtures were studied. One was a mixture of methanol and hexane, doped with 5% ethanol, which phase separated under cooling. The second was a mixture of butoxyethanol and water, doped with 2% decane, which phase separated under heating. The dopants were added to bring down the phase transition temperature around room temperature.rnrnAlthough much work has been done already on classical hydrodynamic instabilities, not much has been done in the understanding of the coupling between phase separation and hydrodynamic instabilities. This work aimed at understanding the influence of phase separation in initiating any hydrodynamic instability, and also vice versa. Another aim was to understand the influence of the applied temperature protocol on the emergence of patterns characteristic to hydrodynamic instabilities. rnrnOn slowly cooling the system continuously, at specific cooling rates, patterns were observed in the first mixture, at the start of phase separation. They resembled the patterns observed in classical Rayleigh-Bénard instability, which arises when a liquid continuously is heated from below. To suppress this classical convection, the cooling setup was tuned such that the lower side of the sample always remained cooler by a few millikelvins, relative to the top. We found that the nature of patterns changed with different cooling rates, with stable patterns appearing for a specific cooling rate (1K/h). On the basis of the cooling protocol, we estimated a modified Rayleigh number for our system. We found that the estimated modified Rayleigh number is near the critical value for instability, for cooling rates between 0.5K/h and 1K/h. This is consistent with our experimental findings. rnrnThe origin of the patterns, in spite of the lower side being relatively colder with respect to the top, points to two possible reasons. 1) During phase separation droplets of either phases are formed, which releases a latent heat. Our microcalorimetry measurements show that the rise in temperature during the first phase separation is in the order of 10-20millikelvins, which in some cases is enough to reverse the applied temperature bias. Thus phase separation in itself initiates a hydrodynamic instability. 2) The second reason comes from the cooling protocol itself. The sample was cooled from above and below. At sufficiently high cooling rates, there are situations where the interior of the sample is relatively hotter than both top and bottom of the sample. This is sufficient to create an instability within the cell. Our experiments at higher cooling rates (5K/h and above) show complex patterns, which hints that there is enough convection even before phase separation occurs. Infact, theoretical work done by Dr.Hayase show that patterns could arise in a system without latent heat, with symmetrical cooling from top and bottom. The simulations also show that the patterns do not span the entire height of the sample cell. This is again consistent with the cell sizes measured in our experiment.rnrnThe second mixture also showed patterns at specific heating rates, when it was continuously heated inducing phase separation. In this case though, the sample was turbid for a long time until patterns appeared. A meniscus was most probably formed before the patterns emerged. We attribute the reason of patterns in this case to Marangoni convection, which is present in systems with an interface, where local differences in surface tension give rise to an instability. Our estimates for the Rayleigh number also show a significantly lower number than that's required for RB-type instability.rnrnIn the first part of the work, therefore, we identify two different kinds of hydrodynamic instabilities in two different mixtures. Both are observed during, or after the first phase separation. Our patterns compare with the classical convection patterns, but here the origins are from phase separation and the cooling protocol.rnrnIn the second part of the work, we focused on the kinetics of phase separation in a polymer solution (polystyrene and methylcyclohexane), which is cooled continuously far down into the two phase region. Oscillations in turbidity, denoting material exchange between the phases are seen. Three processes contribute to the phase separation: Nucleation of droplets, their growth and coalescence, and their subsequent sedimentation. Experiments in low molecular binary mixtures had led to models of oscillation [43] which considered sedimentation time scales much faster than the time scales of nucleation and growth. The size and shape of the sample therefore did not matter in such situations. The oscillations in turbidity were volume-dominated. The present work aimed at understanding the influence of sedimentation time scales for polymer mixtures. Three heights of the sample with same composition were studied side by side. We found that periods increased with the sample height, thus showing that sedimentation time determines the period of oscillations in the polymer solutions. We experimented with different cooling rates and different compositions of the mixture, and we found that periods are still determined by the sample height, and therefore by sedimentation time. rnrnWe also see that turbidity emerges in two ways; either from the interface, or throughout the sample. We suggest that oscillations starting from the interface are due to satellite droplets that are formed on droplet coalescence at the interface. These satellite droplets are then advected to the top of the sample, and they grow, coalesce and sediment. This type of an oscillation wouldn't require the system to pass the energy barrier required for homogenous nucleation throughout the sample. This mechanism would work best in sample where the droplets could be effectively advected throughout the sample. In our experiments, we see more interface dominated oscillations in the smaller cells and lower cooling rates, where droplet advection is favourable. In larger samples and higher cooling rates, we mostly see that the whole sample becomes turbid homogenously, which requires the system to pass the energy barrier for homogenous nucleation.rnrnOscillations, in principle, occur since the system needs to pass an energy barrier for nucleation. The height of the barrier decreases with increasing supersaturation, which in turn is from the temperature ramp applied. This gives rise to a period where the system is clear, in between the turbid periods. At certain specific cooling rates, the system can follow a path such that the start of a turbid period coincides with the vanishing of the last turbid period, thus eliminating the clear periods. This means suppressions of oscillations altogether. In fact we experimentally present a case where, at a certain cooling rate, oscillations indeed vanish. rnrnThus we find through this work that the kinetics of phase separation in polymer solution is different from that of a low molecular system; sedimentation time scales become relevant, and therefore so does the shape and size of the sample. The role of interface in initiating turbid periods also become much more prominent in this system compared to that in low molecular mixtures.rnrnIn summary, some fundamental properties in the kinetics of phase separation in binary mixtures were studied. While the first part of the work described the close interplay of the first phase separation with hydrodynamic instabilities, the second part investigated the nature and determining factors of oscillations, when the system was cooled deep into the two phase region. Both cases show how the geometry of the cell can affect the kinetics of phase separation. This study leads to further fundamental understandings of the factors contributing to the kinetics of phase separation, and to the understandings of what can be controlled and tuned in practical cases. rn

Modeling and simulations of nanoparticles in liquid crystalline systems

The aim of this work is to investigate, using extensive Monte Carlo computer simulations, composite materials consisting of liquid crystals doped with nanoparticles. These systems are currently of great interest as they offer the possibility of tuning the properties of liquid crystals used in displays and other devices as well as providing a way of obtaining regularly organized systems of nanoparticles exploiting the molecular organization of the liquid crystal medium. Surprisingly enough, there is however a lack of fundamental knowledge on the properties and phase behavior of these hybrid materials, making the route to their application an essentially empirical one. Here we wish to contribute to the much needed rationalization of these systems studying some basic effects induced by different nanoparticles on a liquid crystal host. We investigate in particular the effects of nanoparticle shape, size and polarity as well as of their affinity to the liquid crystal solvent on the stability of the system, monitoring phase transitions, order and molecular organizations. To do this we have proposed a coarse grained approach where nanoparticles are modelled as a suitably shaped (spherical, rod and disk like) collection of spherical Lennard-Jones beads, while the mesogens are represented with Gay-Berne particles. We find that the addition of apolar nanoparticles of different shape typically lowers the nematic–isotropic transition of a non-polar nematic, with the destabilization being greater for spherical nanoparticles. For polar mesogens we have studied the effect of solvent affinity of the nanoparticles showing that aggregation takes places for low solvation values. Interestingly, if the nanoparticles are polar the aggregates contribute to stabilizing the system, compensating the shape effect. We thus find the overall effects on stability to be a delicate balance of often contrasting contributions pointing to the relevance of simulations studies for understanding these complex systems.

Geological and hydrogeological features affecting migration, multi-phase partitioning and degradation of chlorinated hydrocarbons through unconsolidated porous media.

Chlorinated solvents are the most ubiquitous organic contaminants found in groundwater since the last five decades. They generally reach groundwater as Dense Non-Aqueous Phase Liquid (DNAPL). This phase can migrate through aquifers, and also through aquitards, in ways that aqueous contaminants cannot. The complex phase partitioning to which chlorinated solvent DNAPLs can undergo (i.e. to the dissolved, vapor or sorbed phase), as well as their transformations (e.g. degradation), depend on the physico-chemical properties of the contaminants themselves and on features of the hydrogeological system. The main goal of the thesis is to provide new knowledge for the future investigations of sites contaminated by DNAPLs in alluvial settings, proposing innovative investigative approaches and emphasizing some of the key issues and main criticalities of this kind of contaminants in such a setting. To achieve this goal, the hydrogeologic setting below the city of Ferrara (Po plain, northern Italy), which is affected by scattered contamination by chlorinated solvents, has been investigated at different scales (regional and site specific), both from an intrinsic (i.e. groundwater flow systems) and specific (i.e. chlorinated solvent DNAPL behavior) point of view. Detailed investigations were carried out in particular in one selected test-site, known as “Caretti site”, where high-resolution vertical profiling of different kind of data were collected by means of multilevel monitoring systems and other innovative sampling and analytical techniques. This allowed to achieve a deep geological and hydrogeological knowledge of the system and to reconstruct in detail the architecture of contaminants in relationship to the features of the hosting porous medium. The results achieved in this thesis are useful not only at local scale, e.g. employable to interpret the origin of contamination in other sites of the Ferrara area, but also at global scale, in order to address future remediation and protection actions of similar hydrogeologic settings.

Atomistic simulations of liquid crystals in the bulk and at their interfaces

In this thesis, we have dealt with several problems concerning liquid crystals (LC) phases, either in the bulk or at their interfaces, by the use of atomistic molecular dynamics (MD) simulations. We first focused our attention on simulating and characterizing the bulk smectic phase of 4-n-octyl-4'-cyanobiphenyl (8CB), allowing us to investigate the antiparallel molecular arrangement typical of SmAd smectic phases. A second topic of study was the characterization of the 8CB interface with vacuum by simulating freely suspended thin films, which allowed us to determine the influence of the interface on the orientational and positional order. Then we investigated the LC-water and LC-electrolyte water solution interface. This interface has recently found application in the development of sensors for several compounds, including biological molecules, and here we tried to understand the re-orientation mechanism of LC molecules at the interface which is behind the functioning of these sensors. The characterization of this peculiar interface has incidentally led us to develop a polarizable force field for the pentyl-cyanobiphenyl mesogen, whose process of parametrization and validation is reported here in detail. We have shown that this force field is a significant improvement over its previous, static charge non polarizable version in terms of density, orientational order parameter and translational diffusion.

Liquid crystalline elastomers as stimuli-responsive microactuators

Liquid crystalline elastomers (LCEs) are known to perform a reversible change of shape upon the phase transition from the semi-ordered liquid crystalline state to the chaotic isotropic state. This unique behavior of these “artificial muscles” arises from the self-organizing properties of liquid crystals (mesogens) in combination with the entropy-elasticity of the slightly crosslinked elastomer network. In this work, micrometer-sized LCE actuators are fabricated in a microfluidic setup. The microtubular shear flow provides for a uniform orientation of the mesogens during the crosslinking, a perquisite for obtaining actuating LCE samples. The scope of this work was to design different actuator geometries and to broaden the applicability of the microfluidic device for different types of liquid crystalline mesogens, ranging from side-chain to main-chain systems, as well as monomer and polymer precursors. For example, the thiol-ene “click” mechanism was used for the polymerization and crosslinking of main-chain LCE actuators. The main focus was, however, placed on acrylate monomers and polymers with LC side chains. A LC polymer precursor, comprising mesogenic and crosslinkable side-chains was synthesized. Used in combination with an LC monomer, the polymeric crosslinker promoted a stable LC phase, which allowed the mixture to be isothermally handled in the microfluidic reactor. If processed without the additional LC components, the polymer precursor yielded actuating fibers. A suitable co-flowing continuous phase facilitates the formation of a liquid jet and lowers the tendency for drop formation. By modification of the microfluidic device, it was further possible to prepare core-shell particles, comprised of an LCE shell and filled with an isotropic liquid. In analogy to the heart, a hollow muscle, the elastomer shell expels the inner liquid core upon its contraction. The feasibility of the core-shell particles as micropumps was demonstrated. In general, the synthesized LCE microactuators may be utilized as active components in micromechanical and lab-on-chip systems.

Liquid crystalline cellulose derivatives for mirrorless lasing

In this thesis cholesteric films made of liquid crystalline cellulose derivatives with improved optical properties were prepared. The choice of the solvent, hydrogen bond influencing additives, the synthetic realization of a very high degree of substitution on the cellulosic polymer and the use of mechanical stirring at the upper concentration limit of the liquid crystalline range were the basis for an improved alignment of the applied cellulose tricarbamates. In combination with a tuned substrate treatment and film preparation method, cholesteric films were obtained, with optical properties that were theoretically predicted and only known from low molecular weight liquid crystals so far. Subsequent polymerization allowed a permanent fixing of the alignment and the fabrication of free standing and insensitive films.rnThe incorporation of inorganic nanorods into the cholesteric host material was mediated with tailored block copolymers, available via controlled radical polymerization methods. In addition to the shape match between the rodlike mesogens of the host and the nanorods it was possible to increase the miscibility of both materials. Nevertheless, the size of the nanorods, in comparison to the mesogens, in these densely packed liquid crystalline phases as well as their long equilibration times were the reasons for phase separation. Nanorods are, in principle, valuable substitutes for organics, but their utilization in cellulosic CLC was not to be combined with a high quality alignment of the cholesteric structure.rnA swelling process of polymerized films in a dye solution or dissolving dyes in non-polymerized CLC was used for incorporation of the organic chromophores. With the first method the CLC could be aligned and polymerized without any disturbance due to dye molecules. The optical properties of dye and CLC were matched, with regard to mirrorless lasing devices. The dye was optically excited and laser emission supported by the cholesteric cavity was obtained. The polarization and wavelength of the emitted radiation as well as its bandwidth, the obtained interference pattern and threshold behavior of the emission proofed the feedback mechanism that was not believed to be realizable in liquid crystalline polymers. rnUtilization of a microfluidic co-flow injection device enabled us to transfer the properties of cellulosic CLC from the planar film shape to spherical micrometer sized particles. The pure material yielded particles with distorted mesogen alignment similar to films prepared by capillary flow. Dilution of the CLC with a solvent that migrated into the carrier phase during particle preparation provided the basis for particles with well ordered areas. rnAlthough cellulose derivatives were known for their liquid crystalline behavior for decades and synthesized in mass production, their application as feedback material was affected by bad optical properties. In comparison to low molar mass compounds, the low degree of order in the CLC phase was the cause. With the improved material, defined lasing emission was shown and characterized. Derivatives of cellulose are desirable materials, because, as a renewable resource, they are available in large amounts for a low price and need only simple derivatization reactions. The fabrication of CLC films with tunable lasing emission, for which this thesis can provide a starting point, is in good agreement with today's requirements of modern technology and its miniaturization.rn

Cerium oxide nanoparticle uptake kinetics from the gas-phase into lung cells in vitro is transport limited

Nowadays, aerosol processes are widely used for the manufacture of nanoparticles (NPs), creating an increased occupational exposure risk of workers, laboratory personnel and scientists to airborne particles. There is evidence that possible adverse effects are linked with the accumulation of NPs in target cells, pointing out the importance of understanding the kinetics of particle internalization. In this context, the uptake kinetics of representative airborne NPs over 30 min and their internalization after 24 h post-exposure were investigated by the use of a recently established exposure system. This system combines the production of aerosolized cerium oxide (CeO(2)) NPs by flame spray synthesis with its simultaneous particle deposition from the gas-phase onto A549 lung cells, cultivated at the air-liquid interface. Particle uptake was quantified by mass spectrometry after several exposure times (0, 5, 10, 20 and 30 min). Over 35% of the deposited mass was found internalized after 10 min exposure, a value that increased to 60% after 30 min exposure. Following an additional 24 h post-incubation, a time span, after which adverse biological effects were observed in previous experiments, over 80% of total CeO(2) could be detected intracellularly. On the ultrastructural level, focal cerium aggregates were present on the apical surface of A549 cells and could also be localized intracellularly in vesicular structures. The uptake behaviour of aerosolized CeO(2) is in line with observations on cerium suspensions, where particle mass transport was identified as the rate-limiting factor for NP internalization.

A Simple Visual Demonstration of Phase Diagram Concepts for an Introductory Materials Science Course Using Colored Solvents

A simple and effective demonstration to help students comprehend phase diagrams and understand phase equilibria and transformations is created using common chemical solvents available in the laboratory. Common misconceptions surrounding phase diagram operations, such as components versus phases, reversibility of phase transformations, and the lever rule are addressed. Three different binary liquid mixtures of varying compatibility create contrastive phase equilibrium cases, where colorful dyes selectively dissolved in each of corresponding phases allow for quick and unambiguous perceptions of solubility limit and phase transformations. Direct feedback and test scores from a group of students show evidence of the effectiveness of the visual and active teaching tool.

On-line SPE LC-MS/MS for the quantification of Δ9-tetrahydrocannabinol (THC) and its two major metabolites in human peripheral blood by liquid chromatography tandem mass spectrometry

A universal and robust analytical method for the determination of Δ9-tetrahydrocannabinol (THC) and two of its metabolites Δ9-(11-OH)-tetrahydrocannabinol (11-OH-THC) and 11-nor-Δ9-carboxy-tetrahydrocannabinol (THC-COOH) in human whole blood was developed and validated for use in forensic toxicology. Protein precipitation, integrated solid phase extraction and on-line enrichment followed by high-performance liquid chromatography separation and detection with a triple quadrupole mass spectrometer were combined. The linear ranges used for the three cannabinoids were from 0.5 to 20 ng/mL for THC and 11-OH-THC and from 2.5 to 100 ng/mL for THC-COOH, therefore covering the requirements for forensic use. Correlation coefficients of 0.9980 or better were achieved for all three analytes. No relevant hydrolysis was observed for THC-COOH glucuronide with this procedure--in contrast to our previous GC-MS procedure, which obviously lead to an artificial increase of the THC-COOH concentration due to the hydrolysis of the glucuronide-conjugate occurring at high pH during the phase-transfer catalyzed methylation step.

Efficient Liquid Liquid Extraction By Emulsion Formation and Separation in Robust Milli-Fluidic Devices

Conventional liquid liquid extraction (LLE) methods require large volumes of fluids to achieve the desired mass transfer of a solute, which is unsuitable for systems dealing with a low volume or high value product. An alternative to these methods is to scale down the process. Millifluidic devices share many of the benefits of microfluidic systems, including low fluid volumes, increased interfacial area-to-volume ratio, and predictability. A robust millifluidic device was created from acrylic, glass, and aluminum. The channel is lined with a hydrogel cured in the bottom half of the device channel. This hydrogel stabilizes co-current laminar flow of immiscible organic and aqueous phases. Mass transfer of the solute occurs across the interface of these contacting phases. Using a y-junction, an aqueous emulsion is created in an organic phase. The emulsion travels through a length of tubing and then enters the co-current laminar flow device, where the emulsion is broken and each phase can be collected separately. The inclusion of this emulsion formation and separation increases the contact area between the organic and aqueous phases, therefore increasing the area over which mass transfer can occur. Using this design, 95% extraction efficiency was obtained, where 100% is represented by equilibrium. By continuing to explore this LLE process, the process can be optimized and with better understanding may be more accurately modeled. This system has the potential to scale up to the industrial level and provide the efficient extraction required with low fluid volumes and a well-behaved system.

Simultaneous determination of 3-hydroxyanthranilic and cinnabarinic acid by high-performance liquid chromatography with photometric or electrochemical detection

A convenient and rapid method for the simultaneous determination by HPLC of 3-hydroxyanthranilic acid and the dimer derived by its oxidation, cinnabarinic acid, is described. Buffers or biological samples containing these two Trp metabolites were acidified to pH 2.0 and extracted with ethyl acetate with recoveries of 96.5 +/- 0.5 and 93.4 +/- 3.7% for 3-hydroxyanthranilic and cinnabarinic acid, respectively. The two compounds were separated on a reversed-phase (C18) column combined with ion-pair chromatography and detected photometrically or electrochemically. The method was applied successfully to biological systems in which formation of either 3-hydroxyanthranilic or cinnabarinic acid had been described previously. Thus, interferon-gamma-treated human peripheral blood mononuclear cells formed and released significant amounts of 3-hydroxyanthranilic acid into the culture medium and mouse liver nuclear fraction possessed high "cinnabarinic acid synthase" activity. In contrast, addition of 3-hydroxyanthranilic acid to human erythrocytes resulted in only marginal formation of cinnabarinic acid. We conclude that the method described is specific, sensitive, and suitable for the detection of the two Trp metabolites in biological systems.

Thickness measurement of dynamic thin liquid films generated by plug-annular flow in non-wetting microchannels

Surface tension forces are significant at millimeter length-scales, causing profoundly different flow morphologies in microchannels than in macroscale flows. The existence and morphology of thin liquid films is particularly relevant for predicting performance and operational stability of devices containing microscale two phase flows. Analytical, computational, and experimental methods previously employed in the study of thin liquid films are discussed. Thicknesses before and after a novel film morphology, referred to as a `shock,' are measured with a novel film thickness measurement technique that uses confocal microscopy. Film thicknesses predicted by previous work are compared to experimental results. Methods for increasing the accuracy of the confocal film thickness measurement technique are discussed.

A validated assay by liquid chromatography-tandem mass spectrometry for the simultaneous quantification of elvitegravir and rilpivirine in HIV positive patients

Because of the large variability in the pharmacokinetics of anti-HIV drugs, therapeutic drug monitoring in patients may contribute to optimize the overall efficacy and safety of antiretroviral therapy. An LC–MS/MS method for the simultaneous assay in plasma of the novel antiretroviral agents rilpivirine (RPV) and elvitegravir (EVG) has been developed to that endeavor. Plasma samples (100 μL) extraction is performed by protein precipitation with acetonitrile, and the supernatant is subsequently diluted 1:1 with 20-mM ammonium acetate/MeOH 50:50. After reverse-phase chromatography, quantification of RPV and EVG, using matrix-matched calibration samples, is performed by electrospray ionization–triple quadrupole mass spectrometry by selected reaction monitoring detection using the positive mode. The stable isotopic-labeled compounds RPV-13C6 and EVG-D6 were used as internal standards. The method was validated according to FDA recommendations, including assessment of extraction yield, matrix effects variability (<6.4%), as well as EVG and RPV short and long-term stability in plasma. Calibration curves were validated over the clinically relevant concentrations ranging from 5 to 2500 ng/ml for RPV and from 50 to 5000 ng/ml for EVG. The method is precise (inter-day CV%: 3–6.3%) and accurate (3.8–7.2%). Plasma samples were found to be stable (<15%) in all considered conditions (RT/48 h, +4°C/48 h, −20°C/3 months and 60°C/1 h). Selected metabolite profiles analysis in patients' samples revealed the presence of EVG glucuronide, that was well separated from parent EVG, allowing to exclude potential interferences through the in-source dissociation of glucuronide to parent drug. This new, rapid and robust LCMS/MS assay for the simultaneous quantification of plasma concentrations of these two major new anti-HIV drugs EVG and RPV offers an efficient analytical tool for clinical pharmacokinetics studies and routine therapeutic drug monitoring service.