Picosecond and Terahertz Perturbation of Interfacial Water and Electropermeabilization of Biological Membranes

A new horizon for neuroscience: terahertz biotechnology in brain research

Terahertz biotechnology has been increasingly applied in various biomedical fields and has especially shown great potential for application in brain sciences. In this article, we review the development of terahertz biotechnology and its applications in the field of neuropsychiatry. Available evidence indicates promising prospects for the use of terahertz spectroscopy and terahertz imaging techniques in the diagnosis of amyloid disease, cerebrovascular disease, glioma, psychiatric disease, traumatic brain injury, and myelin deficit. In vitro and animal experiments have also demonstrated the potential therapeutic value of terahertz technology in some neuropsychiatric diseases. Although the precise underlying mechanism of the interactions between terahertz electromagnetic waves and the biosystem is not yet fully understood, the research progress in this field shows great potential for biomedical noninvasive diagnostic and therapeutic applications. However, the biosafety of terahertz radiation requires further exploration regarding its two-sided efficacy in practical applications. This review demonstrates that terahertz biotechnology has the potential to be a promising method in the field of neuropsychiatry based on its unique advantages.

Progress of the Impact of Terahertz Radiation on Ion Channel Kinetics in Neuronal Cells

In neurons and myocytes, selective ion channels in the plasma membrane play a pivotal role in transducing chemical or sensory stimuli into electrical signals, underpinning neural and cardiac functionality. Recent advancements in biomedical research have increasingly spotlighted the interaction between ion channels and electromagnetic fields, especially terahertz (THz) radiation. This review synthesizes current findings on the impact of THz radiation, known for its deep penetration and non-ionizing properties, on ion channel kinetics and membrane fluid dynamics. It is organized into three parts: the biophysical effects of THz exposure on cells, the specific modulation of ion channels by THz radiation, and the potential pathophysiological consequences of THz exposure. Understanding the biophysical mechanisms underlying these effects could lead to new therapeutic strategies for diseases.

Membrane-mediated modulation of mitochondrial physiology by terahertz waves

Extensive studies have demonstrated the diverse impacts of electromagnetic waves at gigahertz and terahertz (THz) frequencies on cytoplasmic membrane properties. However, there is little evidence of these impacts on intracellular membranes, particularly mitochondrial membranes crucial for mitochondrial physiology. In this study, human neuroblast-like cells were exposed to continuous 0.1 THz radiation at an average power density of 33 mW/cm2. The analysis revealed that THz exposure significantly altered the mitochondrial ultrastructure. THz waves enhanced the enzymatic activity of the mitochondrial respiratory chain but disrupted supercomplex assembly, compromising mitochondrial respiration. Molecular dynamics simulations revealed altered rates of change in the quantity of hydrogen bonds and infiltration of water molecules in lipid bilayers containing cardiolipin, indicating the specific behavior of cardiolipin, a signature phospholipid in mitochondria, under THz exposure. These findings suggest that THz radiation can significantly alter mitochondrial membrane properties, impacting mitochondrial physiology through a mechanism related to mitochondrial membrane, and provide deeper insight into the bioeffects of THz radiation.

Recent advances and research progress on microsystems and bioeffects of terahertz neuromodulation

Terahertz waves can interact with the nervous system of organisms under certain conditions. Compared to common optical modulation methods, terahertz waves have the advantages of low photon energy and low risk; therefore, the use of terahertz waves to regulate the nervous system is a promising new method of neuromodulation. However, most of the research has focused on the use of terahertz technology for biodetection, while relatively little research has been carried out on the biological effects of terahertz radiation on the nervous system, and there are almost no review papers on this topic. In the present article, we begin by reviewing principles and objects of research regarding the biological effects of terahertz radiation and summarizing the current state of related research from a variety of aspects, including the bioeffects of terahertz radiation on neurons in vivo and in vitro, novel regulation and detection methods with terahertz radiation devices and neural microelectrode arrays, and theoretical simulations of neural information encoding and decoding. In addition, we discuss the main problems and their possible causes and give some recommendations on possible future breakthroughs. This paper will provide insight and assistance to researchers in the fields of neuroscience, terahertz technology and biomedicine.

