# Modeling and Simulation of Reaction Environment in Photoredox Catalysis: A Critical Review Gabriela Xavier de Oliveira, et al.

Jan 14, 2022

## Highlights

• Photoredox catalysis is a powerful tool in the field of organic synthesis.

• Overview of applications for the synthesis of added-value chemicals using photoredox catalysis.

• Photoredox catalysis in microreactors: CFD modeling as a pivotal approach towards a better understanding of the process.

• Future perspectives for the optimization of photoredox catalysis process design by photon transport simulation and CFD.

## 1 Introduction

The future development of synthetic organic chemistry involves environmentally friendly and sustainable alternatives in terms of energy to promote chemical transformations, using “light” as an abundant and renewable energy source (Narayanam and Stephenson, 2011; Shaw et al., 2016). Photochemical reactions are based on the use of light (ultraviolet, visible or solar) to provide the activation energy to induce the synthesis of a target molecule from a simple substrate. In this way, researchers explore ways to efficiently collect and use light energy for the activation of organic molecules (Noël, 2017; Crisenza and Melchiorre, 2020).

In the last decade, photoredox catalysis has emerged as a versatile platform for catalytic methodologies for organic synthesis, using photons of visible light as a reagent without tracking. In a general sense, this approach is based on the ability of metallic complexes and organic dyes to get involved in single electron transfer (SET) processes with the organic substrate after photoexcitation with visible light (Nicholls et al., 2016). This technology has become an indispensable tool in the search for new biologically active compounds for applications in medicine, fine chemical and pharmaceutical industries, as well as in many other fields of science, materials and the environment (Aillet et al., 2013; Loubière et al., 2016). Moreover, recent progress in the use of photoredox catalysis as a new conceptual approach to mediate synthetic organic reactions using visible light has been promoted by the combined efforts of three research groups: MacMillan (Shaw et al., 2016), Stephenson (Narayanam and Stephenson, 2011), Yoon (Schultz and Yoon, 2014) and many other researchers.

GRAPHICAL ABSTRACT

However, this growing interest in photochemistry has also brought to light some old and unresolved challenges such as reproducibility, scale and efficiency (Corrigan et al., 2016). Most of these issues are associated with the complexity of photochemical processes. The scale increase is hampered by the photon transport attenuation effect (Bouguer-Lambert-Beer law). Excessive radiation can be a relevant issue in using larger reactors, as reaction times would increase substantially and would also result in the formation of undesired by-products (Cambié et al., 2016). However, these problems associated with batch photochemistry can be overcome with the use of continuous flow systems.

Continuous flow microreactors have recently emerged as promising instruments for chemical synthesis. Due to the small size of the reaction channel, the microreactors favor the heat and mass transport, allow precise control of the reaction and also allow extensive penetration of light (even in concentrated solutions) becoming quite advantageous in micro (photo)chemical reactions. Several studies report that this technology is also suitable and has many benefits for photochemical transformations, increasing interest in flow photochemistry (Coyle and Oelgemöller, 2008; Knowles et al., 2012; Jensen et al., 2014; Su et al., 2014; Loubière et al., 2016; Straathof et al., 2016; Su et al., 2016).

Despite many studies aimed at the application of microreactors for organic synthesis, this approach is mostly laboratory research (Weissman and Anderson, 2015). The experimental methods provide an overview of the process and do not allow a local analysis inside the reactor. To achieve higher levels of accuracy and more rigorous and detailed analysis, microscopic models based on transport phenomena must be employed. In this sense, computational fluid dynamics (CFD) can be highlighted as a useful simulation tool.

These methods simulate the reactive flow within the entire reactor, considering the coupling of all phenomena that occur (fluid flow, mass and heat transport, and photon transport) (Oliveira de Brito Lira et al., 2021). Modeling organic reactions can be challenging due to complex reaction profiles (Armstrong et al., 2019). Therefore, once the model is validated, relevant information regarding the reactive flow in the microchannels and results regarding the yield and selectivity of the photoredox-type reaction can be easily obtained. Furthermore, prediction of the effect of different operating conditions and evaluation of different geometric configurations of the microfluidic device can also be estimated.

Thus, this article aims to assist researchers in the fundamentals of photoredox catalysis and CFD-based technology for simulation of the reaction environment in photoredox catalysis, from the compilation of data in the literature as research gaps were identified on this theme, justifying further investigation.

## 2 Photoredox Catalysis: Fundamentals and Application in the Synthesis of High Added-Value Chemicals

Photoredox catalysis is a branch of catalysis that uses light energy to accelerate a chemical reaction through simple electron transfer events. As the same name implies, “photo” refers to light and “redox” is an expression that represents the chemical processes of reduction and oxidation.

