Formation of Metallic Ag on AgBr by Femtosecond Laser Irradiation

: Laser irradiation of materials induces changes in their structure and functional properties. In this work, lattice heating and electronic excitation on silver bromide (AgBr), provoked by femtosecond laser irradiation, have been investigated by ﬁnite-temperature density functional theory and ab initio molecular dynamics calculations by using the two-temperature model. According to our results, the electronic temperature of 0.25 eV is enough to excite the electrons from the valence to the conduction band, whereas 1.00 eV changes the structural properties of the irradiated AgBr material. Charge density simulations also show that an Ag clustering process and the formation of Br 3 − complexes take place when the electronic temperature reaches 2.00 eV and 5.00 eV, respectively. The present results can be used to obtain coherent control of the extreme nonequilibrium conditions due to femtosecond laser irradiation for designing new functional materials.


Introduction
Femtosecond laser irradiation (FLI) on semiconductor structures has attracted the attention of the scientific community due to its ability to produce and functionalize new materials by engineering their structure, and the electronic and optical properties [1][2][3][4], with applications in medicine and the electronic industry [5][6][7]. FLI induces structural changes and electronic excitation on the irradiated material [8]. The high kinetic energy of the excited atoms leads to the loss of crystallinity and the irradiated material tends to transform into amorphous [9]. Amorphous regions present a high degree of structural disorganization. Under the effect of FLI, both reduction and oxidation processes are expected, forming metallic clusters constituted by reduced species in specific lattice regions [10]. However, understanding how these chemical and physical processes occur and their nature remains a challenge.
Previously, FLI-based techniques have been employed in different applications, such as nonlinear lithography [11], the crystallization of amorphous semiconductors [12], and the increase in the performance of microsuper capacitors [13], as well as to improve the functionalities of existing materials [14] and to perform intelligent real-time control of pulses [15]. These experimental results boost theoretical studies based on real-time ab initio calculations [16]. Due to mass and charge differences between the electrons and nuclei, the FLI first provides a faster electronic response and an immediate electron excitation. As the energy of the laser pulses is higher than the band gap of materials, the closest electrons from the Fermi energy populate the empty levels of the conduction band (CB). Increasing the electronic temperature, electrons from the inner levels of the valence band (VB) will also