Calcium Imaging Characterize the Neurobiological Effect of Terahertz Radiation in Zebrafish Larvae

(1) Objective: To explore the neurobiological effects of terahertz (THz) radiation on zebrafish larvae using calcium (Ca2+) imaging technology. (2) Methods: Zebrafish larvae at 7 days post fertilization (dpf) were exposed to THz radiation for 10 or 20 min; the frequency was 2.52 THz and the amplitude 50 mW/cm2. The behavioral experiments, neural Ca2+ imaging, and quantitative polymerase chain reaction (qPCR) of the dopamine-related genes were conducted following the irradiation. (3) Results: Compared with the control group, the behavioral experiments demonstrated that THz radiation significantly increased the distance travelled and speed of zebrafish larvae. In addition, the maximum acceleration and motion frequency were elevated in the 20 min radiation group. The neural Ca2+ imaging results indicated a substantial increase in zebrafish neuronal activity. qPCR experiments revealed a significant upregulation of dopamine-related genes, such as drd2b, drd4a, slc6a3 and th. (4) Conclusion: THz radiation (2.52 THz, 50 mW/cm2, 20 min) upregulated dopamine-related genes and significantly enhanced neuronal excitability, and the neurobiological effect of THz radiation can be visualized using neural Ca2+ imaging in vivo.

High Frequency Electromagnetic Radiation Stimulates Neuronal Growth and Hippocampal Synaptic Transmission

Terahertz waves lie within the rotation and oscillation energy levels of biomolecules, and can directly couple with biomolecules to excite nonlinear resonance effects, thus causing conformational or configuration changes in biomolecules. Based on this mechanism, we investigated the effect pattern of 0.138 THz radiation on the dynamic growth of neurons and synaptic transmission efficiency, while explaining the phenomenon at a more microscopic level. We found that cumulative 0.138 THz radiation not only did not cause neuronal death, but that it promoted the dynamic growth of neuronal cytosol and protrusions. Additionally, there was a cumulative effect of terahertz radiation on the promotion of neuronal growth. Furthermore, in electrophysiological terms, 0.138 THz waves improved synaptic transmission efficiency in the hippocampal CA1 region, and this was a slow and continuous process. This is consistent with the morphological results. This phenomenon can continue for more than 10 min after terahertz radiation ends, and these phenomena were associated with an increase in dendritic spine density. In summary, our study shows that 0.138 THz waves can modulate dynamic neuronal growth and synaptic transmission. Therefore, 0.138 terahertz waves may become a novel neuromodulation technique for modulating neuron structure and function.

All-atom molecular dynamics simulations of the combined effects of different phospholipids and cholesterol content on electroporation

The electroporation mechanism could be related to the composition of the plasma membrane, and the combined effect of different phospholipid molecules and cholesterol content on electroporation has rarely been studied nor conclusions drawn. In this paper, we applied all-atom molecular dynamics (MD) simulations to study the effects of phospholipids and cholesterol content on bilayer membrane electroporation. The palmitoyloleoylphosphatidylcholine (POPC) model, palmitoyloleoylphosphatidylethanolamine (POPE) model, and a 1 : 1 mixed model of POPC and POPE called PEPC, were the three basic models used. An electric field of 0.45 V nm-1 was applied to nine models, which were the three basic models, each with three different cholesterol content values of 0%, 24%, and 40%. The interfacial water molecules moved under the electric field and, once the first water bridge formed, the rest of the water molecules would dramatically flood into the membrane. The simulation showed that a rapid rise in the Z-component of the average dipole moment of the interfacial water molecules (Z-DM) indicated the occurrence of electroporation, and the same increment of Z-DM represented a similar change in the size of the water bridge. With the same cholesterol content, the formation of the first water bridge was the most rapid in the POPC model, regarding the average electroporation time (t ep), and the average t ep of the PEPC model was close to that of the POPE model. We speculate that the differences in membrane thickness and initial number of hydrogen bonds of the interfacial water molecules affect the average t ep for different membrane compositions. Our results reveal the influence of membrane composition on the electroporation mechanism at the molecular level.