In photoredox catalysis, the absorption of a photon by the photoredox catalyst in the visible or UV region is necessary. The catalyst excited by light allows electron transfer processes to occur with the agents, and oxidation and reduction reactions may be possible. When the photoredox catalyst is excited, an electron moves from the HOMO orbital (abbreviation of the highest occupied molecular orbital) to the LUMO orbital (abbreviation of the lowest unoccupied molecular orbital). As a result, there is a lack of an electron in HOMO and an electron available in LUMO (Anslyn and Dougherty, 2006; Romero and Nicewicz, 2016; Bonardi et al., 2018). That is why the excited photoredox catalyst is both a stronger oxidizer and a stronger reducer than its ground state. Therefore, the photocatalyst (PC) can easily react with an external oxidizer or reducer. Most of the photoredox catalytic reactions follow the catalytic cycle described in Figure 1.

FIGURE 1. Oxidative and reducing quenching cycle typical for a photoredox catalyst.

Following light absorption, the photoexcited species perform the redox reactions in two different ways, depending on the reaction conditions employed and the properties of the sensitizer (Corrigan et al., 2016). As illustrated, the excited photocatalyst (PC*) reacts first with an oxidizing agent (A1, electron acceptor), which leads to PC•+ and A1•−. Then, the PC can be regenerated by a reduction reaction with a reducing agent (D2, electron donor), featuring the oxidative quenching cycle. On the other hand, in a reductive cycle, the PC* reacts first with the reducing agent (D1) which leads to the PC•− and D1•+. Again, the PC can be regenerated from a reaction with an oxidizing agent (A2) (Koike and Akita, 2014; Romero and Nicewicz, 2016; Bonardi et al., 2018). Such a single electron transfer (SET) events have led to a better radical generation and give new life to a radical transformation in synthesis. The terms reductive and oxidative can be confusing and require explanations. Oxidative means oxidation of the photo-excited species PC* concomitant with reduction of the external electron acceptor A. Reductive refers to the reduction of photo-excited species PC*, while the external electron donor D is oxidized in the same process. In either case, by adding external redox agents, the photocatalyst can be regenerated.

Despite a large number of studies on photoredox catalysis, a very small set of photoredox catalysts are employed. They can be represented by transition metals complexes or by organic dyes capable of absorbing visible light, resulting in a kind of excited state. Among the catalyzed reactions based on transition metals, they typically employ metallic polypyridyl complexes, such as Ru (II) polypyridine complexes and cyclometalated Ir (III) complexes, offering charge transfer properties for metal binders (Koike and Akita, 2014; Wang et al., 2018; Glaser and Wenger, 2020). On the other hand, metal-free conditions have been developed using organic dyes such as Eosin, Rhodamine, rose Bengal, 9-fluorenone, xanthone, methylene blue, acridinium, etc., capable of absorbing light in or near the visible region (Nicewicz and Nguyen, 2014; Wang et al., 2018; Zhou et al., 2019).

For a compound to be efficient as a photoredox catalyst, there are some desirable characteristics (Arias-Rotondo and Mccusker, 2016; Romero and Nicewicz, 2016):

1) Strong light absorption on the spectral region of the agents of interest;

2) Stability in the solution;

3) The lifetime of the excited state long enough to react with the agents;

4) The oxidation or reduction of the photoredox catalyst must be reversible to be regenerated;

5) Reversible redox behavior so that the photoredox catalyst is regenerated (to avoid any side reaction from PC•+ and PC•−).

These catalysts were successfully implemented for the discovery of new organic reactions and the synthesis of added-value chemicals, with excellent control of selectivity and stereoregularity, creating a new approach to organic transformations.

In recent years, the field of photoredox catalysis has experienced exponential growth, resulting in the development of a wide variety of new synthesis methodologies (Shaw et al., 2016; McAtee et al., 2019). In this scenario, the successful application of processes mediated by visible light to pharmaceutical products and biologically relevant molecules is essential for its continued expansion. Here, we will discuss some applications of the utility of photoredox catalysis in the synthesis of added-value chemicals. Table 1 presents a simplified summary of works reported in the literature for photoredox catalysis in the synthesis of added-value chemicals.

TABLE 1. Summary of the application of photoredox catalysis for the synthesis of added-value chemicals.

As a result, it is evident that the photoredox catalysis has become a powerful research area for chemical synthesis, with an increasing number of publications. This technology has found applications in a variety of synthetically useful reactions of organic compounds, even using metallic complexes or organic dyes as catalysts.

## 3 Batch Versus Flow Photochemistry

Photochemical reactions occur when light provides energy to trigger a reaction (Plutschack et al., 2017). It includes electron transfer chemistry initiated by the excitation of a chromophore (photoredox catalysis), an attractive method for organic synthesis. Thus, it is evident that photochemistry, in general, is heavily dependent on the efficient irradiation of the reaction mixture. There are currently two main types of reactor systems that researchers use for photochemistry: batch and flow (Figure 2).

FIGURE 2. Schematic representation and comparisons of batch and flow systems.