Theoretical Approach and Computational Details
FT-DFT and AIMD calculations were carried out using the generalized gradient approximation with the exchange correlation energy functional as proposed by Perdew-Burke-Erzenhof (PBE) [32]. The projector augmented-wave (PAW) method was used to describe the electron-ion interactions [33,34], as implemented in the Vienna ab initio simulation package (VASP) [35], version 6.2.1, where the following PAW valence states 4d 10 , 5s 1 and 4s 2 , 4p 5 for Ag and Br were employed, respectively. The equilibrium geometries for AgBr were obtained by stress-tensor minimization with a plane wave expansion cutoff of 500 eV, while the force optimization cutoff was 281 eV, which is 12.5% larger than the cutoff energy recommended by the VASP package. The difference between the stress tensor and the atomic forces calculation is the number of plane waves, increasing (decreasing) the computational cost with the higher (smaller) number of plane waves used.
To project the maximally localized Wannier functions (MLWF), the Wannier90 code was used, as implemented by Marzari and Vanderbilt [36]. In these calculations, the inter-polated band structure was compared with DFT calculations. For the radial distribution function (RDF), the Pymatgen (Python Materials Genomics) code was employed [37], which is an open-source Python library for materials analysis. In our molecular dynamics investigations, the mean square displacement (MSD) was computed through the VASPKIT code [38], version 1.2.5. Additionally, the Phonopy code, version 2.11, allowed us to investigate the phonon dispersion along the high-symmetry q-points by using the atomic forces displacement and the force constants calculations [39].
To simulate the FLI, T e within the FT-DFT approach was employed using a Gaussianlike function for the Brillouin zone smearing [40].
where ε, µ, and σ = K B T e are the energy band, chemical potential, and the smearing method, written as a function of the Boltzmann constant and electronic temperature, as implemented in the VASP package. In our calculations, a low electronic temperature of T e = 0.01 eV was used to minimize the stress tensor and obtain the equilibrium volume of AgBr, and after that, the width of the Brillouin zone integration was varied from T e = 0.01 eV until 6.00 eV, fixing the AgBr lattice parameter. The TTM was adopted using AIMD simulations and the FT-DFT approach [18,[40][41][42]. In these calculations, a supercell was built from the primitive unit cell multiplied by two in the three space directions (resulting in 32 Ag and 32 Br atoms). Only a reduced gammapoint was used for the Brillouin zone integration to minimize the computational efforts. The microcanonical ensemble (NVE) was adopted, the time step in the integration of the equations of motion was 0.5 fs, and simulations spanned for 5.0 ps.
For the calculation of the structural and vibrational properties, a 6 × 6 × 6 k-mesh was used for the Brillouin zone integration and 12 × 12 × 12 for the electronic and charge calculations using a Gaussian smearing method. In addition, the electronic convergence criterion was 10 −6 eV, and the atoms were allowed to relax until all forces were smaller than 10 −2 eV·A −1 on every atom.

Results
The mechanical and electronic properties of AgBr in the zinc-blend bulk phase were investigated within the FT-DFT approach to elucidate the faster electron response to laser irradiation. Furthermore, ab initio molecular dynamics simulations were carried out to investigate the effects on the slower excited ion response and the initial diffusion processes.

FT-DFT Investigation: Limits for Electronic Excitation
To ensure that the choice of the convergence criteria to solve the Kohn-Sham scheme are satisfied for a given k-mesh, T e → 0 was adopted. Smaller electronic temperatures are generally used to avoid partial occupations of electronic states near the Fermi energy. However, within the ab initio DFT calculations, a very small T e provides the discontinuity of the Brillouin zone integration, and therefore, T e must be tested for each particular system. Here, T e = 0.01 eV was used to obtain the ground state properties.
The converged lattice parameter for AgBr was 6.30 Å, which is in good agreement with the reported value of 6.33 Å [43], with a deviation of only 0.03 Å (0.47%). The AgBr unit-cell has a cubic face-centered geometry with 4 Ag and 4 Br atoms, as represented in Figure 1. Additionally, the calculated band gap energy with PBE was 1.10 eV, also close to the experimental value of 1.17 eV [44], which is a direct band gap in the high symmetry Γ-point. These calculations show that the valence band maximum (VBM) is triply degenerated in the Γ-point corresponding to d xy , d xz , and d yz Ag levels (t 2g states) and p x , p y , and p z Br states. Γ-point. These calculations show that the valence band maximum (VBM) is triply degenerated in the Γ-point corresponding to dxy, dxz, and dyz Ag levels (t2g states) and px, py, and pz Br states. The dz2 and dx2-y2 electronic levels (eg-states) in the VBM are empty, and the conduction band minimum (CBM) only shows the nondegenerate s states. For Te = 0.01 eV, the electrons in the VBM preferably occupy the dxy-Ag state and the px and py Br states. The CBM is only composed of s-Ag and s-Br states. Figure 1 displays the agreement between the electronic band structure obtained by DFT and Wannier interpolation, the high symmetry points of the Brillouin zone calculation, the Wannier d-Ag orbitals, and the summed electronic density of states.
FT-DFT was employed to understand the electronic excitation due to femtosecond laser irradiation. In these investigations, comparative calculations increasing Te were performed to obtain the structural parameters of the AgBr lattice converged at Te = 0.01 eV. Analyzing the electron density distribution at Te = 0.25 eV shows that electrons are already in the conduction band. This phenomenon becomes more evident for Te = 0.75 eV, which varies the electronic chemical potential and, consequently, the chemical environment of the AgBr semiconductor. In addition, an increase in the electronic entropy was observed at Te = 0.75 eV, resulting in band gap shrinkage. However, for Te = 1.00 eV, AgBr presents metallic characteristics, i.e., the Fermi energy and CBM are at close energy levels, the electronic occupation increases in the CBM, and the electrons begin to occupy higher levels of the conduction band. In addition, the electronic disorder increases significantly at Te = 1.00 eV, indicating that this is a limit for the Te approach for studying AgBr with semiconducting properties. As electronic disorder represents the energy difference between the total energy obtained with Te → 0 and the free energy; an increase in the DFT-total energy was observed when increasing the electronic excitation from Te > 1.00 eV. The d z2 and d x2-y2 electronic levels (e g -states) in the VBM are empty, and the conduction band minimum (CBM) only shows the nondegenerate s states. For T e = 0.01 eV, the electrons in the VBM preferably occupy the d xy -Ag state and the p x and p y Br states. The CBM is only composed of s-Ag and s-Br states. Figure 1 displays the agreement between the electronic band structure obtained by DFT and Wannier interpolation, the high symmetry points of the Brillouin zone calculation, the Wannier d-Ag orbitals, and the summed electronic density of states.
FT-DFT was employed to understand the electronic excitation due to femtosecond laser irradiation. In these investigations, comparative calculations increasing T e were performed to obtain the structural parameters of the AgBr lattice converged at T e = 0.01 eV. Analyzing the electron density distribution at T e = 0.25 eV shows that electrons are already in the conduction band. This phenomenon becomes more evident for T e = 0.75 eV, which varies the electronic chemical potential and, consequently, the chemical environment of the AgBr semiconductor. In addition, an increase in the electronic entropy was observed at T e = 0.75 eV, resulting in band gap shrinkage. However, for T e = 1.00 eV, AgBr presents metallic characteristics, i.e., the Fermi energy and CBM are at close energy levels, the electronic occupation increases in the CBM, and the electrons begin to occupy higher levels of the conduction band. In addition, the electronic disorder increases significantly at T e = 1.00 eV, indicating that this is a limit for the T e approach for studying AgBr with semiconducting properties. As electronic disorder represents the energy difference between the total energy obtained with T e → 0 and the free energy; an increase in the DFT-total energy was observed when increasing the electronic excitation from T e > 1.00 eV.
The occupation in the CBM increases until T e = 3.00 eV and then decreases. Thus, T e = 3.00 eV is our proposed limit for the electronic excitation on AgBr within the FT-DFT due to band structure stability in the presence of hot electrons. For T e > 3.00 eV, the energy levels decrease with increasing T e . More details of the band structure obtained by DFT calculations and Wannier interpolation for different electronic temperatures are reported in Figure S1 of the Supplementary Materials (SM). Figure S2 of the Supplementary Materials represents the charge density for CBM, the next unoccupied state of the conduction band (CB), and the VBM levels as a function of the electronic temperature. For T e = 0.01 eV, the major charge contributions in the VBM are given by the Ag atoms, while the Br atoms provide the electronic density in the CBM, and the CB is composed of the contributions of Ag atoms. By increasing the electronic temperature, the charge density in the electronic states changes, and AgBr loses its semiconductor properties at T e = 1.00 eV. A charge transfer process between Ag and Br atoms at this electronic temperature is observed. For T e > 2.00 eV, the charge density of the VBM, CBM, and CB in the AgBr ground state is mainly provided by Br. Then, the crystallinity and electronic conductivity are modified by the structural and electronic charges in the electronic states provoked by FLI. Figure  S3 displays the increase (decrease) in the electron (hole) effective mass with the increase in T e , denoting the possibility of engineering the transport properties of AgBr. Similarly, Figure S4 shows the calculation of the d-Ag band center as a function of the electronic excitation as a signature of how FLI can improve the photocatalytic activity of AgBr.
The structural stability of the zinc blend AgBr phase was investigated by calculating the elastic constants and the vibrational phonon dispersion as a function of T e , as shown in Figures 2 and 3, respectively. The structural stability criteria for the cubic face-centered phase are met when the following elastic constants C 11 − C 12 > 0, C 11 + 2C 12 > 0, and C 44 > 0 [45]. The increase in the electronic temperature from 0.01 eV up to 1.00 eV does not alter the structural properties, but for T e > 1.00 eV, the elastic constants decrease smoothly until the system becomes unstable at T e = 3.00 eV. The bulk (B 0 = 1/3(C 11 + 2C 12 )) and shear (S 0 = 1/2(C 11 − C 12 )) moduli also show the limit of T e = 3.00 eV. This limit is also observed for phonon dispersion, where the phonon band structure remains practically the same in the range 0.01 eV < T e < 1.00 eV; however, for T e = 2.00 eV, some divergences appear in the acoustic modes, while for T e = 3.00 eV, the acoustic modes become imaginary, demonstrating that the structure becomes unstable.

Two Temperature Model: Ab Initio Molecular Dynamics with Excited Atoms
An analysis of the results shows that until T e = 1.00 eV, the increase in kinetic energy does not promote significant atomic displacements in the heated lattice. At higher T e , in the first stage, the strong Coulomb repulsion between excited atoms, which shows expanded orbitals, causes the atoms to oscillate around their equilibrium positions. This effect introduces noncrystalline areas and decreases the overlap region between the Ag and Br atoms. However, new energy levels may be generated within the band gap. This is the starting point for symmetry breaking, causing polarization of the AgBr structure. In this way, the system that was polarized at short and medium ranges produces a long-range polarization. This fact leads to a significant structural rearrangement, with the formation of new absorption centers of photons and the creation of metallic Ag. In the second stage, there are mobile electrons in the CB and holes in the VB. Due to the reduction of the Ag cations and the oxidation of the Br anions, AgBr, which is a semiconductor at low electronic temperatures, evolves to metallic behavior. At this stage, clusters of Ag and Br atoms are formed. Geometrical alterations in these clusters, such as distortions, breathing, and tilt, create a large variety of local structures associated with AgBr quasicrystals. In addition, an increase in the conductivity can be related to localized surface plasmon resonance effects producing metallic Ag.
At T e = 2.00 eV, the Ag atoms tend to agglomerate in certain regions of the lattice, indicating the beginning of the Ag clustering process. The simulation time for AgBr was extended to 10.0 ps for an electronic excitation of 2.0 eV and the results obtained follow the same trends presented at 5.0 ps. Higher electronic temperatures (T e > 2.00 eV) accelerate the Ag clustering process in the heated lattice, removing the quasicrystallinity previously obtained with lower electronic temperatures and turning it into an amorphous solid. In the silver diffusion process, the Ag clusters formed present different types of symmetry, which depend on the simulation time and the intensity of the electronic excitation. Due to the intensity of the laser, we can infer that the Br 2 molecules evaporate after the beginning of the Ag clustering and, after intense laser irradiation, Br 3 − species are formed [46]. The laserinduced fluorescence properties of Br 2 were investigated in Ref. [47]. Figure 4 illustrates snapshots from the simulations at different electronic temperatures, where Ag clustering regions and Br 3 − complexes are depicted at T e = 2.00 eV.
At Te = 2.00 eV, the Ag atoms tend to agglomerate in certain regions of the lattice, indicating the beginning of the Ag clustering process. The simulation time for AgBr was extended to 10.0 ps for an electronic excitation of 2.0 eV and the results obtained follow the same trends presented at 5.0 ps. Higher electronic temperatures (Te > 2.00 eV) accelerate the Ag clustering process in the heated lattice, removing the quasicrystallinity previously obtained with lower electronic temperatures and turning it into an amorphous solid. In the silver diffusion process, the Ag clusters formed present different types of symmetry, which depend on the simulation time and the intensity of the electronic excitation. Due to the intensity of the laser, we can infer that the Br2 molecules evaporate after the beginning of the Ag clustering and, after intense laser irradiation, Br3 -species are formed [46]. The laser-induced fluorescence properties of Br2 were investigated in Ref. [47]. Figure 4 illustrates snapshots from the simulations at different electronic temperatures, where Ag clustering regions and Br3 -complexes are depicted at Te = 2.00 eV.  In an analysis of the random distribution of the formed clusters of Ag and Br, the cluster sizes increase as the electronic temperature increases. Additionally, the shape, geometry, and symmetry of the formed clusters vary as a function of the irradiation time and T e . These, in turn, connect to each other to form superclusters with low symmetry. Furthermore, for T e = 5.00 eV and T e = 6.00 eV, the system exhibits a high degree of disorder typical of amorphous materials. Figure 5 displays the evolution of the values of T i as a function of time for the applied T e . The different electronic temperatures used in each simulation result in different kinetic energies for the atoms and, consequently, different lattice temperatures. At the lowest T e (1.00 eV), the lattice is lightly heated, and atoms stay at equilibrium positions. However, at 2.00 eV, the lattice temperature is still low (ca. 170 K); however, significant structural changes occur. A significant increase in the lattice temperature was observed at 3.00 eV and 4.00 eV, namely, 600 K and 1660 K, respectively.
To investigate the structural arrangements, the radial distribution functions (RDFs) were computed at different electronic temperatures after 5.00 ps of simulation, and they are depicted in Figure 6. For lower electronic temperatures up to 1.00 eV, the peaks for the Ag-Ag distances are approximately equally distributed, as expected from a crystalline structure. For T e = 2.00 eV, a single and well-defined Ag-Ag peak emerges at 2.8 Å, indicating the formation of small Ag clusters. At this temperature, the Br-Br profile loses the characteristic crystalline behavior, and the shortest Br-Br distances are approximately 3.5 Å. For higher temperatures, the distinct Ag-Ag peak is slightly shifted to 3.0 Å and increases, reflecting certain growth in the formed Ag clusters. Additionally, the appearance of Br-Br distances at 2.28 Å confirms the formation of Br 3 − species [47]. Additional simulations at T e = 2.00 eV with initial system temperatures set to 10 K, 100 K, and 300 K were conducted. Figure S5 depicts the corresponding RDFs showing similar features and demonstrating that the structural rearrangements produced at T e = 2.00 eV are a consequence of electronic excitation and are not caused by the kinetic energy associated with the lattice temperature. and Te. These, in turn, connect to each other to form superclusters with low symmetry. Furthermore, for Te = 5.00 eV and Te = 6.00 eV, the system exhibits a high degree of disorder typical of amorphous materials. Figure 5 displays the evolution of the values of Ti as a function of time for the applied Te. The different electronic temperatures used in each simulation result in different kinetic energies for the atoms and, consequently, different lattice temperatures. At the lowest Te (1.00 eV), the lattice is lightly heated, and atoms stay at equilibrium positions. However, at 2.00 eV, the lattice temperature is still low (ca. 170 K); however, significant structural changes occur. A significant increase in the lattice temperature was observed at 3.00 eV and 4.00 eV, namely, 600 K and 1660 K, respectively. To investigate the structural arrangements, the radial distribution functions (RDFs) were computed at different electronic temperatures after 5.00 ps of simulation, and they are depicted in Figure 6. For lower electronic temperatures up to 1.00 eV, the peaks for the Ag-Ag distances are approximately equally distributed, as expected from a crystalline structure. For Te = 2.00 eV, a single and well-defined Ag-Ag peak emerges at 2.8 Å, indicating the formation of small Ag clusters. At this temperature, the Br-Br profile loses the characteristic crystalline behavior, and the shortest Br-Br distances are approximately 3.5 Å. For higher temperatures, the distinct Ag-Ag peak is slightly shifted to 3.0 Å and increases, reflecting certain growth in the formed Ag clusters. Additionally, the appearance of Br-Br distances at 2.28 Å confirms the formation of Br3 -species [47]. Additional simulations at Te = 2.00 eV with initial system temperatures set to 10 K, 100 K, and 300 K were conducted. Figure S5 depicts the corresponding RDFs showing similar features and demonstrating that the structural rearrangements produced at Te = 2.00 eV are a consequence of electronic excitation and are not caused by the kinetic energy associated with the lattice temperature.  Figure 4. In the top, center and bottom panels, the radial distribution function is calculated for Ag-Ag, Br-Br and AgBr, respectively, as a function of the electronic temperature.
Furthermore, the diffusion processes induced by the laser irradiation were monitored by computing the MSD, as depicted in Figure 7. For low electronic temperatures (Te < 1.00 eV), the MSD for both Ag and Br atoms present small values, indicating that they oscillate around their equilibrium positions. However, as mentioned before, at 2.00 eV, the  Figure 4. In the top, center and bottom panels, the radial distribution function is calculated for Ag-Ag, Br-Br and AgBr, respectively, as a function of the electronic temperature. Furthermore, the diffusion processes induced by the laser irradiation were monitored by computing the MSD, as depicted in Figure 7. For low electronic temperatures (T e < 1.00 eV), the MSD for both Ag and Br atoms present small values, indicating that they oscillate around their equilibrium positions. However, as mentioned before, at 2.00 eV, the electronic excitation is enough to promote the diffusion of both atom species, with higher mobility for Ag and allowing the consequent formation of Ag clusters. For excitation at 3.00 eV, the MSD slopes are pronounced and similar for both Ag and Br, indicating the loss of crystalline structure and enhancing the growth of the Ag clusters and the formation of Br 3 − species. Figure 6. Radial distribution function for the geometries shown in Figure 4. In the top, center and bottom panels, the radial distribution function is calculated for Ag-Ag, Br-Br and AgBr, respectively, as a function of the electronic temperature.
Furthermore, the diffusion processes induced by the laser irradiation were monitored by computing the MSD, as depicted in Figure 7. For low electronic temperatures (Te < 1.00 eV), the MSD for both Ag and Br atoms present small values, indicating that they oscillate around their equilibrium positions. However, as mentioned before, at 2.00 eV, the electronic excitation is enough to promote the diffusion of both atom species, with higher mobility for Ag and allowing the consequent formation of Ag clusters. For excitation at 3.00 eV, the MSD slopes are pronounced and similar for both Ag and Br, indicating the loss of crystalline structure and enhancing the growth of the Ag clusters and the formation of Br3 -species.

Ag and Br Clustering Processes: Charge Flow
The analysis of the effective Bader charges allows us to investigate the charge transfer processes and understand the attractive or repulsive Coulomb interactions between the atoms during the diffusion mechanisms, leading to the formation of metallic clusters. Figure 8 depicts the calculated effective Bader charges for each T e . The effective Bader charge (Q eff B ) (in |e| units) is given by the difference between the number of valence electrons (Z val ) and the calculated Bader (Q Bader ) charge, as Q eff B = Z val − Q Bader .

Ag and Br Clustering Processes: Charge Flow
The analysis of the effective Bader charges allows us to investigate the charge transfer processes and understand the attractive or repulsive Coulomb interactions between the atoms during the diffusion mechanisms, leading to the formation of metallic clusters. Figure 8 depicts the calculated effective Bader charges for each Te. The effective Bader charge (Qeff B ) (in |e| units) is given by the difference between the number of valence electrons (Zval) and the calculated Bader (Q Bader ) charge, as Qeff B = Zval − Q Bader . For Te = 0.01 eV, after 5.00 ps of TTM simulations, as expected from the ionic nature of AgBr, the Ag atoms are positively charged (ca. + 0.25 |e|), whereas the Br species show negative charges. However, when the system is excited at 1.00 eV, the average charges still maintain similar behavior. At 2.00 eV, when diffusion processes of Ag atoms are For T e = 0.01 eV, after 5.00 ps of TTM simulations, as expected from the ionic nature of AgBr, the Ag atoms are positively charged (ca. + 0.25 |e|), whereas the Br species show negative charges. However, when the system is excited at 1.00 eV, the average charges still maintain similar behavior. At 2.00 eV, when diffusion processes of Ag atoms are initiated, the modulus of the average charges decreases, indicating the reduction and oxidation of Ag cations and Br anions, respectively. This process continues at the higher T e investigated up to 6.00 eV, reflecting the formation of metallic Ag clusters, Br 2 molecules, and Br 3 − species.

Conclusions
The energy provided by FLI on AgBr in the zinc blend structure changes the structure and electronic properties under extreme nonequilibrium conditions. Both the FT-DFT and AIMD methods are employed using the two-temperature model (TTM). The AgBr equilibrium bulk volume and the band structure were simulated with the smaller T e = 0.01 eV, and after that, total energy calculations were performed to study the electronic excitation. The main results of the present study can be summarized as follows: (i) the electronic properties are practically the same until T e = 0.10 eV and the electrons begin to occupy the CBM for T e = 0.25 eV; (ii) the structural and vibrational properties do not alter up to T e = 1.00 eV, while for electronic excitation T e = 3.00 eV, it was noted that the system becomes unstable; (iii) for smaller electronic temperatures than T e = 1.00 eV, the atoms oscillate near their equilibrium positions, while for higher electronic excitation, the atomic displacements increased as a function of the T e ; (iv) for T e = 2.00 eV, the kinetic energy of Ag, in comparison to Br, induces the formation process of Ag clusters and the subsequent diffusion process, while the Br atoms form Br 2 molecules and evaporate instantly; (v) the increase in T e , favors the Ag clustering process; (vi) for T e = 5.00 eV and T e = 6.00 eV, the formation of Br 3 − complexes is observed; and (vii) RDF, MSD and electronic charge investigations for each electronic temperature confirm the appearance of the Ag clustering process, and the formation of Br complexes in specific regions of the material.
We hope that the present results contribute to the understanding and rationalization of the structural and electronic effects that FLI induces on materials and can be employed as a model for the rational design of new functional materials that can be applied in the microelectronics industry.
Supplementary Materials: The following supporting information can be downloaded at https: //www.mdpi.com/article/10.3390/physchem2020013/s1, Figure S1: Electronic band structure for AgBr as a function of T e . Figure S2: Charge density for AgBr. Figure S3: Electron and hole effective mass for AgBr as a function of T e . Figure S4: d-Ag band center for AgBr. Figure S5: Radial distribution function for AgBr at 2.0 eV for 10 K, 100 K and 300 K.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.