Continuous wave irradiation at 0.1 terahertz facilitates transmembrane transport of small molecules

The bioeffects of terahertz (THz) radiation received growing attention because of its influence on the interactions between biomolecules. Our work aimed to investigate the effects of THz irradiation on cell membrane, especially cell membrane permeability. We found that 0.1 THz irradiation promoted the endocytosis of FM4-64-labeled cells and the inhibition of dynamin attenuated but did not fully abolish the THz promoted endocytosis. Moreover, 0.1 THz irradiation also promoted the transmembrane of the rhodamine, as well as the chemical compounds GDC0941 and H89, evidenced by the confocal microscope observation and the western blotting analysis, respectively. These findings demonstrated 0.1 THz irradiation facilitated the transmembrane transport of small molecules by promoting both the cellular endocytosis and the diffusion process. Our study provided direct evidence that THz could affect the cell membrane permeability, broadened the THz affected cellular physiological processes, and implied its potential application in regulating the cell membrane functions.

•

THz irradiation enhances endocytic activity of neuronal cells

•

THz irradiation increases the permeation of rhodamine derivatives into cells

•

THz irradiation promotes the cellular uptake of small drugs GDC0941 and H89

Electric-driven membrane poration: A rationale for water role in the kinetics of pore formation

Electroporation is a well-established technique used to stimulate cells, enhancing membrane permeability by inducing reversible membrane pores. In the absence of experimental observation of the dynamics of pore creation, molecular dynamics studies provide the molecular-level evidence that the electric field promotes pore formation. Although single steps in the pore formation process are well assessed, a kinetic model representing the mathematical description of the electroporation process, is lacking. In the present work we studied the basis of the pore formation process, providing a rationale for the definition of a first-order kinetic scheme. Here, authors propose a three-state kinetic model for the process based on the assessed mechanism of water defects intruding at the water/lipid interface, when applying electric field intensities at the edge of the linear regime. The methodology proposed is based on the use of two robust biophysical quantities analyzed for the water molecules intruding at the water/lipid interface: (i) number of hydrogen bonds; (ii) number of contacts. The final model, sustained by a robust statistical sampling, provides kinetic constants for the transitions from the intact bilayer state to the hydrophobic pore state.

On the molecular mechanisms implicated in the bipolar cancellation of membrane electroporation

Bipolar cancellation is the phenomenon in which the permeability of cell membranes subjected to high intensity short pulsed electric field (ns-μs range) is reduced or eliminated when the system is subjected to bipolar instead of monopolar pulses. Although several studies have tried to explain bipolar cancellation, the underlying mechanisms remain unclear. Very few articles study bipolar cancellation by means of molecular dynamics (MD) simulation. In this paper, we investigated the molecular mechanisms underlying the difference in electroporation induced by bipolar and monopolar picosecond electric pulses (EPs) using MD simulation. The electric field gradients and electric forces on water molecules of the two pulses were analyzed in detail for the first time. For a certain pulse width, when the field intensity is relatively small, the direction of bipolar electric force on the interfacial water molecule reverses as the bipolar EPs reverse, while the electric force on interfacial water molecules of the cathode side remains in the same direction as that of applied monopolar EPs. The bipolar electric force reversal delays the water protrusion and increases the pore formation time. Therefore, this phenomenon could correspond to bipolar cancellation. When the field intensity is relatively large, although the bipolar electric force direction still reverses, half of the total time of the monopolar EPs has no electric fields. The electric forces of monopolar no-field half-cycles are much smaller than those of the bipolar EPs. Therefore, the pore formation time of bipolar EPs reduces, and this phenomenon is called bipolar enhancement. The occurrence of bipolar cancellation or bipolar enhancement depends on conditions such as the width and intensity of the pulse.