Continuous flow (micro) reactors have been highlighted as a technology of choice for photochemical reactions to overcome the bottleneck of dependence on efficient irradiation of the reaction mixture (Cambié et al., 2016; Reischauer and Pieber, 2021). The main benefits of continuous-flow microreactors over batch reactors are high surface-to-volume ratio and small scale, favorable characteristics for improving mass and photon transfer (Van Gerven et al., 2007; Lira et al., 2021). In addition, other advantages are short molecular diffusion distances, laminar flow, and fast and efficient heat dissipation (Vesborg et al., 2010; Padoin et al., 2016; Heggo and Ookawara, 2017; Yusuf et al., 2018). However, to achieve higher reaction rates, the distance required for the interaction between the components must be shorter. The smaller dimensions of the channels favor the decrease of sample volumes, which leads to a reduction in the amount of reagents used for tests and analyzes (Nge et al., 2013). Therefore, recent advances in this field mean better efficiency in terms of kinetics, safety and cost (Wang et al., 2002; Plutschack et al., 2017). Because of these advantages, the synthesis of high added-value compounds in continuous flow microreactors has become a promising area in recent years (Wiles and Watts, 2012; Rashmi Pradhan et al., 2019).

Additionally, when it comes to light-dependent systems (photochemical reactions), the use of continuous microreactors leads to improved irradiation on the reaction mixture due to better light penetration throughout the reactor depth and greater homogeneity of space lighting (Matsushita et al., 2007; Krivec et al., 2013). As revealed by the well-known Beer-Lambert-Bouguer law, the radiation distribution will not be uniform in a reactor due to absorption effects (Cambié et al., 2016). This equation establishes the light absorption

$(A)$

is dependent on the molecule’s molar attenuation coefficient

$(ε)$

, their concentration

$(c)$

, and the path length of the light propagation

$(l)$

(Eq. 1).

Unfortunately for batch reactors, this compensation has serious implications, especially for sizing a photochemical reaction. The decrease in radiation intensity with distance from the light source in batch reactors becomes one of the main disadvantages for photochemistry (Knowles et al., 2012; Oelgemoeller, 2012). In addition to avoiding this disadvantage of conventional reactors, microfluidic devices are also attractive for catalytic reactions due to their large surface-to-volume ratio and the excellent control over the reaction conditions they offer. Consequently, photochemical reactions in continuous flow reactors can be substantially accelerated and scaled to larger quantities (Reischauer and Pieber, 2021).

Cambié et al. (2016) reported two important reasons for photochemistry to reach so much repercussion in academic and industrial areas. The first one is the exposure of visible photo-redox catalysis to organic synthetic chemistry. The second reason is the use of continuous flow reactors. Additionally, the same authors list nine good reasons for using photoflow: 1) better irradiation of the reaction mixture; 2) reliable increase; 3) better reaction selectivity and greater reproducibility; 4) fast mixing; 5) rapid heat exchange; 6) multiphase chemistry; 7) reaction sequences in several stages; 8) immobilized catalysts and 9) greater operational safety.

Noël et al. (2011) compared the results of the microfluidic system with those obtained in batch experiments. The authors noted that there was no loss of activity when performed inflow. In fact, in the flow reaction, they observed greater activity in very short residence times. This result demonstrated one of the main advantages of flow chemistry. The authors stated that the increase in heat and mass transfer allowed the microreactor to heat more quickly and efficiently. Therefore, these properties make it ideal for quick reactions. Similarly, Nair et al. (2019) demonstrated that continuous-flow reactors have a competent configuration compared to batch reactors due to the high conversion rates and selectivity for the selective organic synthesis of aromatic alcohols. In addition, Abdiaj and Alcázar (Abdiaj and Alcázar, 2017) demonstrated an improved efficacy for photoredox and nickel catalysis for C(sp2)-C(sp3) cross-coupling using a continuous-flow reactor instead of a batch system. The authors mentioned that the space-time yield was 430 times better than the batch method, thus resulting in a more suitable procedure for drug discovery in terms of productivity per unit of time. Interestingly, Oelgemöller and Shvydkiv (Oelgemöller and Shvydkiv, 2011) summarized advances in microflow photochemical transformations. The authors presented comparative details between the batch system vs. microreactor for several examples of homogeneous and heterogeneous reactions, including photoadditions, photorearrangements, photoreductions, photodecarboxylations, photooxygenations and photochlorinations. In general, microfluidic devices have resulted in reduced irradiation times, greater selectivity and greater light efficiency. Therefore, the microreactors demonstrated their superiority over conventional techniques.

Although, in theory, the use of flow reactors for photochemical reactions has numerous benefits, as shown in this section, compared to batch reactors; however, developing a flow process can be time-consuming and there are still several limitations. One of the biggest limitations of flow chemistry is the handling of solids (Gisbertz and Pieber, 2020). Plutschack et al. (2017) developed a very interesting decision diagram to try to facilitate decision-making for flow versus batch analysis with the aim of finding the balance between convenience and achieving the overall goal. Therefore, it is expected that new equipment concepts applied to organic synthesis are expected to emerge as an important future area of research and development (R&D) and production methodology to improve synthetic efficiency through process automation and optimization.

## 4 Modeling and Simulation of Photoredox Catalytic Processes in Flow Reactor

As technologies become more developed and commercialized, they can help and speed up research based on experimental laboratory tests (Plutschack et al., 2017). Simulation and numerical modeling have become a more widely used tool to gain deeper insights into the system and to examine features that are difficult in experiments (Dong et al., 2021). Computational methods can be employed in the prediction and design of photochemical reactions, and also it allows scientists to gain deeper physical insights that are not easily accessible experimentally (Lira et al., 2021). Thus, computational fluid dynamics (CFD) has been widely and consistently used in reactive systems modeling. This development was made possible by significant advances in computing power and numerical models due to the high dimensionality, complexity, and highly non-linear dependencies on the properties and responses of continuous-flow devices. This area requires a greater computational effort, but even so, it shows promise in reducing the time of analysis and optimization of the photochemical reaction.

Photoredox catalytic processes in continuous-flow reactors are significantly more technologically complex, i.e., not as trivial as traditional batch techniques, in which the reactants are simply added into a vessel and stirred at a defined temperature until that the limiting reagent is consumed (Plutschack et al., 2017). Therefore, the limitations of photon and mass transfer can be addressed through the development of mathematical models that characterize the effect of fluid dynamics on photochemical reactions and represent the transport of light through the flow reactor. Thus, a set of partial nonlinear differential equations must be solved for a complete description of coupling, the fluid flow, mass transport, and photon transport (irradiation distribution) for photoredox catalytic processes in continuous flow reactors (Lira et al., 2021). The three-dimensional governing equations for continuity, momentum, and mass conservation of chemical species in Cartesian coordinates can be given by Eqs 24, respectively, according to Bird et al. (2007).

$∂ρ∂t+∇⋅(ρv)=0,(2)$
$ρDDtv=−∇p+μ∇2v+ρg,(3)$
$∂Ci∂t=−∇⋅(vCi)+∇⋅(Di∇Ci)+ri({Ci}),(4)$

where,

$ρ$

is the density of the mixture,

$v$

is the velocity vector,

$p$

is the pressure,

$μ$

is the dynamic viscosity of the mixture,

$ρg$

is the gravitational force,

$Ci$

is the concentration of species i,

$Di$

is the diffusion coefficient of the species i and

$ri({Ci})$

term represents the consumption or production rate of species i.

However, the main challenge for reaction engineering in photoredox catalysis processes is that, in addition to all the parameters that govern a catalytic reaction, the photocatalyst must first be activated by light. The photons distribution becomes more problematic with increasing reactor dimensions because they have a non-homogeneous light distribution. This scenario presents a non-trivial problem in modeling and optimizing the reaction parameters since the reaction rate must take into account the local volumetric rate of photon absorption (Bloh and Marschall, 2017), as the problems associated with the radioactive and chemical reaction are well interconnected. In a photoreaction, these two problems are coupled and the solution of the governing equations (mass balance of reactive species and radiative transfer) ends up being a very complex problem (Martin et al., 1999). Besides that, in the last decades, the study of photochemical processes has mainly involved research for the formation of new products, reaction kinetics and mechanisms; however, the implementation of the energy source is often neglected (Su et al., 2014). This is the main reason why very little is known today about how photon transport influences photoredox catalysis and opens new horizons for the development of studies, modeling, simulation and optimization for this type of reaction.

As mentioned earlier, one of the greatest challenges related to continuous-flow photochemistry is the participation of radiation in chemical reactions. Thus, the conversion and reaction rates are critically dependent on the efficiency of the radiation field. The lack of homogeneity of the irradiation and the corresponding variations in the reaction conditions can reduce the yield and selectivity. Therefore, in this section, we will highlight the radiation transport modeling, and later, we will list the main types of light sources that can be used in photochemical processes.

Currently, many studies report photochemical processes, either for the synthesis of added-value chemicals or for photodegradation; however, few of them are the studies that evaluate the effect of photon transport (Visan et al., 2014; Van Gerven et al., 2007). This can be solved by associating simulation software with experimental tests of chemical transformations in microreactors. Rigorous modeling of the radiation field in photochemical processes requires solving the radiation transport equation (RTE) in the appropriate coordinate system. This equation can be associated with the reactive flow inside microfluidic devices from CFD codes. The general form of the RTE can be written according to Eq. 5, and it describes the movement of the beam of light rays traveling in the

$s$

direction and having a solid angle

$Ω$

, intensity

$I$

, and wavelength

$λ$

where,

$κa,λ$