Disassembly of microtubules by intense terahertz pulses

The biological effects of terahertz (THz) radiation have been observed across multiple levels of biological organization, however the sub-cellular mechanisms underlying the phenotypic changes remain to be elucidated. Filamentous protein complexes such as microtubules are essential cytoskeletal structures that regulate diverse biological functions, and these may be an important target for THz interactions underlying THz-induced effects observed at the cellular or tissue level. Here, we show disassembly of microtubules within minutes of exposure to extended trains of intense, picosecond-duration THz pulses. Further, the rate of disassembly depends on THz intensity and spectral content. As inhibition of microtubule dynamics is a mechanism of clinically-utilized anti-cancer agents, disruption of microtubule networks may indicate a potential therapeutic mechanism of intense THz pulses.

Effects of electromagnetic fields on neuronal ion channels: a systematic review

Many aspects of chemistry and biology are mediated by electromagnetic field (EMF) interactions. The central nervous system (CNS) is particularly sensitive to EMF stimuli. Studies have explored the direct effect of different EMFs on the electrical properties of neurons in the last two decades, particularly focusing on the role of voltage-gated ion channels (VGCs). This work aims to systematically review published evidence in the last two decades detailing the effects of EMFs on neuronal ion channels as per the PRISM guidelines. Following a predetermined exclusion and inclusion criteria, 22 papers were included after searches on three online databases. Changes in calcium homeostasis, attributable to the voltage-gated calcium channels, were found to be the most commonly reported result of EMF exposure. EMF effects on the neuronal landscape appear to be diverse and greatly dependent on parameters, such as the field's frequency, exposure time, and intrinsic properties of the irradiated tissue, such as the expression of VGCs. Here, we systematically clarify how neuronal ion channels are particularly affected and differentially modulated by EMFs at multiple levels, such as gating dynamics, ion conductance, concentration in the membrane, and gene and protein expression. Ion channels represent a major transducer for EMF-related effects on the CNS.

The interplay of excitation and electroporation in nanosecond pulse stimulation

Conventional electric stimuli of micro- and millisecond duration excite or activate cells at voltages 10-100 times below the electroporation threshold. This ratio is remarkably different for nanosecond electric pulses (nsEP), which caused excitation and activation only at or above the electroporation threshold in diverse cell lines, primary cardiomyocytes, neurons, and chromaffin cells. Depolarization to the excitation threshold often results from (or is assisted by) the loss of the resting membrane potential due to ion leaks across the membrane permeabilized by nsEP. Slow membrane resealing and the build-up of electroporation damages prevent repetitive excitation by nsEP. However, peripheral nerves and muscles are exempt from this rule and withstand multiple cycles of excitation by nsEP without the loss of function or signs of electroporation. We show that the damage-free excitation by nsEP may be enabled by the membrane charging time constant sufficiently large to (1) cap the peak transmembrane voltage during nsEP below the electroporation threshold, and (2) extend the post-nsEP depolarization long enough to activate voltage-gated ion channels. The low excitatory efficacy of nsEP compared to longer pulses makes them advantageous for medical applications where the neuromuscular excitation is an unwanted side effect, such as electroporation-based cancer and tissue ablation.

The Effect of KcsA Channel on Lipid Bilayer Electroporation Induced by Picosecond Pulse Trains