means the volumetric absorption coefficient,

$σs,λ$

is the volumetric scattering coefficient, and

$p(Ω′→Ω)$

is the phase function that describes the directional distribution of scattered radiation.

As RTE is an integral-differential equation, its analytical solutions are only possible under specific conditions; for most practical reactors, numerical methods should be used (Boyjoo et al., 2013). Different approaches for the RTE numerical solution have been developed and used in recent years to describe the radiant energy field in photochemistry processes. Examples are analytical approaches such as the Two- and Six-Flux Model (TFM and SFM, respectively) (Li Puma and Brucato, 2007; Ramírez-Cabrera et al., 2017), Monte Carlo (MC) stochastic method (Zazueta et al., 2013; Akach and Ochieng, 2018), discrete ordinates method (DOM) (Boyjoo et al., 2014), among others.

The two-flux model (TFM) and the six-flux model (SFM) were developed for analytical approximations to the RTE solution in an annular slurry photoreactor (Li Puma et al., 2004; Li Puma and Brucato, 2007). The TFM model assumes that the photons travel in two directions (forward and backwards) and the SFM model represents the external dispersion in the six directions of the Cartesian coordinates. The advantage of these models is their ability to estimate the local volumetric rate photon absorption

$(LVRPA)$

using simple algebraic equations (Peralta Muniz Moreira and Li Puma, 2021). Different configurations were also used to describe the radiation field using mathematical models, such as multi-annular reactor (Imoberdorf et al., 2007) and corrugated plate reactor (Passalía et al., 2011). However, no work was found using these radiation field models for application in micro (photo) reactors.

The P1 approximation to solving the RTE, also called the differential approximation, is the most general case for solving the radioactive problem (Arancibia-Bulnes and Cuevas, 2004; Cuevas et al., 2007; Boyjoo et al., 2013). This model assumes that the intensity field depends only weakly on the direction. A linear dependence is assumed, which transforms the RTE into a second degree linear partial differential equation (Ramírez-Cabrera et al., 2017). The P1 model can be solved for simpler cases, obtaining satisfactory results with low absorption coefficients and zero-dispersion; however, it has also been shown that the same model can be very imprecise for high absorption coefficients (Yu et al., 2008).

Another fact that requires further consideration is that photons are not entirely directed to the reaction medium after impinging the system, which directly contributes to the system′s irradiation efficiency. In this sense, Monte Carlo method can be used to solve the RTE solution, probabilistically simulating the possible trajectory of photons in a medium, including absorption, scattering and reflection (Boyjoo et al., 2013; Asadollahfardi et al., 2018; Ayoub, 2020; Akach et al., 2020). This method involves statistical sampling to model complex chemical and physical phenomena and has shown good results when applied to model the radiation field in photochemistry systems (Singh et al., 2007). Cambié et al. (2017) used Monte Carlo ray-tracing simulations for a detailed description of the photon path within the LSC-PM microdevice [Figure 3 (i-a)]. The relative contributions of the different photon fates were obtained, generating spectral distribution of the photons per each fate, resulting in a complete photon-balance graph [Figure 3 (i-b and c)]. It was found that the largest fraction in the photon balance, 43% of the photons reaching the device, is constituted by the transmitted photons. These simulation results were then used to investigate the expected efficiency of the extended devices. The authors noted that it would be possible to obtain accurate information about the instantaneous flow of photons reaching the reaction channel just by monitoring the variation in the photons emitted by the edges. Matiazzo et al. (2022) explored the impact of the LEDs-to-LEDs distance and also the evaluation of the LEDs-reactor distance from a non-sequential ray tracer based on Monte Carlo. Therefore, it is worth noting that the great advantage of algorithms based on Monte Carlo is the ability to deal with complex geometries without excessive computational demands.

FIGURE 3. (i) (a) LSC-PM device, (b) and (c) relative contribution of the different photon fates to the overall photon balance [reproduced from (Cambié et al., 2017)] and (ii) incident radiation flux under FSI and BSI mechanisms [reproduced from (Passalía et al., 2020)].

The discrete ordinates method (DOM) consists of transforming the RTE differential-integral form into a system of algebraic equations that can be solved by a computer (Jensen et al., 2007). This system of equations describes the transport of photons in a way that can be solved following the direction of propagation, starting from the values provided by the boundary conditions (Alvarado-Rolon et al., 2018). This method solves the RTE for a finite number of discrete solid angles, each associated with a vector direction (Casado et al., 2017a). Currently, the DOM model is used by CFD software to model the radiation field and can be coupled with the reaction kinetics. The ability of this model to provide more accurate results for problems with a radiation field coupled with a chemical reaction has made it widely applied for photochemical application (Coelho, 2014; Alvarado-Rolon et al., 2018; Moreno et al., 2019; Lira et al., 2020; Peralta Muniz Moreira and Li Puma, 2021). For example, Passalía et al. (2020) simulated an array of 18 LEDs to light a milli-reactor known as NETmix employing a DOM method, and they used the oxidation of trivalent arsenate to pentavalent arsenate in aqueous solutions as a reaction model [Figure 3 (ii)]. Two lighting settings were compared: front-side (FSI) and back-side (BSI) illumination. And therefore, from the radiation modeling, the results showed that although the FSI configuration presents a larger surface area, this configuration is worse or less illuminated since most of the catalytic surface corresponds to the side walls. Similarly, Casado et al. (2017b) used DOM to solve the RTE, allowing the evaluation of the radiation field at any point inside the reactor to evaluate and optimize the irradiation distribution with LEDs.

##### 4.1.1 Types of Radiation Sources

One of the most important factors in the study of photochemical processes is the choice of the light source. Some considerations must be made while choosing a light source. The optimal wavelength for the process must be determined, and this stage includes in most of the times a UV/vis characterization of the reaction substrate or photosensitizer (Williams and Kappe, 2020). The spectra distribution of some of the most used light sources can be seen in Figure 4.

FIGURE 4. Spectra of the different light-emitting sources. (A) High pressure mercury lamp, (B) sunlight, (C) monochromatic LEDs, and (D) white LED (blue LED with phosphor).

Another aspect that needs attention is the overlap between the emission spectrum of the light source and the absorption spectrum of absorbent species or photoactive catalysts. Moreover, for photochemical synthesis, the choice of a suitable light source depends on the volume of synthesis and the intensity of light required. However, for safety reasons, protection needs to be considered to prevent damage from the light source and the ozone generated if UV light is used (Liu, 2017). Lamps for photochemical applications are characterized by their spectral purity and radiant emission power. The physical principle behind all radiation sources is the luminescence of excited atoms or molecules, in which the electronic transition from excited states to ground state leads to the emission of radiation (Pareek et al., 2008). Based on this, light sources can be divided into the following categories: mercury discharge lamps, light-emitting diodes (LEDs) and sunlight.

• Mercury discharge lamps

Currently, the most popular light sources are mercury lamps due to their convenience because they are easy to use and provide light with wavelengths ranging from ultraviolet to visible (200–750 nm) (Figure 4A). These light sources are usually located outside the microreactors, and the capillary version is usually wrapped around the lamps (Su et al., 2014). The incorporation of these light sources inside the microreactors is not possible since the dimensions of these lamps are much larger than the microreactors.

In mercury gas discharge lamps, the light originates from an excited electron in the mercury atom, returning to its ground state energy level and emitting radiation (Kuipers, 2014). This light source has a long life, has a relatively efficient light emission and can emit over a wide wavelength range. Based on the mercury vapor pressure in the lamps, they can be classified as low pressure (0.6665–13.33 Pa), medium pressure (1.013 × 104–1.013 × 105 Pa) and high pressure (2.026 × 107 Pa) (Liu, 2017). Regarding use, medium and low-pressure lamps are often used for medium and low UV wavelengths (UV-A and UV-B, 350–250 nm), while high lamps emit a length from where up to 600 nm (Cambié et al., 2016).

Despite the advantages associated with this type of lamp, they have some disadvantages, of which one of the most pressing is the use of mercury. The use of mercury is severely restricted due to its high toxicity to humans and the environment.

• Solar light

The use of sunlight as a light source for chemical transformations is motivated, in part, by the perception that sunlight is effectively a source of clean and inexhaustible energy (Hoffmann, 2008). Sunlight and other polychromatic light sources are used due to their universality, as it is possible to choose wavelengths suitable for the study process (Figure 4B). However, in photochemical transformations, most organic molecules absorb only ultraviolet (UV) light, and they cannot be activated by the visible wavelengths that comprise most of the solar energy that reaches the Earth′s surface (Schultz and Yoon, 2014). Incompatible light energy can allow unwanted reaction pathways, leading to substantial by-product formation. Another disadvantage of this type of light source is attributed to the low and transient intensity of solar radiation and its diffuse nature (dispersion by clouds, buildings, etc.), significantly reducing the performance of photochemical processes (Cambié et al., 2016).

• Light-emitting diodes (LEDs)

Recently, light-emitting diodes (LEDs) have become more popular for microfluidic applications; however their narrow wavelength range may limit their use (±20 nm) (Su et al., 2014) (Figures 4C,D). The LEDs are semiconductors (PN junction) capable of emitting light when activated, due to the electroluminescence produced via recombination of electrons and gaps. The wavelength emitted by this type of light source can be adapted (from UV to Vis light) and combined with the photochemical application (Cambié et al., 2016). With LEDs, it is possible to produce artificial light with selected peak wavelength emission that closely matches the absorption peak of a known important photoreceptor (Dutta Gupta et al., 2017).