Membrane proteins are the major component of plasma membranes, and they play crucial roles in all organisms. To understand the influence of the presence of KcsA channel on cell membrane electroporation induced by picosecond pulse trains (psPT), in this paper, the electroporation of KcsA membrane protein system and bare lipid bilayer system (POPC) with the applied psPT are simulated using molecular dynamics (MD) method. First, we find that the average pore formation time of the KcsA system is longer than the bare system with the applied psPT. In the KcsA system, water protrusions appear more slowly. Then, the system size effects of psPT in the MD simulations are investigated. When the system size decreases, the average pore formation time of small KcsA membrane protein system is shorter than the bare system with the applied psPT. It is found that the psPT makes the protein fluctuation of small system increase greatly; meanwhile the instability of protein disturbs the water and then affects the water protrusion appearance time. Furthermore, it shows that the protein fluctuation of constant electric field is smaller than that of psPT and no field, and protein fluctuation increases with the psPT repetition frequency increasing.

Interpretation of the molecular mechanism of the electroporation induced by symmetrical bipolar picosecond pulse trains

Picosecond pulse trains (psPTs) are emerging as a new characteristic diagnostic and therapeutic tool in biomedical fields. To specifically determine the stimulus provided to cells, in this article, we use a molecular dynamics (MD) model to show the molecular mechanisms of electroporation induced by symmetrical bipolar psPTs and predict a bipolar cancellation for the studied picosecond pulses. Electric field conditions that do not cause electroporation reveal that the interfacial water molecules continuously flip and redirect as the applied bipolar psPT reverses, and the molecules cannot keep moving in one direction or leave the lipid-water interface. Based on our simulation results, we determine the threshold for electroporation with symmetrical bipolar psPTs. For a fixed electric field intensity, a lower repetition frequency leads to more rapid electroporation. For a fixed repetition frequency, a higher electric field intensity leads to more rapid electroporation. We found that the water dipole relaxation time decreases as the electric field magnitude increases. Additionally, the influences of the symmetrical bipolar psPT intensity and frequency on the pore formation time are presented. Discrete nanoscale pores can form with the applied psPT at terahertz (THz) repetition frequency. When the psPT amplitude increases or the frequency decreases, the number of water bridges will increase. Moreover, for the first time, the molecular mechanism of bipolar cancellation for the studied picosecond pulse is discussed preliminarily. Our results indicate that the influence of the unipolar picosecond pulse on the interfacial water dipoles will accumulate in one direction, but the bipolar picosecond pulse does not cause this effect.

A primary model of THz and far-infrared signal generation and conduction in neuron systems based on the hypothesis of the ordered phase of water molecules on the neuron surface I: signal characteristics

In this paper, we use the theory of quantum optics and electrodynamics to study the electromagnetic field problem in the nervous system based on the assumption of an ordered arrangement of water molecules on the neuronal surface. Using the Lagrangian of the water molecule-field ion, the dynamic equations for neural signal generation and transmission are derived. Perturbation theory and the numerical method are used to solve the dynamic equations, and the characteristics of high-frequency signals (the dispersion relation, the time domain of the field, the frequency domain waveform, etc.) are discussed. This model predicts some intrinsic vibration modes of electromagnetic radiation on the neuronal surface. The frequency range of these vibration modes is in the THz and far-infrared ranges.

Terahertz Electric Field-Induced Membrane Electroporation by Molecular Dynamics Simulations

In this paper, the membrane electroporation induced by the terahertz electric field is simulated by means of the molecular dynamics method. The influences of the waveform and frequency of the applied terahertz electric field on the electroporation and the unique features of the process of the electroporation with the applied terahertz electric field are given. It shows that whether the electroporation can happen depends on the waveform of the applied terahertz electric field when the magnitude is not large enough. No pore appears if the terahertz electric field direction periodically reverses, and dipole moments of the interfacial water and the bulk water keep reversing. The nm-scale single pore forms with the applied terahertz trapezoidal electric field. It is found that the average pore formation time is strongly influenced by the terahertz electric field frequency. An abnormal variation region that shows decline exists on the correlation curve of the average pore formation time and the trapezoidal electric field frequency, whereas the overall trend of the curve is increasing. The decrease of the water oriented polarization degree results in the increase of the electroporation time, and the abnormal variation region appearance may be related to the drastic change of average water hydrogen bond number that is resulted from the resonance of water hydrogen bond network and the applied electric field. Compared to the nanosecond electric pulse and constant electric field, the numbers of the water protrusions and the water bridges are smaller and the pore formation time is relatively longer with the applied terahertz electric field.