The advantages of LEDs are high energy efficiency, small size, low cost, durability, long life and can be made of non-toxic materials, different from mercury discharge lamps. With these advantages, this type of light source is ideal for combining with continuous flow reactors. Additionally, the LEDs have low energy loss, and only a certain amount of heat is generated, i.e., practically all electrical energy can be transformed into light monochrome energy (Preto et al., 2018). A disadvantage of LEDs is their light emission as a point source and the intensity of a single LED is quite low; however, this problem can be reversed by combining several LEDs to form a matrix, thus providing efficient irradiation (Sender–E et al., 2017).

### 4.2 Computational Fluid Dynamics Towards an Optimized Photoredox Catalysis Process Design

As progress continues for the photoredox catalysis, there is a constant need for a deeper understanding of how the process for refining chemical methods takes place to enable improvements in this technology. This is due to the fact that a photoredox catalysis is mainly an experimental approach, and, therefore, the method only offers an overview and does not allow a detailed study inside the reactor. In particular, a complete understanding of hydrodynamics, chemical reaction, and the influence of light in the photoredox catalysis process would be of paramount importance for the design of microreactors and optimization of the added-value chemicals synthesis. Therefore, the numerical modeling of these processes can be utilized in conjunction with the experimental component to assess the viability of the process and reach higher degrees of precision. Thus, computational simulation (CFD) becomes essential for the design and optimization of microfluidic devices applied to the photoredox catalysis, as it avoids the costly time and money expenses, in addition to providing increasingly efficient microdevices (Padoin et al., 2014; Padoin and Soares, 2017; Lira et al., 2019; Lira et al., 2020; Cao et al., 2022).

CFD stands out as a useful tool as it considers the coupling of all phenomena that occur (fluid flow, mass transport, and photon transport). Therefore, once the CFD-based model has been validated, investigations of velocity profiles, radiation and concentration distribution, prediction of the effect of different operating conditions, and geometric configurations can be easily performed with virtual prototypes (Mohajerani et al., 2012). Currently, several commercials software are used to solve the numerical equations that describe the microreactor systems. Examples of such software are ANSYS Fluent®, COMSOL Multiphysics®, OpenFOAM®, among others. Despite the vast number of software available on the market, the basic principle is the same for everyone, where the CFD code is divided into three main steps: pre-processing, solver and post-processing. With this, it is possible to evaluate fluid flows, evaluate the influence of channel geometry only with virtual prototypes and optimize the micro-reactive system (Santana et al., 2020).

It is worth noting that many studies have been carried out to simulate CFD-based micro-flow devices. Environmental protection, power plants, and added-value chemicals can be listed among the various sectors in which this technology has already been implemented, performing single-phase or multiphase flows (Cao et al., 2022). An et al. (An et al., 2012) examined eight different configurations of a microreactor in the form of a square cross-section channel via computational fluid dynamics to determine which results in the best performance. Similarly, Odiba et al. (Odiba and Olea, 2016) explored the behavior of the fluid in the microchannels of four different geometries of a microreactor from CFD simulation applied in the catalytic oxidation of volatile organic compounds. Ortega-Casanova (Ortega-Casanova, 2017) applied CFD in the optimization process of a chemical microreactor for the mixing efficiency evaluation. On the other hand, Alfano and Cassano (Alfano and Cassano, 2009) developed a general procedure to increase the scale of photoreactors of any type using modeling with the mass balances, at the moment, of radiation with the kinetic model.

However, when it comes to photoredox catalysis, few works have been reported in the literature involving CFD simulations. de Oliveira et al. (2020) implemented a CFD-based virtual model of the luminescent solar concentrator-based photomicroreactor (LSC-PM). The CFD model was validated with the experimental data based on the [4 + 2] cycloaddition reaction of 9,10-diphenylanthracene, and in the CFD code, a feedforward control system was implemented [Figure 5 (i-a and b)]. Furthermore, an investigation of the geometric characteristic of the microchannel was carried out, highlighting the ability to use the computational model for design purposes [Figure 5 (i-c)]. The formation of dead zones in the square-channel microreactor was verified when analyzing the flow effect for the analyzed geometries. It corroborated with the decrease in conversion for this channel setting compared to the curved channel. Debrouwer et al. (2020) used CFD simulation to support the beneficial effect of pulsation on the residence time distribution for MacMillan′s light-driven C(sp2)–C(sp3) cross-electrophile coupling [Figure 5 (ii-a)]. Moreover, by calculating tracer transport over time, these simulations allow comparison of vertical mixing times in the absence and presence of pulsation based on the concentration field in the reactor [Figure 5 (ii-b and c)]. It was found that the time required for the species to be transported to half the channel height is less than 1 s when the pulse is applied. Aillet et al. (2015) used the CFD code to predict the performance (conversion and photonic efficiency) of an

photochemical reaction, where A and B compete for the absorption of incident photons. Interestingly, the model of this work (Aillet et al., 2015) is expressed using classical dimensionless numbers such as Damköhler numbers I and II, absorbance, and competitive absorption factor. Differently, although Sezen-Edmonds et al. (2020) did not use CFD, they demonstrated a modeling tool to predict the performance of a photochemical reaction as a function of easily measurable parameters. Experimental data were combined with computational modeling for two case studies: 1) [2 + 2] photocycloaddition reaction and 2) decarboxylative intramolecular arene alkylation mediated by photoredox catalysis.

FIGURE 5. (i) (a) [4 + 2] cycloaddition of 9,10-diphenylanthracene (DPA) using photogenerated singlet oxygen, (b) process with and without the control system during the power disturbances, and (c) non-ideal flow effects on velocity in three different channel configurations [reproduced from (de Oliveira et al., 2020)]. (ii) (a) Ir/Ni Photoredox C(sp2)–C(sp3) Cross-Electrophile Coupling of Methyl 4-Bromobenzoate and 4-Bromotetrahydropyran in the continuous flow regime, concentration field in the reactor t = 1.0 s after tracer injection with (b) no pulsation and (c) with pulsation frequency of 3 Hz and amplitude of 33% [reproduced from (Debrouwer et al., 2020)].

It is worth noting that, in addition to CFD-based simulations, some authors are using ray-tracing modeling to visualize photon transport to optimize the design of reactors applied to photochemical reactions (Buglioni et al., 2021). It can be realized using a Monte Carlo ray-trace algorithm from the code development (Zhu et al., 2020) or through the use of commercial software such as LightTools® (Vishwanathan et al., 2015), Ray Optics Module (an add-on to COMSOL Multiphysics®) (Jacobs et al., 2020; Zheng et al., 2020), ANSYS SPEOS® (Matiazzo et al., 2022), among others.

Thus, the use of computational software, whether CFD or based on a ray-tracing algorithm, is essential for the development and improvement of reactors for photochemical reactions where there is a coupling of hydrodynamics, mass transport, and photon transport. However, more studies still need to be conducted and published for photoredox catalysis reactions so that researchers understand more deeply the phenomena that occur in this type of reaction. Therefore, the modeling and simulation of reaction environments in photoredox catalysis still has a lot to be explored to establish itself as a reliable technology.

## 5 Discussions on Future Directions and Concluding Remarks

Photoredox catalysis has been a transforming force in synthetic organic chemistry in the last decade, mainly for the pharmaceutical sector. This means that academic discoveries have a better chance of being promptly converted into future technology. There are no doubt that the photoredox catalysis processes will continue to use (micro) flow reactors to perform organic photochemical reactions. This growing popularity of the microreactors can be attributed to the high yields (in terms of space/time), high selectivity, and few side reactions (when compared to batch reactors) achieved by the majority of these reactors.

However, as the field of photoredox catalysis advances, a better understanding of how the process of refining chemical techniques works is required to enable advancements in this technology. For the design of microreactors and optimization of the synthesis of added-value compounds, a comprehensive understanding of the principles between the reagents, as well as a perception of the influence of light on the photoredox catalysis process, would be crucial. Nevertheless, the challenge is to integrate the complexity of continuous flow photochemical processes in order to elucidate some gaps found when approaching this subject.

For that, computational tools can provide, in the research and development of microfluidic devices for photoredox catalysis, an immediate qualitative and quantitative analysis of the species′ behavior and light distribution. Furthermore, knowledge of fluid flows is also essential to adequate the design and optimization of micro(photo)reactors. The great advantages of following the computational path are the decrease in the time-consuming and money expenses and the possibility of providing increasingly efficient models. In this sense, CFD can stand out as a powerful tool for solving and analyzing transport equations. However, difficulties linked to light efficiency, such as dissipation, penetration, and transmission of light, are still not taken into account in those CFD simulations. Other computational methods capable of simulating this phenomenon are being sought as a solution to this problem. Therefore, the framework consisting of the combination of photon fate modeling and CFD can be crucial to improve the performance of the photoredox catalysis process.

In this present literature review, we summarized the advances made in photoredox catalysis, highlighting the computational numerical approach. It was noteworthy that great efforts have been made to improve the understanding of the fundamentals and application of this technology for the added-value chemicals synthesis. Furthermore, the advantages of working with a continuous flow system were also understood. Despite efforts, photoredox catalysis is still a challenging issue as it involves the consideration of several simultaneous phenomena. Therefore, computer simulation of these reaction systems is a scientific hotspot for researchers to clarify improves the interaction between reagents and the influence of light on the reaction medium. However, the simulation of photoredox catalysis for the added-value chemicals synthesis still has a lot to be explored to establish itself as a reliable approach.

## Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Acknowledgments

The authors acknowledge CAPES (Coordination for the Improvement of Higher-Level Personnel and CNPq (National Council for Scientific and Technological Development) for the financial support.