Monopole patch antenna for in vivo exposure to nanosecond pulsed electric fields

To explore the promising therapeutic applications of short nanosecond electric pulses, in vitro and in vivo experiments are highly required. In this paper, an exposure system based on monopole patch antenna is reported to perform in vivo experiments on newborn mice with both monopolar and bipolar nanosecond signals. Analytical design and numerical simulations of the antenna in air were carried out as well as experimental characterizations in term of scattering parameter (S ₁₁) and spatial electric field distribution. Numerical dosimetry of the setup with four newborn mice properly placed in proximity of the antenna patch was carried out, exploiting a matching technique to decrease the reflections due to dielectric discontinuities (i.e., from air to mouse tissues). Such technique consists in the use of a matching dielectric box with dielectric permittivity similar to those of the mice. The average computed electric field inside single mice was homogeneous (better than 68 %) with an efficiency higher than 20 V m^-1 V^-1 for the four exposed mice. These results demonstrate the possibility of a multiple (four) exposure of small animals to short nanosecond pulses (both monopolar and bipolar) in a controlled and efficient way.

Electroporation by subnanosecond pulses

Electropermeabilization of cell membranes by micro- and nanosecond-duration stimuli has been studied extensively, whereas effects of picosecond electric pulses (psEP) remain essentially unexplored. We utilized whole-cell patch clamp and Di-8-ANEPPS voltage-sensitive dye measurements to characterize plasma membrane effects of 500 ps stimuli in rat hippocampal neurons (RHN), NG108, and CHO cells. Even a single 500-ps pulse at 190 kV/cm increased membrane conductance and depolarized cells. These effects were augmented by applying brief psEP bursts (5-125 pulses), whereas the rate of pulse delivery (8Hz - 1 kHz) played little role. psEP-treated cells displayed large inward current at negative membrane potentials but modest or no conductance changes at positive potentials. A 1-kHz burst of 25 pulses increased the whole-cell conductance in the range (-100) - (-60) mV to 22-26 nS in RHN and NG108 cells (from 3 and 0.7 nS, respectively), but only to 5 nS in CHO (from 0.3 nS). The conductance increase was reversible within about 2 min. Such pattern of cell permeabilization, with characteristic inward rectification and slow recovery, was similar to earlier reported effects of 60- and 600-ns pulses, pointing to the similarity of structural membrane rearrangements in spite of a different membrane charging mechanism.

A molecular insight into the electro-transfer of small molecules through electropores driven by electric fields

The transport of chemical compounds across the plasma membrane into the cell is relevant for several biological and medical applications. One of the most efficient techniques to enhance this uptake is reversible electroporation. Nevertheless, the detailed molecular mechanism of transport of chemical species (dyes, drugs, genetic materials, …) following the application of electric pulses is not yet fully elucidated. In the past decade, molecular dynamics (MD) simulations have been conducted to model the effect of pulsed electric fields on membranes, describing several aspects of this phenomenon. Here, we first present a comprehensive review of the results obtained so far modeling the electroporation of lipid membranes, then we extend these findings to study the electrotransfer across lipid bilayers subject to microsecond pulsed electric fields of Tat11, a small hydrophilic charged peptide, and of siRNA. We use in particular a MD simulation protocol that allows to characterize the transport of charged species through stable pores. Unexpectedly, our results show that for an electroporated bilayer subject to transmembrane voltages in the order of 500mV, i.e. consistent with experimental conditions, both Tat11 and siRNA can translocate through nanoelectropores within tens of ns. We discuss these results in comparison to experiments in order to rationalize the mechanism of drug uptake by cells. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Membrane pore formation in atomistic and coarse-grained simulations

Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (Lα) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the Lα phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the Lα phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.